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FOREWORD 


The  work  reported  herein  was  performed  during  the  period  from  October 
1985  to  May  1989  under  Contract  F33615-85-C-3610  from  the  Air  Force  Wright 
Research  and  Development  Center  (formerly  Air  Force  Wright  Aeronautical 
Laboratories),  Air  Force  Systems  Command.  Mr.  Thomas  Gentry  of  WRDC  was 
the  original  Contract  Technical  Monitor;  this  duty  was  later  transferred 
to  Captain  Mark  Detroit  (USAF) .  Mr.  Duane  McRuer  was  the  STI  Technical 
Director.  The  STI  Project  Engineers  were  Messrs.  Roger  H.  Hoh  and  David 
G.  Mitchell. 

The  project  was  initiated  to  explore  the  modern  nature  of  minimum 
flying  qualities  in  the  presence  of  modern  aircraft  and  multi -redundant 
flight  control  system  technology.  It  had  several  phases  including:  11  an 
intensive  effort  to  develop  and/or  elaborate  existing  pilot  modeling  anal¬ 
ysis  techniques  to  apply  to  situations  associated  with  minimum  flying 
qualities,  divided  attention  pilot  operations,  and  multi-axis  control 
tasks;  2)  preliminary  analyses  and  associated  fixed  base  simulations  to 
expand  the  meager  multi -axis  data  base  and  to  serve  as  pilot  studies  for 
more  extensive  simulations  on  the  LAMARS;  3)  an  extensive  simulation 
program  on  LAMARS  to  investigate  minimum  flying  qualities  and  related 
situations;  and  4)  analysis  and  interpretation  of  both  the  early  and 
LAMARS  simulation  efforts  in  the  context  of  the  pilot  modeling  advances. 
The  project  documentation  appears  in  three  volumes.  Volumes  II  and  III 
present  the  results  of  1)  above  while  Volume  I  covers  2)  through  4)  above. 
This  Volume  (II)  is  a  stand-alone  monograph  on  pilot  modeling,  including 
procedures  for  estimating  pilot  workload  as  "measured"  by  pilot  ratings. 
Volume  III  is  a  stand-alone  monograph  which  presents  a  detailed  implemen¬ 
tation  of  a  much  expanded  version  of  the  human  optimal  control  model  on 
Program  CC.  This  permits  detailed  analyses  using  algorithmic  pilot  models 
on  personal  computers  with  commercially  available  controls  analysis  soft¬ 
ware.  It  is  expected  to  make  pilot -vehicle  analyses  by  flying  qualities 
engineers  possible  on  a  more  routine  basis. 

The  authors  wish  to  acknowledge  the  significant  contributions  of 
Professor  David  Schmidt,  now  at  Arizona  State  University,  and  Professor 
Ronald  A.  Hess  of  University  of  California  at  Davis  for  assistance  with 

ix 


the  optimal  control  module.  Dr.  Schmidt  accomplished  an  important  break¬ 
through  in  pilot  rating  assessment  using  the  OCM  which  is  summarized  in 
the  appendix  to  this  report.  Dr.  Hess  was  of  great  help  in  increasing  our 
understanding  of  the  OCM  to  the  point  where  the  new  implementation  of 
Volume  III  could  be  undertaken.  The  authors  are  also  grateful  to  the  STI 
publications  staff,  Ms.  Bess  Shields,  Ms.  Dorie  Taylor,  Ms.  Laura 
Sherard,  and  Mr.  Charles  Reaber  for  their  careful  work  in  preparing  the 
report  for  publication. 
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SECTION  I 


INTRODUCTION 


A .  GENERAL  BACKGROUND 

The  human  operator  in  a  man-machine  system  is  the  archetype  hier¬ 
archical,  adaptive,  optimalizing,  decision-making  controller.  In  accom¬ 
plishing  these  functions  the  human  exhibits  a  bewildering  variety  of 
behavior  which  defy  quantitative  description  when  considered  in  the  large. 
Nonetheless,  since  World  War  II  scientists  and  engineers  have  attempted  to 
describe  specific  elements  within  this  functional  li%t  in  terms  of  quan¬ 
titative  models.  Because  we  are  speaking  of  human  activities  in  a  control 
context,  the  field  of  control  theory  has  been  a  major  source  of  paradigms 
for  quantitative  descriptions.  Because  control  theories  can  also  be  clas¬ 
sified  using  similar  adjectives,  it  is  not  surprising  that  almost  every 
new  advance  in  control  theory  has  led  to  attempts  to  better  understand 
additional  aspects  of  human  behavior  in  the  perspective  of  this  advance. 
When  these  attempts  have  been  fruitful  a  control  theory  paradigm  has 
emerged  which  is  useful  in  quantifying  the  human's  operations.  Just  as 
theory  has  been  used  to  "explain"  experiment,  so  unexplained  experimental 
results  beget  new  theory.  The  results  of  this  widespread  synergistic 
activity  have  been  documented  in  hundreds  of  research  papers  and  in  a 
series  of  sum  mary  surveys  which  have  appeared  aperiodically .  (A  chrono¬ 
logical  listing  of  surveys  is  given  at  the  end  of  this  report,  succeeding 
the  reference  list).  As  a  consequence,  much  of  the  successful  art  is  now 
mature.  Furthermore,  it  has  become  a  fundamental  mode  of  thinking  on  the 
part  of  technical  practitioners  in  the  fields  of  operator/vehicle  control 
system  integration,  vehicle  handling  qualities  and,  indeed,  all  aspects  of 
interactive  man-machine  systems. 

Besides  the  technological  aspects  of  manual  control,  interdiscip¬ 
linary  activities  between  control  engineers,  physiologists,  and  experimen¬ 
tal  psychologists  have  led  to  control  theory  descriptions  of  human  subsys¬ 
tem  behavior  and  to  the  interpretation  of  the  human's  psychophysiological 
outputs  in  control  engineering  terms.  These  interdisciplinary  areas  have 
been  especially  productive  in  building  psychophysiological  models  of  those 
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human  subsystems  involved  in  the  human  controller,  in  understanding  bio¬ 
dynamics  as  affected  by  environmental  variables,  in  establishing  connec¬ 
tions  between  behavior  and  subjective  workload  indices,  and  in  interpre¬ 
ting  objectively  the  effects  of  alcohol,  drugs,  fatigue,  etc.,  as  operator 
impairments . 

From  this  rich  variety  of  man-machine  control  aspects  that  have  been 
addressed,  the  emphasis  here  will  be  confined  to  models  particularly  per¬ 
tinent  to  ordinary  and  minimum  flying  qualities  situations.  Although  this 
will  not  oe  an  exhaustive  cross-section  of  the  field,  the  models  treated 
will  be  quite  comprehensive  (as  encapsulations  of  experimental  data)  and 
representative  of  useful  theory  (in  that  both  classical  and  modern  control 
viewpoints  are  presented) . 

Flying  qualities  in  general  can  be  divided  into  "unattended"  (and 
trim),  large  amplitude  maneuvering,  and  "closed- loop"  operations.  All 
three  categories  have  some  degree  of  pilot  interaction.  The  differences 
in  pilot  behavior  exhibited  in  this  wide  range  of  operations  are 
"explained"  in  the  most  general  control -oriented  pictures  of  pilot -vehicle 
systems  by  signal  "pathways"  internal  to  the  human.  We  shall  end  this 
introductory  section  with  an  explanation  of  this  general  paradigm  for 
human  control  behavior  as  the  underlying  basis  for  the  simpler  models 
(which  are  components  of  the  general  paradigm)  commonly  used  for  flying 
qualities  analysis.  The  general  paradigm  also  serves  to  introduce  some  of 
the  mysterious  complexities  which  face  researchers  who  wish  to  describe 
human  dynamics  in  general. 

From  this  starting  point  the  scope  will  be  contracted  to  manageable 
proportions  by  reducing  the  types  of  behavior  which  are  to  be  character¬ 
ized  by  pilot  models.  The  man-machine  systems  most  relevant  for  the 
exposure  of  critical  flying  qualities  involve  operations  in  which  the 
pilot  controls  the  effective  aircraft  dynamics  in  a  closed-loop  fashion 
"Closed- loop"  in  this  sense  means  operations  wherein  at  least  part  of  the 
pilot's  control  actions  are  conditioned  by  the  differences  between  the 
aircraft's  desired  and  actual  outputs.  The  kinds  of  piloting  covered 
include  precision  control,  regulation,  and  stabilization  tasks;  the  types 
of  flying  qualities  tests  represented  include  "flying  qualities  while 


tracking".  For  these  cases  the  man-machine  system,  and  the  associated 
human  pilot  behavior,  are  referred  to  as  "compensatory".  Fortunatel} , 
compensatory  operations  are  the  most  definitive  in  disclosing  critical 
flying  qualities  deficiencies  an''  che  associated  pilot  models  are  the  most 
extensive  and  advanced. 

There  are  currently  three  predominant  types  of  human  operator  models 
used  to  describe  compensatory  behavior.  By  far  the  simplest  model  the 
human  pilot-vehicle  system  dynamics  in  the  crossover  frequency  region. 
(The  crossover  frequency  occurs  where  the  open- loop  amplitude  ratio  of  the 
pilot-vehicle  system  is  unity.)  It  is  often  sufficient  for  flying  quali¬ 
ties  analyses  intended  to  elicit  the  governing  vehicle  parameters, 
key  variations,  and  basic  relationships.  Section  II  will  treat  this  fun¬ 
damental  model. 

The  most  elaborate  description  of  human  dynamic  properties  r  -  a  con¬ 
troller  is  the  structural  -  isomorphic  model.  This  is  an  expansion  of  the 
crossover  model  which  attempts  to  account  for  many  of  the  subsystem 
aspects  of  the  human  controller  as  well  as  the  total  input-output  beha¬ 
vior.  It  is  covered  in  Section  III.  The  fourth  section  brings  a  quite 
different  perspective  to  dynamic  pilot  modeling  in  reviewing  an  algorith¬ 
mic  model,  also  called  the  optimal  control  model  (OCM) .  The  primary  pur¬ 
pose  of  this  model  is  to  mimic  the  human  operator's  total  response  by 
appropriate  specialization  of  modern  control  computational  procedures. 
Because  the  "Crossover  Model"  is  the  most  broadly  applicable  and  best 
understood  of  human  dynamic  descriptions,  the  behavior  predicted  by  either 
the  structural - isomorphic  or  the  algorithmic  models  must  "reduce"  to  this 
form  in  the  crossover  frequency  region.  Thus  the  more  elaborate  models 
must  inevitably  return  to  the  crossover  model  as  a  necessary  limiting  case 
"consequence" . 

In  order  to  exercise  the  algorithmic  model  a  new  formulation  of  the 
computational  steps  involved  in  the  optimal  control  model  has  been  devel¬ 
oped  in  the  context  of  a  commercially  available  control  system  analysis 
program  (Program  CC) .  Besides  the  PC  compatible  format  this  new  formula¬ 
tion  includes  additional  sequences  which  allow  the  analyst  to  determine 
the  actual  estimated  pilot  characteristics  from  the  optimal  controller 
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solution.  These  steps  should  improve  the  understanding  and  interpretation 
of  algorithmic  model-based  estimates,  and  should  broaden  the  use  of  the 
OCM  by  making  it  available  as  a  PC  compatible  routine.  The  new  formula¬ 
tion  is  documented  in  Volume  III. 

In  modern  high  performance  aircraft  the  pilot  is  no  longer  primarily 
a  controller.  Instead  communications,  monitoring  and  management  of  auto¬ 
mated  equipment,  planning,  re -adjusting  to  adapt  to  changing  circum¬ 
stances,  etc.  place  increasingly  arduous  demands  on  the  pilot.  Thus  fly¬ 
ing  qualities  for  stressful,  high  workload  mission  phases  requires  the 
pilot  to  divide  his  attention  between  control  and  managerial  tasks.  The 
dynamic  models  for  the  pilot  must  take  these  divided  attention  operations 
into  account.  This  is  done  for  both  classical  and  algorithmic  models  in 
the  later  parts  of  their  respective  sections. 

The  first  four  sections  of  the  report  focus  on  the  dynamics  of  pilots 
in  flying  qualities  tasks  for  both  full  and  divided  attention  conditions. 
The  fifth,  and  last,  section  and  two  appendices  address  the  workload 
associated  with  generating  these  dynamics  and  accomplishing  the  control 
task.  This  is  measured  by  subjective  impressions  delivered  as  pilot 
ratings.  Ideally  we  would  like  to  predict  the  workload  (pilot  ratings) 
along  with  the  prediction  of  the  underlying  pilot  dynamics.  To  some 
extent  this  can  be  done,  although  much  is  left  for  the  future.  Thus  Sec¬ 
tion  V  presents  only  a  status  current  at  the  start  of  the  present  program. 
New  results,  based  on  the  new  experimental  data  obtained  as  part  of  this 
project,  are  presented  in  Volume  I. 

B.  THE  SEVERAL  NATURES  OF  MAN  MACHINE  CONTROL 

--  A  CATALOG  OF  BEHAVIORAL  COMPLEXITIES 

The  human  pilot  is  complicated  to  describe  quantitatively  because  cf 
his  enormous  versatility  as  an  information  processing  device.  Figure  1 
shows  the  general  pathways  required  to  describe  human  behavior  in  an 
interactive  man-machine  system  wherein  the  human  operates  on  visually 
sensed  inputs  and  communicates  with  the  machine  via  a  manipulative  output. 
This  control  system  block  diagram  indicates  the  minimum  number  of  the 
major  functional  signal  pathways  internal  to  the  human  operator  which  must 
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Major  Human  Pilot  Pathways  in  a  Pilot-Vehicle  System 


be  present  functionally  in  order  to  characterize  different  behavioral 
features  of  the  human  controller  The  constituent  sensing,  data  process 
ing,  computing,  and  actuating  elements  are  connected  as  internal  signal 
processing  pathways  which  can  be  reconfigured  as  the  situation  changes. 

Functional  operations  on  internal  signals  within  a  given  pathway  may  also 
be  modified.  Thus,  we  have  adaptation  both  of  the  pathways  involved  and 
of  the  functions  performed.  The  specific  internal  signal  organizational 
possibilities  shown  have  been  discovered  by  manipulating  experimental 
situations  (e.g.,  by  changing  system  inputs  and  machine  dynamics)  to  iso¬ 
late  different  combinations  of  the  specific  blocks  shown. 

To  describe  the  components  of  the  figure  start  at  the  far  right  with 
the  controlled  element;  this  is  the  machine  being  controlled  by  the 
human.  To  its  left  is  the  actual  interface  between  the  human  and  the 
machine  --  the  neuromuscular  actuation  system,  which  is  the  human's  output 
mechanism.  This  in  itself  is  a  complicated  feedback  control  system 
capable  of  operating  as  an  open-loop  or  combined  open- loop/closed- loop 
system,  although  that  level  of  complication  is  not  explicit  in  the  simple 
feedback  control  system  shown  here.  The  neuromuscular  system  comprises 
limb,  muscle,  and  manipulator  dynamics  in  the  forward  loop  and  muscle 
spindle  and  tendon  organ  ensembles  as  feedback  elements.  All  these 
elements  operate  within  the  human  at  the  level  from  the  spinal  cord  to  the 
periphery . 

There  are  other  sensor  sources,  such  as  joint  receptors  and  peri¬ 
pheral  vision,  which  indicate  limb  output  position.  These  operate  through 
higher  centers  and  are  subsumed  in  the  proprioceptive  feedback  loop  incor¬ 
porating  a  block  at  the  perceptual  level  further  to  the  left  in  the  dia¬ 
gram.  If  motion  cues  are  present  these  too  can  be  associated  in  similar  * 

proprioceptive  blocks  with  feedbacks  from  the  controlled  element  output. 

The  three  other  pathways  shewn  at  the  perceptual  level  correspond  to 
three  different  types  of  control  operations  on  the  visually  presented 
system  inputs.  Depending  on  which  pathway  is  effectively  present,  the 
control  structure  of  the  man-machine  system  can  appear  to  be  open- loop,  or 
combination  open-loop  closed-loop,  or  totally  closed-loop  with  respect  to 
visual  stimuli. 
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When  the  compensatory  block  is  appropriate  at  the  perceptual  level, 
the  human  controller  acts  in  response  to  errors  or  controlled  element 
output  quantities  only.  With  this  pathway  operational,  continuous 
closed- loop  control  is  exerted  on  the  machine  so  as  to  minimize  system 
errors  in  the  presence  of  commands  and  disturbances.  Compensatory  beha¬ 
vior  will  be  present  when  the  commands  and  disturbances  are  random¬ 
appearing  and  when  the  only  information  displayed  to  the  human  controller 
consists  of  system  errors  or  machine  outputs. 

When  the  command  inputs  can  be  distinguished  from  the  system  outputs 
by  virtue  of  the  display  (e.g.,  i  and  m  are  shown  or  are  detectable  as 
separate  entities  relative  to  a  reference)  or  preview  (e.g.,  as  in  follow¬ 
ing  a  curved  pathway),  the  pursuit  pathway  joins  the  compensatory.  This 
new  pathway  provides  an  open- loop  control  in  conjunction  with  the  compen¬ 
satory  closed-loop  error-correcting  action.  The  quality  of  the  overall 
control  can,  in  principle,  be  much  superior  to  that  where  compensatory 
acts  alone. 

An  even  higher  level  of  control  is  possible.  When  complete  familiar¬ 
ity  with  the  controlled  element  dynamics  and  the  entire  perceptual  field 
is  achieved,  the  operator  can  generate  neuromuscular  commands  which  are 
deft,  discrete,  properly  timed,  scaled,  and  sequenced  so  as  to  result  in 
machine  outputs  which  are  exactly  as  desired.  These  neuromuscular  com¬ 
mands  are  selected  from  a  repertoire  of  previously  learned  control 
movements.  They  are  conditioned  responses  which  may  be  triggered  by  the 
situation  and  the  command  and  control  quantities,  but  they  are  not  con¬ 
tinuously  dependent  on  these  quantities.  This  pure  open-loop  programraed- 
control-like  behavior  is  called  precognitive .  Like  the  pursuit  pathway, 
it  often  appears  in  company  with  the  compensatory  operations  as  a  dual¬ 
mode  control  --  a  form  where  the  control  exerted  is  initiated  and  largely 
accomplished  by  the  precognitive  action  and  then  may  be  completed  with 
compensatory  error  -  reduction  operations. 

The  above  description  of  pathways  available  for  human  control  activi¬ 
ties  has  emphasized  the  visual  modality.  Similar  behavior  patterns  are 
present  in  the  other  modalities  as  well.  Thus,  man’s  interactions  with 
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machines  can  be  even  more  extraordinarily  varied  than  described  here,  and 
can  range  completely  over  the  spectrum  from  open- loop  to  closed- loop  in 
character  in  one  or  more  modalities.  Just  what  pathways  of  the  overall 
system  are  present  at  a  particular  time  depends  on  the  detailed  nature  of 
the  specific  task  at  hand  and  the  corresponding  perceptual  situation.  All 
of  the  fundamental  pathways  are  involved  in  vor^us  piloted-aircraf  t 
maneuvers.  Thus  all  these  features  are  potentially  significant  in  vehicle 
flying  qualities.  In  the  sequel  we  shall,  however,  consider  only  the 
simplest  form  of  closed- loop  behavior  --  compensatory  operations. 
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SECTION  II 


COMPENSATORY  OPERATION  AND  THE  CROSSOVER  MODEL 

A.  COMPENSATORY  MODELS  IN  GENERAL 

The  compensatory  pathways  in  the  visual  modality  have  been  by  far  the 
most  extensively  studied  in  man-machine  systems.  Thousands  of  experiments 
have  been  performed,  and  most  of  the  adaptive  features  of  the  human  opera¬ 
tor  associated  with  these  kinds  of  operations  are  well  understood.  Both 
classical  control  and  optimal  control  theoretical  formulations  are  avail¬ 
able  to  predict  steady- state  and  dynamic  performance. 

Figure  2  illustrates  in  vector  block  diagram  form  a  general  system 
configuration  appropriate  to  closed- loop  man-machine  control.  The  diagram 
shows  the  human  operating  on  a  number  of  perceived  quantities,  y(t) ,  and 
exerting  control  over  an  aircraft  ("controlled  element")  by  actuating  a 
number  of  controls,  u^t)  .  The  response  of  the  controlled  element  to 
actuation  of  the  controls  and  to  disturbances  is  presented  on  a  "display." 
As  used  here,  display  includes  dynamic  geometrical  perspectives  of  the 
visual  field,  other  visual  stimuli  present  on  physical  display  elements 
either  on  the  aircraft  or  in  the  surround,  and  proprioceptive,  tactile, 
aural,  and  other  information  impinging  on  the  pilot.  From  the  display  the 
human  separates  the  information  needed  for  monitoring  from  that  required 
for  control  purposes.  Only  the  latter  directly  affects  the  human's  opera¬ 
tions  as  a  controller,  although  both  present  attentional  demands  and 
thereby  affect  workload. 

After  receiving  the  displayed  information  the  pilot  internally 
selects  and  equalizes  appropriate  signals  and  sends  the  results  on  to  the 
neuromuscular  actuation  subsystem  for  control  action.  The  equalization 
and  neuromuscular  properties  depend  on  the  task  variables  (effective  air¬ 
craft  dynamics,  display,  and  inputs);  they  constitute  the  pilot’s  adaptive 
features  whereby  he  attempts  to  offset  any  dynamic  deficiencies  of  the 
remaining  system  elements.  In  the  process  of  accomplishing  control  the 
human  introduces  observation,  scanning,  divided  attention,  equalization, 
and  motor  noises  which,  joined  together  are  referred  to  as  "remnant". 
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Figure  2.  A  Generalized  Man-Machine  System  Structure 


These  unwanted  components  of  the  operator's  signals  are  functions  of  the 
task  and  the  qualities  of  the  display. 

As  noted  in  the  Introduction,  two  types  of  human  operator  models  are 
available  to  handle  the  details  in  Fig.  2.  The  first  is  a  multiloop, 
multi-modality  model,  based  on  describing  functions,  which  is  structurally 
isomorphic  in  that  its  component  dynamics  are  intended  to  parallel  the 
dynamics  of  more  or  less  identifiable  human  operator  subsystems.  The 
emphasis  is  on  cause  and  effect  relationships  having  similarity  in  form 
and  structural  connections  with  those  of  the  human  operator.  The  second 
type  of  model  is  algorithmic.  It  uses  linear-quadratic-guassian  optimal 
control  theory,  modified  to  permit  a  pure  time  delay  and  operator- induced 
noises  to  be  given  quantities  along  with  the  machine  characteristics. 

Both  types  of  models  represent  the  man-machine  system  as  quasilinear 
in  the  sense  that  the  response  to  a  given  input  is  divided  into  two  parts 
--  a  component  which  corresponds  to  the  responses  of  equivalent  linear 
elements  driven  by  that  input  and  a  "remnant"  or  noise  component  which 
represents  the  difference  between  the  response  of  the  actual  system  and  an 
equivalent  system  based  on  the  linear  element.  Verbal -analytical  instruc¬ 
tions  which  express  the  adaptation  of  the  human  population  to  the  task 
variables  are  an  important  formal  feature  of  the  structural  isomorphic 
model  and  have  counter  parts,  such  as  the  specification  of  the  performance 
index,  in  the  algorithmic  model  form.  For  limited  situations,  both  repre¬ 
sentations  can  be  used  to  predict  human  operator  dynamic  behavior  (in  some 
sense),  operator- induced  noise  (remnant),  workload  indices,  visual  scan¬ 
ning  effects,  and  overall  system  performance  such  as  mean-squared  system 
errors  and  control  activities. 

The  major  fundamental  differences  between  the  models  are  their  con¬ 
ceptual  bases,  i.e.,  causal  and  structural  isomorphic  as  contrasted  to 
algorithmic  and  (potentially)  teleologic;  the  computational  techniques 
associated  with  the  exercise  of  the  model;  and  the  nature  of  model  iden¬ 
tification  processes.  At  the  present  time  there  are  other  differences 
between  the  structural  isomorphic  and  algorithmic  models  relating  to  their 
regimes  of  application  and  their  validated  capabilities  for  prediction. 
These  latter  differences  are  not,  however,  fundamental;  instead,  they 
reflect  the  relative  maturity  and  extent  of  application. 
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B.  THE  CROSSOVER  MODEL  DESCRIBING  FUNCTION 


Both  the  structural  isomorphic  and  the  algorithmic  model  approaches 
will  be  described  below.  As  a  preliminary  let  us  first  examine  some  of 
the  general  characteristics  of  human  pilot  dynamic  response  in  compensa¬ 
tory  man-machine  systems  by  considering  an  elementary  example.  Figure  3a 
shows  a  display  and  functional  block  diagram  of  a  simple  single-loop  man- 
machine  system.  The  controlled  element  dynamics  are  given  by: 

Kc 

Yc  "  s (Ts  +  1)  (1) 

This  could  represent,  for  example,  the  idealized  roll  angle  to  aileron 
transfer  function.  The  compensatory  display  presents  the  pilot  with  a 
visual  stimulus  which  shows  only  the  difference  between  the  system  forcing 
function  and  the  system  output.  (Historically  this  is  the  definition  of 
compensatory :  modern  usage  applies  the  word  compensatory  to  the  situations 
wherein  the  human  operates  on  errors  regardless  of  the  display  details.) 
The  pilot's  task  is  to  minimize  the  presented  error  signal  by  attempting 
to  keep  it  superimposed  on  a  stationary  point  or  line  on  the  display. 
This  is  accomplished  by  the  manipulative  control  action  c(t)  which  affects 
the  controlled  element,  and  gives  rise  to  the  system  output  m(t)  being 
controlled.  The  usual  purpose  of  a  system  of  this  nature  is  to  make  the 
system  output  closely  resemble  the  system  forcing  function  or,  in  other 
words,  to  make  the  output  follow  the  input.  The  quality  of  the  following 
is  indicated  by  the  system  error,  which  is,  of  course,  the  operator's 
visual  stimulus. 

Figure  3b  (Ref.  1)  presents  typical  time  histories  in  this  system 
when  a  random-appearing  forcing  function  is  applied.  The  first  thing  to 
notice  about  the  time  histories  is  that  the  system  output,  m,  does  indeed 
follow  the  forcing  function,  i,  very  closely.  Only  a  slight  time  lag 
keeps  the  output  from  being  a  nearly  identical  duplicate  of  the  forcing 
function,  although  there  are  some  small,  random  wiggles  here  and  there  on 
the  output.  On  the  other  hand,  the  operator's  output  does  not  correspond 
at  all  well  with  the  system  error,  even  if  the  error  is  delayed.  However, 
the  operator  output  lagged  by  (s  +  1/T)  is  approximately  proportional  to 
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the  error  signal  delayed  by  0.16  sec.  Thus,  as  an  approximation,  the 
operator's  transfer  characterist  : can  be  inferred  to  be: 

Yp  i  Kp(Ts  +  l)e-TS  (2) 


This  result  states  that  the  operator  develops  a  lead  which  is  approxi¬ 
mately  equal  to  the  first-order  lag  component  of  the  controlled  element 
dynamics  and  that  the  operator's  response  lags  his  stimulus  by  r  sec.  The 
open-loop  man-machine  transfer  characteristic  appears  as: 

[Kp(Ts  +  le-TS ]KC 

Vc  "  s  (Ts  +  1)  (3) 


or 


KpKce~TS 


G  *  YpYc 


(4) 


The  data  of  Fig.  4  illustrate  how  well  Eq.  4  is  obeyed  for  a  particular 
Yc  and  a  variety  of  subjects.  The  agreement  with  the  amplitude  ratio  is 
excellent  over  a  broad  range  of  frequencies.  The  phase  agreement  is  good 
in  the  region  of  the  crossover  frequency,  wc ,  but  departs  somewhat  at 
lower  frequencies. 

If  now  a  large  variety  of  controlled  element  forms  are  used  and  simi¬ 
lar  measurements  are  taken,  the  human  transfer  characteristics  will  be 
different  for  each  controlled  element.  But,  for  a  very  wide  range  of 
controlled  element  dynamics,  the  form  of  the  total  open- loop  transfer 
characteristic  about  the  crossover  frequency  will  remain  substantially 
invariant.  In  other  words,  experiment  shows  that  Eq .  4  has  some  preten¬ 
sion  to  general  applicability.  The  effective  time  delay,  r,  which  is  of 
course  only  a  low-frequency  approximation  to  all  manner  of  high-frequency 
leads  and  lags,  is  not  a  constant.  It  depends  primarily  on  the  amount  of 
lead  equalization  required  of  the  operator,  as  shown  in  Fig.  5  (Ref.  1). 
This  indicates  that  pilot  equalization  to  offset  controlled  element 
dynamic  deficiencies  has  an  associated  computational  time  penalty.  With 
this  proviso  on  r,  the  Eq .  4  relationship  becomes  the  well-known 
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Figure  4.  Data  and  Crossover  Models  for  a  Simple 
Rate-Control-Like  Controlled  Element 
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Figure  5.  Variation  of  Crossover  Model  Dynamic 
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of  Operator  Lead  Equalization 


simplified  crossover  model  of  compensatory  manual  control  theory.  For 
the  vast  majority  of  flying  qualities  applications  the  pilot  lead 
equalization  ranges  from  none  to  first  order.  Thus,  in  order  to  estimate 
the  effective  time  delay,  re,  we  need  only  to  approximate  that  portion  of 
Fig.  5  in  the  interval 


0  < 


d(  Y. 


P  dB' 


d(log  w) 


<  20  dB/decade 


This  is  given  by 


1_ 

r  i 

1  .  ,  0.833  d( 

!yp 

dB) 

ro 

sec 

|  20  d( log 

w) 

The  human  operator's  adaptation  to  controlled  element  dynamics  is 
implicit  in  the  relationship ,  i.e.,  for  a  particular  set  of  controlled 
element  dynamics  defined  by  Yc  the  human  will  adopt  a  crossover  region 
transfer  characteristic  Yp  £  wce~Jwf/j u  Yc .  The  general  form  of  the 
human's  response  would  thus  be  determined  by  the  specifics  of  Yc ,  and 
changes  in  this  task  variable  evoke  changes  in  Yp ,  such  that  the  crossover 
model  open-loop  transfer  characteristic  form  is  preserved.  An  extremely 
important  consequence  is  that  a  duality  principle  is  established  for  those 
cases  where  the  crossover  model  is  valid.  That  is,  since  Yc  and  Yp  are 
related  by  Eq .  U,  properties  of  e i  ther  can  be  used  as  "coordinates"  in 
describing  results  of  experiments.  In  particular  one  can  use  open- loon 
effective  aircraft  dynamics  to  specify  what  is  needed  to  achieve  "good" 
pilot  dynamics  (and  subjective  ratings)  in  closed- loop  operations. 

The  crossover  model  also  applies  when  the  machine  dynamics  are 
smoothly  time  varying  (Ref.  2).  The  crossover  frequency  tends  to  be  con¬ 
stant  for  a  given  set  of  task  variables.  It  increases  slightly  as  forcing 
function  bandwidth  is  increased  and  is  reduced  for  very  small  input  ampli¬ 
tudes.  This  is  a  consequence  of  the  operator's  indifference  threshold, 
which  is  the  most  important  nonlinearity  to  be  considered  in  connection 
with  crossover  model  transfer  characteristics. 

Most  of  the  interesting  open  and  closed  loop  dynamic  characteristics 
of  the  crossover  model  are  summarized  in  Fig.  6  (adapted  from  Ref.  1). 
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Because  only  two  parameters  are  present  many  features  can  be  presented  as 
functions  of  non-dimensional  quantities,  such  as  the  normalized  crossover 
frequency,  wcr ,  which  is  also  the  normalized  open-loop  gain.  The  system 
survey  summarizes  the  key  open  and  closed  loop  features.  For  instance, 
the  most  significant  property  ol  feedback  systems  is  that  which  prevails 
when  the  open- loop  transfer  characteristic  is  much  larger  than  1,  for  then 
the  system  output  is  almost  exactly  equal  to  the  system  input  and  the 
system  error  is  very  small.  From  examination  of  the  G(jw)  Bode  plot  it 
can  be  seen  that  this  occurs  at  low  frequencies  and  is,  of  course,  the 
reason  that  the  system  output  and  system  forcing  function  were  similar  in 
the  example  of  Fig.  4.  At  high  frequencies  |G(ju|  «  1,  so  the  closed- 
loop  relationship  between  M(ju)  and  I(jw)  is  substantially  the  same  as  the 
open-loop,  i.e.,  the  feedback  loop  is  effectively  open.  For  the  crossover 
model,  the  frequency  which  divides  these  two  regimes  of  near- ideal  follow¬ 
ing  of  the  forcing  function  and  little  or  no  feedback  action  is  the  cross¬ 
over  frequency,  u>c .  In  Fig.  6a  this  is  the  intersection  of  the  0  dB  line 
established  for  a  particular  gain  with  the  G(jw)  Bode  plot.  The  name 
"crossover  frequency"  comes  from  this  crossover  intersection  of  the  0  dB 
line  by  the  open- loop  frequency  response  characteristic.  For  stable  oper¬ 
ation  of  the  system  the  normalized  crossover  frequency,  rwc ,  can  range 
from  0  to  w/2 . 

When  tu>c  is  relatively  small  compared  to  n/2 ,  then  rwc  is  also  tanta¬ 
mount  to  the  closed-loop  system  "bandwidth"  (the  frequency  at  which  the 
output  amplitude  is  3  dB  less  than  the  amplitude  of  an  input  sinusoid)  . 
For  higher  crossover  frequencies  this  direct  equivalence  between  ru>c  and 
system  bandwidth  degrades  because  of  a  peaking  in  the  closed-loop  fre¬ 
quency  response  near  the  crossover  frequency,  but  even  for  this  kind  of  a 
system  rwc  is  always  equal  to  or  less  than  the  bandwidth  and  thus  provides 
a  lower  bound. 

The  root  locus  plots  emphasize  the  closed- loop  system  roots.  The 
root  locus  plots  shown  in  Fig.  6a  indicate  that  the  root  which  starts  at 
the  origin  for  open-loop  (zero  gain)  conditions  progresses  further  into 
the  left  half  plane  as  the  gain  tuc  is  increased.  At  the  point  to  -  -1 
(for  which  ruc  =  1/e),  this  branch  of  the  locus  meets  with  the  first  of  an 
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infinity  of  branches  present  because  of  the  e-rs  term  (the  other  branches 
are  not  shown) .  The  quadratic  formed  by  the  two  branches  then  increases 
in  undamped  natural  frequency  and  decreases  in  damping  as  gain  is 
increased  until  neutral  stability  occurs  at  the  gain  rwc  -  jt/2  . 

Some  connections  between  the  dominant  closed- loop  mode  characteris¬ 
tics,  i.e.  f  and  wn,  and  the  open- loop  characteristics,  wc  and  r,  are 
illustrated  in  Fig.  6b.  Other  relationships,  such  as  the  phase  margin, 
and  gain  margin,  G^,  are  recapitulated  in  Fig.  6d. 

An  "extended  crossover  model"  form  is  sometimes  used  to  improve  the 
low  frequency  phase  fit.  It  is, 


Figure  4  illustrates  how  well  the  additional  term,  e~jQ/w,  describes  those 
phase  contributions  in  the  crossover  region  which  arise  from  leads  and 
lags  (in  the  pilot  and/or  the  rest  of  the  system)  which  are  present  well 
below  the  crossover  frequency  band.  The  phase  contribution  represented 
by  e~ J0/00  in  Eq .  5,  is  not  intended  to  extend  to  extremely  low  frequen¬ 
cies  . 

C .  REMNANT 

The  second  component  of  the  operator's  response  is  operator- induced 
noise  or  remnant.  This  can,  in  principle,  result  from  several  sources, 
but  in  single- loop  systems  with  linear  manipulators  the  basic  cause 
appears  to  be  random  time -varying  behavior  within  the  operator  primarily 
associated  with  fluctuations  in  the  effective  time  delay.  This  can  be 
interpreted  as  a  random  change  in  phase,  akin  to  a  random  frequency  modu¬ 
lation,  or  to  variations  of  internal  sampling  rate  in  a  sampled  data 
interpretation  of  the  operator  (Refs.  1,  3-6).  In  any  event,  the  remnant 
is  a  continuous,  relatively  broadband,  power  spectral  density  which,  as 
shown  in  Fig.  7,  scales  approximately  with  the  mean-squared  error 
(Refs.  4,  5). 


u)  (rod /sec) 


Figure  7.  Normalized  Remnant  Spectra 

Task  variables  other  than  the  machine  dynamics,  as  well  as  environ¬ 
mental  and  operator-centered  variables,  can  change  open- loop  gain,  effec¬ 
tive  time  delay,  and  remnant.  Accordingly,  wc  and  r  variations  become  a 
quantification  of  changes  or  differences  in  the  task,  environmental,  and 
operator-centered  variables  expressed  directly  in  terms  of  the  operator's 
control  actions.  In  measuring  the  effects  of  training  for  instance,  u>c 
increases  with  trials  until  stable  conditions  are  obtained  for  that  par 
ticular  subject  and  set  of  constant  task  and  environmental  variables. 
Similarly,  the  remnant  may  also  change  as  a  function  of  the  control  situa¬ 
tions.  For  instance,  comparison  of  Figs.  7a  and  7b  shows  the  change  in 
remnant  bandwidth  and  level  associated  with  the  lead  equalization  required 
to  offset  controlled  element  lags.  As  another  example  Ref.  7  shows  that 
operator  gain  is  decreased  and  remnant  is  increased  as  a  consequence  of 
ingested  alcohol. 
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To  generalize  these  remarks,  the  total  pilot  actions  can  be  thought 
of  as  that  of  an  adaptive  plastic  sensory-motor  link  --  adaptive  in  that 
the  pilot  is  task-adjusted  to  offset  controlled  element  dynamic  deficien¬ 
cies  and  to  respond  to  forcing  function  commands  or  regulate  against  dis¬ 
turbances;  plastic  in  that  the  adaptive  characteristics  are  further  shaped 
by  the  external  and  internal  (pilot-centered)  environments.  These  beha¬ 
vioral  features  must  be  accounted  for  in  either  the  structural  isomorphic 
or  algorithmic  models.  A  general  description  of  these  models  and  some  of 
their  characteristics  follows. 


SECTION  III 


THE  STRUCTURAL  ISOMORPHIC  HUMAN  OPERATOR  MODEL 

A.  BACKGROUND 

The  extensive  analytical  and  experimental  studies  of  closed- loop  man- 
machire  systems  conducted  since  World  War  II  have  had  as  a  principal  goal 
the  mathematical  quantification  of  human  dynamic  behavior  and  the  develop¬ 
ment  of  laws  which  permit  this  behavior  to  be  predicted.  In  general, 
emphasis  has  been  on  the  human  operator  as  a  complete  entity  rather  than 
as  a  summation  of  functional  subsystems. 

In  recent  years,  the  precision  and  dynamic  range  of  measurements 
taken  with  the  total  human  operator  have  increased  greatly  --  to  the  point 
that  certain  of  the  measurements  made  over  certain  frequency  ranges  can  be 
associated  with  the  human  subsystem  dynamics.  Thus,  the  study  of  the 
human  operator  as  a  whole  has  now  arrived  at  the  stage  where  not  only  must 
subsystem  models  sum  up  to  be  compatible  with  the  total  human  dynamic 
model,  but  subsystem  and  total  system  studies  can  be  directly  related. 
Accordingly,  control  engineering  descriptions  of  the  overall  human  (see, 
e.g.,  the  list  of  surveys),  dynamical  descriptions  of  the  human  motor 
coordination  system,  studies  of  predictive  control  conducted  for  physiolo¬ 
gical  understanding,  and  studies  of  neuromuscular  actuation  systems,  which 
were  originally  separated  disciplines,  now  become  united. 

As  described  in  Ref.  8,  the  adaptive  and  plastic  properties  of  the 
operator  permit  the  experimenter  to  set  the  stage  and  write  a  script 
calling  for  a  particular  form  of  action.  Table  1  illustrates  some  of  the 
experimental  procedures  which  can  be  used  to  evoke  various  types  of 
behavior . 

B.  THE  FULL  ATTENTION  MODEL 

By  properly  selecting  combinations  of  these  procedures  and  tech¬ 
niques,  particular  channels  of  human  dynamic  operations  can  be  isolated, 
examined,  and  measured.  Appropriate  models  which  "explain"  each  of  these 
varieties  of  behavior  and  which  are  also  compatible  with  what  is  known 
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TABLE  1.  EXPERIMENTAL  PROCEDURES  TO  EVOKE  HUMAN  OPERATOR  BEHAVIORAL  CHANGES 


ially  affect  various  internal  signal  pathways 


from  other  views  of  experimental  psychology  and  physiology  can  then  be 
constructed  to  form  a  current  version  of  the  structural  isomorphic  model. 
One  such  construction,  which  is  somewhat  simplified,  is  given  in  Fig.  8. 
Here  the  controlled  element  and  display  blocks  constitute  the  machine, 
whereas  all  the  remaining  detail  reflects  the  man. 

Starting  at  the  far  right  is  the  neuromuscular  actuation  system. 
Because  the  man-machine  system  depicted  here  is  operating  on  random¬ 
appearing  signals  which  have  essentially  stationary  statistics,  the  neuro¬ 
muscular  system  is  fluctuating  about  an  operating  point  which  in  general 
corresponds  to  some  steady-state  or  average  tension.  This  is  graphically 
illustrated  by  examination  of  the  average  and  differential  EMG  signals 
shown  in  Fig.  3b.  Consequently,  the  dynamic  operations  of  muscles,  which 
can  act  only  in  contraction,  can  be  treated  as  positive  or  negative  fluc¬ 
tuations  of  many  agonist/antagonist  pairs  about  a  steady  tension  bias 
value.  This  permits  a  great  simplification  in  depicting  the  dynamic 
essentials  in  terms  of  a  block  diagram.  The  forward  path  of  the  neuromus¬ 
cular  system  shown  includes  ensembles  of  muscles  operating  on  coupled 
skeletal  and  manipulator  dynamics.  The  feedback  path  sensors  operating  at 
the  spinal  level  are  primarily  spindle  and  Golgi  tendon  organs.  Because 
the  individual  actions  of  specific  sensors  are  difficult  to  separate  in 
the  intact  human  the  system  shown  has  a  feedback  element  labeled  as 
spindle/tendon  organ  ensembles.  The  spindle  characteristics  may  very  well 
be  predominant  for  the  small  motions  and  relatively  light  forces  involved 
in  most  of  the  measurements  thus  far  accomplished.  The  effective  dynamics 
of  the  closed- loop  neuromuscular  system  from  the  alpha  motor  neuron 
command  signals  to  manipulator  force  can  be  approximated  over  a  wide  fre¬ 
quency  range  by  the  third-order  transfer  function  shown.  This  form  is 
also  compatible  with  small  perturbation  dynamics  based  on  experimentally 
verified  analytical  models  of  muscle  and  manipulator  characteristics 
(Refs.  9,  10).  The  parameter  values  are  strongly  dependent  on  the  steady- 
state  neuromuscular  tension,  70 ,  due  to  the  gamma  motor  system.  The  gamma 
commands  also  affect  the  dynamics  of  the  spindle  ensembles  and,  in  fact, 
provide  another  pathway  (not  shown)  capable  of  actuating  the  neuromuscular 
system  via  the  spindle  ensembles.  These  features  are  pictured  by  the 
arrows  indicating  variation  in  the  ZSp  and  PSp  factors  in  the  neuromus¬ 
cular  system  feedback  block  and  in  the  7^  and  70  inputs. 
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fiure  8.  Structural  Isomorphic  Mode 


This  rudimentary  level  of  neuromuscular  actuation  system  description 
is  a  minimum  to  have  value  even  in  gross  physiological  descriptions.  It 
is  an  essential  feature  in  the  study  of  human  pilot  characteristics  in 
vibratory  environments  (Ref.  11)  and  is  also  often  needed  for  the  study  of 
limb/manipulator  system  dynamics  in  aircraft  control  (e.g.,  Refs.  12,  13). 
For  many  other  man-machine  system  applications,  however,  the  neuromuscular 
actuation  dynamics  are  so  high  in  frequency  as  to  be  relatively  unimpor¬ 
tant  in  their  details.  In  these  cases,  a  pure  time  delay,  r^,  or  a 
first-order  lag  can  be  used  as  a  low-frequency  approximation. 

The  neuromuscular  actuation  system  described  thus  far  is  appropriate 
when  the  manipulator  is  restrained  by  a  stiff  spring  and  the  control 
actions  involve  very  little  joint  movement.  When  significant  joint  move¬ 
ments  are  present,  proprioceptive  pathway  elements  enter  into  the  neuro¬ 
muscular  actuation  system  dynamics.  These  derive  from  several  sources, 
the  most  important  being  peripheral  vision  and  joint  receptors  in  the 
limb.  These  feedbacks  act  through  higher  centers  and  thereby  exhibit 
larger  response  time  delays.  When  they  are  present,  the  neuromuscular 
actuation  system  bandwidth  may  be  reduced  significantly. 

Proceed  now  to  the  sensory  mechanisms  at  the  far  left  of  the  human 
operator.  A  good  deal  of  the  detail  in  the  visual  pathway  is  intended  to 
emphasize  the  parallel  operations  of  parafoveal  and  foveal  vision  and  the 
control  of  eye  movements.  An  important  feature  of  the  visual  pathways  is 
that  essentially  continuous  signals  from  a  particular  display  element  can 
be  available  to  the  operator,  by  virtue  of  the  parallel  foveal  and  para¬ 
foveal  pathways,  even  when  the  eye  is  scanning.  The  essence  of  past  work 
in  man-machine  systems  involving  many  displays  (Refs.  1,  14-18)  shows 

that : 

1.  A  fairly  stationary  scanning  strategy  evolves  for  a  given 
task  and  display  array. 

2.  The  operator's  output  control  motions  are  much  more  con¬ 
tinuous  than  a  discrete  sampling  of  input  signals  coinci¬ 
dent  with  foveal  eye  fixations  would  imply. 

3.  The  first-order  effects  of  scanning  are  to  reduce  gain  and 
increase  remnant  in  the  scanned  channels. 
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The  degree  of  gain  reduction  depends  on  parafoveal  viewing  angle  and  rela¬ 
tive  parafoveal  to  foveal  dwell  times. 

The  other  sensory  elements  are  vestibular  and  kinesthetic  (Refs.  19- 
23),  which  are  present  when  the  pilot  is  moving,  as  in  a  maneuvering  air¬ 
plane  or  a  moving  base  simulator.  The  pilot  contains  neurological 
elements  capable  of  sensing  rotary  and  linear  accelerations.  These  are 
primarily  in  the  vestibular  apparatus,  although  other  sensors  and  pathways 
can  also  be  involved.  The  rotary  motion  feedbacks  usually  associated  with 
the  semicircular  canals  act  like  signals  from  a  highly  overdamped  angular 
accelerometer.  Over  the  frequency  range  from  about  0.2  to  10  rad/sec  the 
output  signal  is  proportional  to  angular  rate,  so  the  sensor  can  function 
as  a  rate  gyro.  For  prolonged  steady  turning  the  sensor  washes  out;  thus, 
spurious  sensations  occur  in  steady  rotations  or  when  the  turning  motion 
stops.  This  pathway  has  a  threshold  on  the  order  of  1-2  deg/sec.  Because 
the  rotary  motion  sensing  apparatus  gives  rise  to  an  angular-rate-like  cue 
directly,  any  need  for  generating  angular  rate  information  by  means  of  a 
lead  equalized  visual  cue  may  be  reduced.  This  feedback  can  also  be 
thought  of  as  an  inner  loop  which  tends  to  reduce  the  effective  operator 
time  delay.  That  is,  in  terms  of  crossover  model  characteristics  measured 
as  if  visual  pathways  only  were  active,  the  presence  of  rotary  motion  can 
reduce  the  effective  time  delay  by  as  much  as  0.1  sec. 

The  other  functional  operation  of  the  vestibular  and  kinesthetic 
pathways  is  the  provision  of  the  "nystagmus  crossfeeds"  to  the  oculomotor 
system.  These  produce  involuntary  eye  motions  as  a  function  of  the  exci¬ 
tation  of  the  vestibular  apparatus.  These  eye  movements  can  be  helpful  in 
properly  directing  the  gaze,  although  many  of  their  most  interesting  pro¬ 
perties  involve  their  effects  in  disorientation  and  illusions.  The  motion 
effects  which  conflict  with  the  visual  modality  can  seriously  distort  the 
operator's  perception  of  the  state  of  affairs  and  can  be  so  severe  as  to 
affect  the  human's  control  capacity. 

Turn  now  to  the  central  elements.  As  shown  there,  the  operator  can 
develop  a  neuromuscular  system  input  command  which  is  the  summation  of  a 
lag,  proportional,  lead,  and  double -lead  function  of  the  system  error. 
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The  lag  and  proportional  channels  have  a  basic  time  delay,  rc,  associated 
with  them.  The  higher  derivative  channels  have  additional  incremental 
delays.  These  incremental  time  delays  constitute  the  dynamic  cost  of  lead 
generation.  They  are  about  1/5  sec  for  rate,  rp,  and  greater  than  1/2  sec 
for  the  acceleration  channel,  r^.  The  proportional,  rate,  and  accelera¬ 
tion  equalizations  are  shown  as  separate  parallel  channels  primarily 
because  of  their  respective  latency  differences.  This  independence  of 
these  channels  is  oversimplified,  for  common  neurological  apparatus  is 
undoubtedly  present  for  each  function.  These  common  elements  are  modeled 
here  by  the  central  processing  and  integration  block  preceding  the  visual 
channel  and  the  motor  command  integrative  mechanisms  succeeding  it. 
Besides  the  different  time  delays,  the  other  evidence  for  parallel  chan¬ 
nels  is  the  difference  in  response  quality  as  a  function  of  the  low- 
frequency  equalization  supplied  by  the  operator.  For  example,  when  very- 
low-frequency  leads  are  present,  as  if  operations  were  through  the  rate  or 
acceleration  channels,  the  operator's  output  tends  to  be  more  discrete  and 
pulselike  than  when  little  or  no  lead  is  required. 

The  channel  gains  and  the  lag  time  constant  Tj  are  all  shown  as  vari¬ 
able  quantities.  These,  in  conjunction  with  the  neuromuscular  system 
variations  with  y0,  constitute  the  principal  adaptive  changes  in  the  oper¬ 
ator  characteristics  as  display,  controlled  element,  and  environmental 
conditions  change.  For  a  given  controlled  element,  these  are  adjusted 
such  that  the  crossover  model  applies  over  its  frequency  range  of  valid¬ 
ity.  Thus,  the  extremely  complicated  structural  isomorphic  model  reduces 
to  the  visual  and/or  vestibular  equalization  actually  present  and  with 
neuromuscular  dynamics  as  pertinent  to  the  task.  When  a  higher  degree  of 
exactitude  is  required,  the  structural  isomorphic  model  is  adjusted  via  a 
series  of  analytical/verbal  rules  which  take  into  account  the  details  of 
the  task  variables.  A  version  of  these  rules  is  summarized  below. 

C.  EQUALIZATION  SELECTION  AND  ADJUSTMENT 

For  essentially  all  flying  qualities  while  tracking  aircraft  applica¬ 
tions  the  form, 

K(TJ  w  +  1) 

(Txj  w  +  1) 
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(6) 


is  used  in  lieu  of  the  multi  channel  parallel  equalization  elements.  This 
form  doesn't  reflect  the  acceleration  channel,  which  is  seldom  used  in 
aircraft  control.  For  a  specific  problem  a  particular  equalization  is 
selected  from  the  general  form  such  that  the  following  properties  obtain: 

(a)  The  system  can  be  stabilized  by  proper  selection  of  gain, 
preferably  over  a  very  broad  region. 

(b)  Over  a  considerable  frequency  range  in  the  unit  gain  cross¬ 
over  region  (that  frequency  band  centered  on  the  crossover 
frequency,  wc) ,  the  open-loop  describing  function 
|YpYc( jw) 1 | dg  has  approximately  a  -20  dB/decade  slope. 

(c)  |YpYc(jw)|  »  1  at  low  frequencies  to  provide  good  low- 

frequency  closed- loop  response  to  system  forcing  functions 
(commands )  . 

Examples  of  form  selection  and  basic  adjustment  are  provided  in  Table  2. 

1.  Time  Delay  Adjustment 

Examples  of  time  delay  adjustment  appropriate  for  aircraft  are  listed 
in  Table  3.  (This  detailed  breakdown  replaces  the  Fig.  5  data  for  the 
lead-lag  equalization  form  of  Eq .  6.)  The  visual  lag  and  proportional 
channels  have  a  basic  (minimum)  time  delay,  rv,  of  0.1  sec  associated  with 
either  or  both  of  them  when  all  other  effects  (e.g.,  motion  sensing,  full 
limb/manipulator  neuromuscular  system,  and  display  computational  lags)  are 
represented  separately;  rv  should  be  increased  to  0.2  sec,  if  fixed-base 
operations  are  being  considered  with  visual  lag  and/or  proportional  equal¬ 
ization,  full  neuromuscular  system  and  separate  display  effects.  If  the 
problem  emphasis  is  on  low  frequency  phenomena,  the  neuromuscular  system 
can  be  approximated  by  a  pure  delay;  then  T^  is  added  to  rv,  where 
examples  of  values  for  T^  are  given  in  Table  3.  Lead  equalization 
developed  on  visual  inputs  gives  rise  to  an  additional  incremental  delay, 
Arv.  This  incremental  time  delay  constitutes  the  dynamic  cost  of  pilot 
lead  generation  in  the  visual  modality . 

2.  Crossover  Frequency  with  Full  Attention 

The  factors  involved  in  estimating  crossover  frequency,  o>c ,  with  full 
attention  to  control  activity  consist  of  the  following: 
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TABLE  2.  TYPICAL  PILOT  EQUALIZATION  CHARACTERISTICS 


TABLE  3.  TIME  DELAY  ADJUSTMENT 


(a)  wc  estimates  for  rectangular  and  quasi-rectangular  forcing 
function  spectra  (discrete  power-spectral  densities  that 
are  essentially  rectangular  and  low-pass  continuous  spectra 
with  a  high-frequency  cutoff  equivalent  to  a  third-or 
higher-order  lag  filter) .  For  these  forcing  functions  an 
effective  forcing  function  bandwidth,  ,  can  be  defined 
as , 


[i 


o  ^ii(w)dw 


•] 


fo  [<i>i i Cw)]  dw 


(7) 


(1)  Basic  crossover  frequency,  wc  .  The  basic  cross¬ 
over  frequency  for  quasi-rectangular  forcing 
function  spectra  is  found  by  adding  the  phase 
angle,  -ur0,  due  to  the  human's  base  effective 
time  delay,  to  the  phase  angles  of  the  controlled 
element  and  the  previously  estimated  Yp  equalizer 
characteristics.  The  net  effect  of  controlled 
element  high-frequency  (well  past  crossover) 
leads  and  lags  is  approximated  by  an  effective 
time  delay,  rc .  Estimates  for  uc  and  the  asso¬ 
ciated  pilot  gain  are  then  made  Prom  the  condi¬ 
tions  for  neutral  stability, 


2(r0  +  Tc) 


(8) 


(2)  Phase  margin,  4m-  The  phase  margin  for  this 
forcing  function  category  corresponds  to  an 
incremental  time  delay,  Are(w^e). 
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(b)  uc  Estimates  for  low-pass  with  a  roll-off  of  less  than 
third-order  and  augmented  (shelf -type)  continuous  input 
spectra.  For  these  types  of  forcing  functions  the  cross¬ 
over  frequency  is  adjusted  to  minimize  the  mean-square 
error  in  the  presence  of  remnant. 
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Continuous  attention  remnant.  Approximations  to 
the  forms  of  injected  remnant,  3>nn  ,  when 
reflected  to  the  pilot's  input  signal  untler  con¬ 
ditions  of  continuous  attention  were  shown  in 
Fig.  7. 

(0.1  to  0.5)  where  integral  and/or  proportional 
(u>2  4  32)  equalization  are  used 


(0.1  to  0.5) 

(w2  +  1) 


where  lead  equalization  is  used 
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(2)  Nominal  crossover  frequency,  wc.  With  equaliza¬ 
tion  and  effective  time  delay,  re,  selected  as 
above,  the  nominal  crossover  frequency,  u>c ,  and 
associated  pilot  gain  is  estimated  from  the  con¬ 
dition  to  provide  minimum  mean-squared  error  in 
the  presence  of  the  appropriate  form  of  contin¬ 
uous  attention  remnant  in  Item  (b)(1)  above.  The 
nominal  cases  (continuous  remnant  magnitude  set 
to  the  geometric  mean  of  the  values  cited  above) 
are : 

^c/^u 

No  Pilot  Lead:  0.783 

Low-Freauencv  Pilot  Lead:  0.662 

where  wu  is  the  maximum  full  attention  crossover  frequency  at  the  maximum 
full  attention  crossover  frequency  at  the  dynamic  stability  limit  corres¬ 
ponding  to  zero  phase  margin  (<^^  -  0).  Thus  wu  -  ?r/2 re. 

(c)  Nominal  crossover  frequency  regression.  When  nears  or 
becomes  greater  than  0.8  wc  for  the  quasi -  rectangular 
forcing  function  case  or  when  /u>c  is  greater  than  1  for 
the  low-pass  and  augmented  low-pass  spectra,  then  the 
cross-over  frequency  regresses  to  values  much  lower  than 
wCo  and  wc ,  respectively. 

(d)  Nominal  crossover  frequency  invariance  properties. 

(1)  wc  -  Kc  independence.  After  initial  adjustment, 
changes  in  controlled  element  gain,  Kc ,  are 
offset  by  changes  in  pilot  gain,  Kp,  i.e., 
nominal  crossover  frequency,  u>c ,  is  invariant 
with  Kc . 

(2)  uc  -  independence.  Nominal  crossover  fre¬ 

quency  increases  only  slightly  with  forcing  func¬ 
tion  bandwidth  until  crossover  frequency  regres¬ 
sion  occurs.  (Phase  margin,  of  course,  increases 
per  Eq .  9 . ) 

(e)  Threshold  properties.  With  very  low  stimulus  amplitudes,  a 
threshold  characteristic  should  be  included  in  series  with 
the  pilot's  describing  function.  Also,  when  full- 
attention,  nearly  continuous  control  actions  are  not 
required,  an  indifference  threshold  is  likely  to  be 
present.  Both  of  these  lower  wc  from  what  would  be  esti¬ 

mated  using  the  above  adjustment  rules. 
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The  wc  regression  phenomenon  mentioned  in  the  adjustment  rules  refers 
to  a  reduction  of  pilot  gain  and,  hence,  of  crossover  frequency  when  the 
forcing  function  bandwidth  becomes  too  large.  The  reason  for  this  is  best 
described  by  referring  to  the  relative  mean-squared  error  plotted  in 
Fig.  9  for  the  crossover  model  subjected  to  a  rectangular  forcing  function 
spectrum  with  bandwidth  .  If  the  ratio  is  less  than  about  0.8,  an 
increase  in  normalized  gain  (fewc)  will  result  in  a  decrease  in  normalized 
mean-squared  error.  When  this  approximate  inequality  is  reversed,  the 
normalized  mean-squared  error  can  become  greater  than  1  as  gain  is 
increased.  The  trend,  therefore,  for  high  forcing  function  bandwidths  is 
to  reduce  gain.  This  regression  effect  has  practical  consequences  when¬ 
ever  the  pilot  is  required  to  track  broadband  signals. 

The  adjustment  rules  given  above  are  generally  adequate  for  the 
pilot's  lower -frequency  dynamics  in  tasks  with  spring-restrained  manipula¬ 
tors.  The  higher- frequency  properties  due  primarily  to  the  neuromuscular 
actuation  system  are  included  only  to  the  extent  that  Tjqn  is  a  component 
of  r  e . 

The  neuromuscular  system  dynamics  will  change  markedly  as  the  manipu¬ 
lator  load  dynamics  are  modified.  One  of  the  most  important  of  these 
possible  modifications  is  reduction  in  stiffness  of  the  spring  restraints. 
This  is  a  common  feature  of  aileron  controls,  as  opposed  to  elevator  and 
rudder  controls.  When  the  spring  forces  are  light,  the  manipulator 
approaches  the  free-moving  (isotonic)  extreme.  In  these  cases,  the  pilot 
must  supply  proprioceptive  feedbacks  that  introduce  into  the  neuromuscular 
system  dynamics  additional  delays  that  are  not  present  with  the  isometric 
situation.  Available  data  from  Refs.  10,  24,  25  indicate  that  the  effect 
of  this  proprioceptive  feedback  required  of  the  pilot  when  the  manipulator 
is  free-moving  is  to  increase  the  effective  time  delay  by  approximately 
0.1  sec.  This  can  be  added  directly  to  the  previously  discussed  time 
delay,  rD.  It  amounts  to  an  additional  time  delay  cost  incurred  by 
forcing  the  pilot  to  close  a  positional  loop  about  the  manipulator. 

For  some  configurations  of  manipulator  and  effective  vehicle 
dynamics,  the  higher  frequency  characteristics  of  the  neuromuscular  system 
can  be  important.  In  particular,  the  peaking  tendency  associated  with  the 
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Stability  Limit 


second-order  mode  in  the  limb/manipulator  block  of  Fig.  7  can  be  suffi¬ 
ciently  large  to  make  a  higher  frequency  gain  margin  (in  the  frequency 
range  from  2  to  3  Hz)  negligible  or  even  negative.  Whether  this  will  lead 
to  an  instability  will  depend  on  the  accompanying  phase.  Such  very  high 
frequency  pilot-effective  vehicle  oscillations  as  "roll  ratchet"  can  be 
caused  by  this  coupling.  The  detailed  nature  of  the  peaking  tendency  is  a 
very  strong  function  of  the  manipulator  and  the  rest  of  the  controlled 
element  dynamics.  The  peak  can  be  "tuned”  to  a  maximum  or  minimum  by  the 
presence  of  just  the  right  amount  of  controlled  element  lag.  Thus,  for 
example,  a  pure  Yc  -  K/s  will  have  little  if  any  peaking  while  a  Yc  - 
K/s(Ts  +  1),  with  T  about  0.1  sec,  will  have  a  great  deal.  The  known 
connections  are  all  empirical;  therefore,  the  reader  is  referred  to 
Refs.  12  and  13,  which  present  all  of  the  available  data. 

Another  "structural  model"  of  the  human  pilot  has  been  fruitfully 
applied  to  flying  qualities  problems  (Refs.  26-29).  This  model  makes  most 
of  the  adjustments  of  the  pilot  equalization  via  feedback  pathways  instead 
of  in  the  forward  loop,  and  the  "isomorphic"  features  and  the  incremental 
time  delays  associated  with  rate  and  acceleration  equalization  are  not 
modeled.  A  good  deal  of  effort  has  been  spent  on  validation  with  the 
existing  data  base,  and  with  developing  connections  with  pilot  ratings  via 
the  theory  of  Ref.  29. 

Having  completed  this  review  of  the  structural- isomorphic  and  cross¬ 
over  pilot  models  for  full  attention  situations,  we  will  next  examine 
relationships  affecting  pilot  dynamics  in  divided  attention  situations 
involving  control  operations. 

D.  DIVIDED  ATTENTION  PILOT -VEHICLE -TASK  MODEL 

The  pilot  is,  in  general,  involved  in  two  types  of  operations  -- 
control  tasks  and  a  diverse  combination  of  monitoring/supervising/ 
communicating/data-gathering/decision-making  activities  referred  to  as 
"managerial"  tasks.  While  the  pilot's  attention  is  "divided"  between  the 
"control"  and  "managerial"  tasks,  these  are  often  performed  nearly  simul¬ 
taneously  as  parallel  processing  operations.  Neither  type  of  task  is 
necessarily  primary  or  secondary. 
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In  the  most  complex  or  demanding  mission  phases,  the  two  task  cate¬ 
gories  may  require  all  of  the  pilot's  available  attention.  These  high 
workload  mission  phases  have  a  major  impact  in  design,  because,  as  tasks 
that  are  critical  for  either  control,  decision  making,  or  human  error 
potential,  they  provide  the  context  in  which  system  roles  are  established 
and  human  and  equipment  resources  are  allocated. 

The  managerial  tasks  often  result  in  discrete  action  sequences.  For 
many  of  these,  the  skilled  and  experienced  pilot  has  developed  a  nearly 
routine,  highly  rehearsed,  response  repertoire  to  meet  normal  and  many 
unusual  demands.  These  types  of  nearly  automatic  action  sequences  are 
subject  to  "slips"  of  intention  or  execution,  also  referred  to  as  "absent- 
minded  errors."  A  commonly-cited  example  of  a  slip  is  the  pilot's  failure 
to  lower  the  landing  gear  or  flaps  due  to  distractions  like  voice  communi¬ 
cations  and  in-cockpit  warning  alarms.  Current  studies  of  cognitive 
behavior,  associated  with  human  error  (e.g.,  Ref.  30),  emphasize  that 
slips  are  most  likely  to  occur  under  divided  attention  conditions. 

For  a  given  situation,  the  minimum  divided  attention  level  will  be 
established  by  the  control  tasks.  Consequently,  we  need  a  divided  atten¬ 
tion  model  for  control  operations.  The  model  should  provide  such  results 
as : 

•  The  nature  of  control  task  performance  degradation  due  to 

divided  attention  3^ 

•  An  indication  of  the  attentional  demands  required  for 
various  levels  of  control  activity  and  the  excess  capacity 
left  for  managerial  tasks. 

An  elementary  model  suitable  for  such  purposes  is  summarized  below.  It 
places  heavy  emphasis  on  both  the  attentional  demands  for  control  tasks 
and  the  excess  capacity  left  for  managerial  tasks.  These  are  quantitative 
indices.  The  attentional  demand,  r) ,  for  control  is  equal  to  the  average 
"control  dwell  fraction"  (0  <  r)  <  1)  ,  while  the  "excess  capacity"  left 
over  for  other  operations  is  the  average  "control  interrupt  fraction" 
(1-fj).  The  control  interrupt  fraction  is  therefore  also  termed  the  "mana¬ 
gerial  dwell  fraction." 
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The  theory  of  divided  attention  operations  can  be  considered  as  an 
extension  to  the  well-established  theory  of  display  scanning  and  signal 
sampling/reconstruction  (Refs.  1,  15,  and  31).  In  the  control  task,  the 
human  pilot's  behavior  can  again  be  characterized  in  mathematical  terms  by 
describing  functions  that  depend  on  the  effective  dynamics  of  the  aircraft 
being  controlled,  the  dynamics  of  the  pilot -vehicle  interfaces  (displays 
and  controllers),  and  a  "remnant."  These  two  components  are  depicted  in 
the  block  diagram  of  Fig.  10,  wherein  the  dynamics  of  the  effective 
pilot-vehicle  system  are  characterized  by  the  crossover  model  described 
previously.  Far  more  elaborate  models  of  the  pilot  can  be  used,  but  the 
crossover  model  is  quite  adequate  to  characterize  matters  at  the  level 
needed  here. 

Recall  that,  when  the  pilot's  full  attention  is  focused  on  the  con¬ 
trol  task,  the  crossover  frequency,  cjc,  of  the  pilot- vehicle  system  is 
maximized  consistent  with  near  minimum  mean- squared  error.  The  closed- 
loop  performance  issue  is  handled  by  a  minimization  process  that  arises 
from  a  compromise  in  following  the  command  input  while  reducing  the  rela¬ 
tive  influence  of  the  remnant.  The  remnant  in  full  attention  operations 
is  a  broadband  random  process  that  can  be  considered  as  a  pilot- induced 
noise . 

When  managerial  tasks  are  also  considered,  both  the  describing  func¬ 
tion  and  the  remnant  characterizing  the  pilot's  control  behavior  will  be 
affected  by  the  divided  attention  nature  of  the  pilot's  total  operations. 
The  describing  function  and  remnant  will  be  modified  to  account  for  the 
additional  signal  processing  or  supplementary  parallel  sensing  needed  to 
continue  control  operations  while  the  pilot  is  attending  to  the  managerial 
tasks.  Depending  on  the  specific  details,  these  modifications  may  reduce 
the  effective  pilot  gain,  add  to  the  effective  time  delay,  and/or  increase 
the  injected  noise.  Thus  the  system  crossover  frequency  will  be  reduced 
simultaneously  with  an  increased  contribution  of  noise  to  the  uncorrelated 
system  error.  Both  effects  will  cause  the  precision  of  control  task  per¬ 
formance  to  be  reduced  from  a  full  attention  baseline.  Similar  modifica¬ 
tions  to  the  pilot -vehicle  dynamics  are  made  even  with  full-attention 
control  operations  when  the  visual  cues  are  modified  to  call  for  divided 


Yc  =  Effective  dynamics  of  vehicle  (eg.,  aircraft  plus  stability 
augmentation  plus  displays) 

Yp  =  Full  attention  pilot  describing  function 

Yh  =  Perceptual  describing  function  to  account  for  divided  attention 

cu c  =  System  crossover  frequency 

re  =  Overall  pilot-vehicle  system  effective  latency 

c/)M  =  System  phase  margin  [v/Z  -  re u/c) 

<f>nn  =  Processing  remnant  spectrum  (n2~  e2 ) 
e2  =  Mean  squared  system  error 

Note:  eft)  is  "error"  and  subscript  "e"  in  re  is  "effective" 

Figure  10.  Pilot-Vehicle  System  for  Divided  Attention  Control  Task 
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visual  attention,  for  example,  in  changing  from  head-up  visual  meteoro¬ 
logical  conditions  (VMC)  to  the  head-down  instrument  panel  scanning  needed 
for  manual  approach  in  IMC  operations . 

In  the  divided  attention  situations  of  primary  interest,  it  is 
assumed  that  the  pilot  has  been  well  trained  in  the  control  and  managerial 
tasks  involved.  His  attention  is  allocated  among  control  and  managerial 
tasks  in  which  information  is  simultaneously  gathered  from  several  "per¬ 
ceptual  fields."  These  fields  may  include: 

Visual  "Segments" 

Foveal 

Parafoveal  •  Parallel  Pathways 

Peripheral 

Proprioceptive  "Segments" 

Vestibular 
Joint  receptors 

Stretch  receptors  •  Parallel  Pathways 

Pressure  receptors 
Etc . 

Aural  "Segments" 

Tactile  "Segments" 

and  others.  The  word  "segment"  is  intended  to  convey  the  properties  of 
extent,  thresholds,  input/output  dynamics,  etc.,  that  characterize  the 
particular  sensory  modalities  involved  as  they  are  integrated  into  useful 
perceptual  signal  sources.  The  easiest  to  describe  are  the  visual  percep¬ 
tual  field  segments,  which  can  be  divided  on  a  physiological  basis  into 
foveal,  parafoveal,  and  peripheral  pathways.  Besides  the  differing  spa¬ 
tial  (geometric)  extent  of  these  segments  there  are  also  differences  in 
threshold,  dynamic  properties,  contrast  background,  etc.  —  all  the  bewil¬ 
dering  complexities  associated  with  vision  in  its  myriad  details.  For  our 
purposes  here,  the  key  point  to  understand  is  that  a  visual  "display"  can 
be  attended  to  not  only  with  the  foveal  segment  but  also  with  the  para¬ 
foveal.  Thus  a  control  task  not  requiring  the  high  acuity  property  of 
foveal  vision  could  involve  sharing  between  the  foveal  and  parafoveal 
pathways  for  control,  with  attentional  adjustments  of  the  foveal  pathways 
between  the  control  operations  and  elsewhere  (e.g.,  reading  information, 
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conducting  visual  search,  etc.).  The  "perceptual  scanning"  process  in 
this  case  is  the  "switching"  of  the  input  signals  for  the  pilot's  control 
task  from  the  foveal  plus  parafoveal  to  the  parafoveal  alone  pathways. 

"Perceptual  scanning"  is,  of  course,  more  general  than  the  simple 
shifts  between  foveal  and  parafoveal  visual  pathways  serving  to  provide 
continuous  information  to  the  pilot  from  a  visual  display.  All  of  the 
other  perceptual  fields  for  each  input  modality  are  also  operating  more  or 
less  continuously  to  provide  signals  that  impinge  on  the  pilot's  sensor- 
ium.  Although  all  of  these  data  inputs  are  present,  they  are  not  neces¬ 
sarily  acted  upon  simultaneously.  However,  in  the  highly  trained,  unim¬ 
paired  pilot,  the  inputs  delivered  from  several  perceptual  fields  may  be, 
in  some  sense,  "operated  on"  in  parallel  all  of  the  time.  One  feature  of 
"impairment"  is  a  reduction  in  this  capacity  of  parallel  or  nearly  simul¬ 
taneous  operations  in  different  input  channels. 

For  example,  alcohol - induced  impairment  significantly  reduces  divided 
attention  capability.  The  operator  tends  to  focus  on  only  one  aspect  of  a 
task  at  a  time,  becoming  a  single  channel  device.  Physical  fatigue  and 
task-induced  stress  can  have  similar  effects,  although  these  are  more 
difficult  to  quantify  in  precise  terms  than  alcohol  effects. 

A  related  concept  needed  here  is  that  of  "attention,"  adding  to  the 
ability  to  sense  and  perceive  stimuli  a  readi ness  to  respond  to  selected 
stimul i .  By  analogy  with  visual  perception,  we  can  conceive  of  an  atten- 
tional  field  having  a  principal  focus  and  borders.  Attentional  fields 
have  both  spatial  and  intensive  aspects.  Thus  inattention  or  impaired 
attention  can  result  in  a  narrowing  of  the  spatial  borders,  an  increase  in 
the  minimum  stimulus  needed  to  cause  an  operator  output,  or  both.  A 
common  example  is  "tunneling"  of  vision  ("gunsight  vision")  wherein,  under 
highly  stressful  conditions,  the  visual  perceptual  field  is  narrowed.  As 
far  as  active  pilot  control  processes  are  concerned,  the  perceptual  scan¬ 
ning  and  attentional  field  features  are  joined;  that  is,  all  kinds  of 
perceptual  inputs  are  impinging  on  the  pilot  at  any  one  time,  but  the 
^  :tentional  foci  serve  to  activate  selected  perceptual  fields  as  sources 
of  control  or  managerial  task  "signals." 


41 


The  pilot's  primary  attention  may  be  shifted  from  one  signal  source 
to  another  in  the  course  of  conducting  a  particular  mission  phase.  Yet, 
when  a  control  task  is  involved,  it  must  be  attended  to  from  time  to  time. 
So,  too,  for  the  managerial  tasks.  In  the  course  of  operational  training, 
the  pilot  learns  to  switch  primary  attention  from  one  task  element  and/or 
perceptual  segment  to  another,  and  then  another,  and  back  to  the  first, 
etc.  This  is  conveniently  thought  of  as  a  perceptual  scanning  process. 
When  the  pilot  finally  becomes  skilled  in  the  operational  scenario,  the 
scanning  behavior  over  the  task  duration  exhibits  certain  stable  proper¬ 
ties  in  a  statistical  sense.  For  instance,  the  proportion  of  the  time 
spent  on  a  particular  input-gathering  chore,  the  dwell  times  on  certain 
instruments,  and  the  total  time  before  prominent  features  of  the  scanning 
process  are  repeated  tend  to  develop  stable  probabilities.  This  is  not  to 
say  that  the  scanning  is  either  periodic  or  uniformly  sequential  (i.e., 
from  "A"  to  "B"  to  "C"  and  back  every  time)  but  rather  that  cyclical 
activity  is  present  in  the  perceptual  scanning  process. 

Control  tasks  conducted  under  divided  attention  conditions  both  in 
flight  and  laboratory  research  have  shown  that  the  coverage  of  elements 
(e.g.,  instruments  or  perceptual  fields)  in  a  given  array  of  input  sources 
has  a  definite  average  frequency  and  corresponding  mean  sampling  interval, 
Ts ,  albeit  with  appreciable  variance.  The  mean  "control  dwell  time,"  T^ 
is  the  time  spent  on  information  sources  needed  for  control  purposes.  Its 
duration  depends  on  what  information  has  to  be  extracted.  The  ratio  of 
these  two  times  gives  the  "control  dwell  fraction,"  q  -  T(j/Ts ,  which  indi¬ 
cates,  on  the  average,  the  proportion  of  the  total  control  plus  managerial 
task  scanning  time  interval  required  by  the  control  task. 

The  information  transfer  characteristics  of  the  divided  attention 
attributes  of  the  human  controller  may  be  modeled  as  a  quasi - linear , 
random- input  "perceptual  describing  function,"  Y^.  This  multiplies  the 
full-attention  (continuous  control)  human  describing  function(s),  Yp,  to 
provide  the  describing  function(s)  for  the  human  pilot's  control  activi¬ 
ties  . 

The  simplest  way  to  develop  an  internal  signal  from  a  finite  duration 
sampled  input  is  to  act  proportionally  to  the  sampled  signal.  Then, 
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during  the  fixation  period,  T^,  the  pilot's  output  would  be  proportional 
to  the  perceptual  input  being  sampled,  while  outside  the  fixation  period, 
it  will  be  zero  (see  Fig.  11  lines  a  and  b) .  The  describing  function  is 
based  on  the  best  linear  fit  of  the  output,  in  the  mean- squared  sense. 
For  this  simple  finite  dwell  time  sampling,  the  perceptual  describing 
function  is  just  the  dwell  fraction  itself,  -  r; .  The  "remnant" 

accounts  for  all  of  the  pilot's  higher  frequency  power  not  linearly  con¬ 
nected  with  the  input.  The  describing  function  and  remnant  are  shown  on 
line  c  of  Fig.  11,  (Ref.  15).  [It  is  important  to  emphasize  that  the 
signals  shown  in  Fig.  11  are  highly  idealized  for  clarity.  Everything  is 
really  much  more  random:  the  signals  themselves,  the  dwell  times  (T<j)  , 
and  the  sampling  intervals  (Ts) . ] 

From  Fig.  11  it  is  easy  to  see,  as  the  divided  attention  level  is 
changed  to  reduce  the  control  dwell  fraction,  rj ,  that 

•  The  describing  function,  Y^,  is  reduced 

•  The  remnant  is  increased 

The  crossover  model  in  Fig.  11  shows  that  a  reduction  in  Y^  will  cause  a 
concomitant  reduction  in  the  pilot-vehicle  system  crossover  frequency,  tac . 
For  the  crossover  model,  uc  is  also  the  pilot -vehicle  system  loop  gain. 
Because  this  is  directly  related  to  the  system  phase  margin,  by 


where  re  is  the  overall  pilot-vehicle  system  latency,  then  the  reduction 
in  wc  will  be  reflected  in  increased  phase  margin. 

As  can  readily  be  appreciated  from  the  above  discussion,  the  effects 
of  divided  attention  can  have  profound  consequences  on  the  pilot-vehicle 
system  performance  in  control  activities.  These  can  be  conveniently  sum¬ 
marized  by  the  illustrative  case  sketched  in  Fig.  12.  As  already  noted, 
divided  attention  results  in  lower  crossover  frequency  and  associated 
increased  phase  margin.  As  far  as  the  pilot-vehicle  system  dynamics  are 
concerned,  a  major  consequence  is  a  significantly  increased  error  in  con¬ 
trol  activities.  As  shown  in  Fig.  12,  divided  attention  penalizes  the 
error  performance  in  two  ways: 
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•  By  reduction  of  the  permissible  crossover  gain,  and 

•  By  a  major  increase  in  the  remnant  due  to  the  divided 
attention  (i.e.,  lack  of  attention  to  the  control  tasks). 

Figure  12  shows  that  the  full  attention  pilot-vehicle  system  error 
begins  to  increase  only  as  the  dynamic  stability  limit  is  approached;  at 
lower  gains,  error  is  reduced  as  gain  increases.  While  a  similar  trend  is 
shown  for  divided  attention,  the  error  may  still  increase  without  bound 
for  circumstances  where  there  is  still  a  large  dynamic  stability  margin. 
This  is  because  the  closed- loop  effect  of  divided  attention  remnant,  the 
power  level  of  which  scales  with  mean-squared  error  as  in  Weber-law  noise, 
causes  error  signal  instability  in  the  mean-squared  sense  (Refs.  15  and 
32).  From  the  analyst's  point  of  view,  this  property  of  control  tasks 
with  divided  attention  requires  a  larger  phase  margin  (even  more  stable 
operation  of  the  control  task  than  with  full  attention)  as  the  control 
dwell  fraction  is  decreased. 

A  Webcr-law  model  of  divided  attention  remnant  has  been  applied  to 
the  error  signal  in  the  "crossover  model"  shown  in  Fig.  10  (Refs.  15,  33 
through  36).  The  model  of  divided  attention  remnant  includes  factors 
representing  average  attentional  dwell  time  fraction  (on  the  control  task) 
and  variability  thereabout.  A  quantitative  example  of  the  effects  of  divi¬ 
ded  attention  on  performance  is  presented  in  Fig.  13.  In  Fig.  13a,  the 
abscissa  is  normalized  crossover  frequency  (analogous  to  Fig.  12),  while 
Fig.  13b  provides  the  same  data  plotted  with  phase  margin  as  the  abscissa. 
The  forcing  function  is  white  noise  passed  through  a  third-order  Butter- 
worth  filter  with  normalized  breakpoint,  w^T-g  -  0.25.  The  full  attention 
condition  is  the  lowest  curve  in  both  portions  of  Fig.  13.  The  divided 
attention  conditions  that  govern  the  remnant  are  shown  as  families  with 
control  task  dwell  fraction,  ij ,  as  the  parameter.  In  this  example,  the 
normalized  control  dwell  interval  is  set  at  T^/re  -  1.5,  and  the  normal¬ 
ized  lower  bound  on  the  scanning  interval,  &  =  0.5. 

Figures  12  and  13  show  the  profound  effects  of  divided  attention  on 
control  system  performance  particularly  emphasizing  the  two  "stability 
limits."  The  first  is  the  full  attention  limit  given  by  wu  =  w/2 re,  which 
is  approached  by  the  full  attention,  r?  -  1 ,  curve.  The  second  is  the 
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"instability  in  the  mean  square."  This  is  associated  with  the  inequality 
constraint 

^2 

F  -  1  -  p  -  <  1  (12) 

2 

e 

The  bases  for  this  phenomenon  and  other  divided  attention  analytical  rela¬ 
tionships  will  be  summarized  below. 

Note  that  mean-square  error  instability  occurs  at  progressively 
increasing  phase  margins  as  the  attentional  dwell  fraction  on  the  control 
task  decreases.  Furthermore,  the  phase  margins  for  minima  in  normalized 
error  variance  are  even  greater,  and  the  minima  are  broad.  Typically,  the 
"blow  up"  phase  margin  is  less  than  the  phase  margin  for  best  performance 
by  10  to  16  deg.  Figure  14  puts  these  points  into  context  by  showing  the 
phase  margins  for  the  blow  up  condition  (F  -  1),  the  phase  margins  for  the 
minima  (from  Fig.  13),  and  the  phase  margins  for  a  value  of  error  coher¬ 
ence,  p|,  of  1/2  (corresponding  to  F  -  0.5).  The  latter  curve  nearly 
coincides  with  the  minimum  mean-squared  error  curve  when  the  control  dwell 
fraction  is  less  than  1/2.  For  larger  control  dwell  fractions,  say  from 
1/2  to  1,  the  phase  margin  for  minimum  mean-squared  error  is  essentially  a 
linear  function  of  dwell  fraction,  as  indicated  by  the  fit  on  Fig.  14. 

Analytical  formulas  (derived  in  Ref.  34),  on  which  constructions  such 
as  Figs.  13  and  14  are  based,  are  summarized  in  Table  4.  The  phase 
margin- dependent  function  K  ($M’  Te/^d)  [or  normalized  crossover 
frequency-dependent  function  A  (re  u>c ,  re/T^)  ]  is  shown  in  Fig.  15.  The 
curves  are  given  as  families  with  two  parameters:  (a)  the  normalized 
control  dwell  time  T d/re,  and  (b)  the  nondimensional  variable  /3re,  where 
fi  -  2/T(j-  The  T^/tg  -  0  (/3re  -  ■»)  curve  is  the  simplified  function  A(^)  . 
As  phase  margin  increases,  this  becomes  a  reasonable  approximate  bound  for 
the  more  complete  function. 

One  of  the  most  interesting  features  provided  by  the  formulas  is  the 
limit  associated  with  the  fundamental  constraint,  i.e.,  the 
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Figure  14.  Effect  of  Divided  Attention  on  Phase  Margins 
for  Minimum  Mean-Squared  Error 
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TABLE  4.  BASIC  DIVIDED  ATTENTION  RELATIONSHIPS 


DIVIDED  ATTENTION  REMNANT  POWER  SPECTRAL  OENSITY 


T  ( 1  -  n )(  1  -  <V  — r  0  <  S  -  T  /T  <  1 

«  (u)  i  -5 - 5 -  e2.  0  * 

/« T  vM  0  <  n  -  T./T  <  1 

4  +  ( -r)\ 

where  e^  &  fm  ♦  («)  du 

-'O  e  e 

Tq  la  Che  lower  bound  on  che  accendonal  scanning  or  saopllng  Interval 
Tfl  la  the  mean  value  of  the  accendonal  scanning  or  sampling  interval 
Td  la  the  mean  value  of  Che  accendonal  dwell  Interval 
♦ee(di)  Is  '■he  error  power  spectral  density,  (units  of  e  r  ror )  2/ (  rad /se  c  ) 
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where  G  Is  the  open-loop  describing  function  of  the  pilot  vehicle  system  (Fig.  /O) 
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(reu>c  ,  Td  /  re  ) 


Normalized  Crossover  Frequency,  wcTe(rad/sec  ) 


Figure  15a.  A  re/Td)  as  a  Function  of  Normalized  Crossover  Frequency 
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Phase  Margin,  <£M(deg) 


Figure  15b.  £  ^($M>  re/^d)  as  3  Function  of  Phase  Margin 
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Divided  Attention  Error  Stability  Limit: 


(13) 


The  curves  of  Fig.  15  can  be  used  directly  in  conjunction  with 
Eq.  (13)  to  determine  the  minimum  phase  margin  or  maximum  normalized 
crossover  frequency  available  for  a  given  level  of  divided  attention. 
The  maximum  value  of  A(^m,  T^/Tg)  must  be  less  than  re/(Ts  -  Te)(l  -  6). 
With  an  appropriate  change  of  labeling  on  the  ordinate,  the  curves  then 
become  boundaries  for  stability  in  the  mean-square,  with  locations  below 
the  curves  corresponding  to  allowable  phase  margins. 

For  some  purposes,  the  inequality  of  Eq.  13  may  be  awkward  to  work 
with  because  of  the  dependence  of  both  sides  on  T^.  The  simpler,  more 
approximate  form  using  the  A(^m)  may  therefore  be  more  useful.  With  this 
approximation,  the  Eq.  13  condition  becomes 


These  last  relationships  emphasize  the  need  to  constrain  the  system 
phase  margin  to  keep  the  error  in  divided  attention  operations  within 
bounds.  This  follows  because  (1  -r\)/r)  increases  as  the  managerial  demands 
increase.  [For  a  given  control  task,  the  overall  system  latency  is  the 
sum  of  the  net  high-frequency  system  lag  and  the  pilot's  effective  delay. 
The  control  task  dwell  time,  T^,  defines  how  long  the  pilot  must  fixate  on 
various  "display"  elements  to  assimilate  the  information  needed  for  con¬ 
trol.  Thus  T d/re  Is  approximately  constant  for  a  given  control-display 
task,  and  (1  -r\)/r\  governs  the  inequality].  Then,  as  the  maximum  allowable 
value  of  A(<£m>  rg/T^j)  is  reduced  to  maintain  the  inequality,  Fig.  13  indi¬ 
cates  that  the  divided  attention  control  task  phase  margin  must  be 
increased.  Because  the  normalized  crossover  frequency,  re  wc ,  is  directly 
’-elated  to  the  phase  margin  by  <j>fl  -  (tt/2)  -  re  wc ,  this  can  also  be 
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interpreted  as  indicating  that  the  control  task  crossover  frequency  is 
reduced . 


The  implications  of  these  statements  include: 

•  The  control  task  error  has  an  extremely  strong  dependence 
on  the  control  task  dwell  fraction.  (The  pilot-vehicle 
system  gain  is  reduced  and  the  system  "remnant"  or  effec¬ 
tive  uncorrelated  input  due  to  lack  of  attention  to  the 
control  task  is  increased  as  control  task  attention 
decreases) . 

•  If  the  task  complex  requires  significant  division  of  pilot 
attention  between  managerial  and  control  tasks,  the  dynam¬ 
ics  of  the  system  being  controlled  by  the  pilot  must  be 
able  to  support  very  large  pilot-vehicle  system  phase  mar¬ 
gins.  As  a  corollary,  the  controlled  system  must  possess 
dynamic  properties  that  require  little  attention  to  con¬ 
trol  . 

These  implications  are,  of  course,  consistent  with  the  conventional  wisdom 
that  attitude  control  and  path  control  functions  are  among  the  highest 
priorities  for  automation  as  means  to  reduce  pilot  control  workload. 
Steps  in  this  direction  cut  down  the  control  dwell  fraction  directly,  and 
increase  the  fraction  of  attention  that  can  be  devoted  to  managerial  task 
sequences . 
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SECTION  IV 


ALGORITHMIC  HUMAN  PILOT  MODEL 

A.  OVERVIEW  OF  THE  MODEL 

An  alternative  approach  to  the  estimation  and  description  of  human 
control  behavior  has  been  the  application  of  modern  optimal  control 
theory.  The  starting  points  in  this  process  are  the  well-founded  theory 
of  the  linear-quadratic-gaussian  stochastic  control  problem,  and  manual 
control  theory  and  data.  To  successfully  marry  these  two  elements  is  not 
easy,  yet  great  progress  has  been  made  (e.g.,  Refs.  37-47)  over  the  past 
twenty  years.  Some  notable  applications  to  flying  qualities  problems  have 
also  been  published  (e.g.,  Refs.  48-50).  The  concept  rests  on  the  pre¬ 
sumption  that  human  operator  responses  can  be  emulated  by  an  analogous 
optimal  control  system.  The  optimal  system  operates  to  minimize  a  qua¬ 
dratic  performance  index  in  the  presence  of  various  system  inputs  and 
noises.  In  doing  so  it  provides  a  representation  for  at  least  some  of  the 
adaptive  characteristics  of  the  human  operator.  The  basic  consideration  in 
this  algorithmic  approach  is  the  provision  of  techniques  for  imposing 
those  characteristics  of  the  human  which  represent  both  favorable  (e.g., 
adaptation)  and  unfavorable  (e.g.,  time  delay  and  remnant)  features  so 
they  are  consonant  with  experiment.  Related  techniques  must  account  for 
certain  very  fundamental  human  characteristics,  such  as  the  effective  time 
delay  and  neuromuscular  delays. 

A  general  picture  of  the  algorithmic  or  optimal  control  (OCM)  model 
is  shown  in  Fig.  16.  At  the  top  are  the  machine  properties  involving  the 
controlled  element  and  display  as  acted  on  by  disturbances.  These  are 
represented  by  linear  state  vector  and  display  vector-matrix  equations. 
The  disturbance,  v(t) ,  is  a  vector  of  white  gaussian  noise  processes. 
This  represents  both  forcing  functions  and  disturbances  impacting  the 
controlled  element.  If  the  forcing  functions  are  colored,  they  are 
represented  by  filtered  white  gaussian  noise.  The  additional  states 
required  to  represent  the  filter  dynamics  are  appended  to  the  controlled 
element  state  vector  and  result  in  expanded  A,  B,  C,  and  E  matrices. 
Deterministic  disturbances  can  be  modeled  by  adding  non-zero  mean 
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components  to  the  disturbance  vector,  with  the  addition  of  still  more 
elements  to  the  state  vector  and  associated  matrices.  The  display  vari¬ 
ables  are  linear  combinations  of  the  system  states  and  the  pilot  output. 

In  the  optimal  control  formulation  the  human  pilot's  characteristics 
can  be  divided  into  two  categories  --  those  which  represent  intrinsic 
human  limitations,  and  thus  which  are  not  subject  to  optimization,  and 
those  properties  which  are  subject  to  adaptation  and  thus  optimization. 
In  the  first  category  are  the  effective  time  delay  and  the  remnant.  To 
some  extent  the  neuromuscular  system  properties  and/or  the  pilot- vehicle 
system  crossover  form  and  bandwidth  also  fall  into  this  category,  although 
their  connections  in  the  OCM  formulation  are  somewhat  obscure.  (This  will 
be  illustrated  later.)  In  the  optimal  control  model  the  remnant  is 
accounted  for  by  observation  noise  and  motor  noise,  shown  respectively  at 
the  pilot's  input  and  neuromuscular  command  output  points.  The  observation 
noise  vector  is  added  to  the  display  output  y(t) .  A  separate  noise  com¬ 
ponent,  Vy^  (t)  is  associated  with  each  display  output  component,  y^(t). 
As  noted  in  Fig.  7a,  the  remnant  added  at  the  operator's  input  is  rela¬ 
tively  wideband,  so  each  component  is  assumed  to  be  an  independent  gaus- 
sian  white  noise  process.  The  spectral  density  is  proportional  to  the 
mean- squared  value  of  the  displayed  component,  with  a  proportionality 
factor  Py^ ,  which  is  a  noise- to-signal  ratio.  In  general,  the  human 
operator  is  presumed  to  obtain  both  displacement  and  rate  information  from 
a  single  display  variable,  and  good  results  have  been  obtained  by  assuming 
that  Py^  for  the  position  and  rate  variables  is  the  same.  In  single-loop 
control  situations  numerical  values  of  Py^  of  about  0.01  are  typical.  As 
can  be  appreciated  from  Fig.  7a,  this  is  relatively  invariant  over  a  wide 
range  of  system  dynamics  and  input  spectra.  To  the  extent  that  this  is 
so,  the  normalized  observation  noise  can  be  considered  to  be  primarily 
operator-dependent . 

The  many  internal  time  delays  associated  with  visual,  central  proces¬ 
sing,  integrative,  and  other  operations  are  combined  into  a  lumped  percep¬ 
tual  delay,  t.  For  simplicity,  it  is  assumed  in  the  current  optimal  con¬ 
trol  model  that  all  outputs  are  delayed  by  the  same  amount.  (As  noted 
previously  in  connection  with  the  structural  isomorphic  model,  there  is  a 
delay  increment  associated  with  rate  perception.) 
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The  "motor  noise,"  like  the  observation  noise,  is  assumed  to  be  a 

zero-mean  gaussian  white  noise  with  spectral  density  proportional  to  the 

mean-squared  operator  output.  An  additional  component,  Vu  ,  is  sometimes 

aQ 

included  to  account  for  the  fact  that  the  human  operator  introduces  noise 

into  an  undisturbed  system.  A  motor  noise/signal  ratio,  pn  ,  of  0.003 

ai 

has  been  found  to  provide  a  good  match  to  some  experimental  data. 


The  neuromuscular  system  is  represented  by  a  lag  matrix,  Tjj.  This 
is  not  explicitly  modeled  as  an  inherent  limitation.  Instead,  it  is 
imposed  by  weighting  control  rate  terms  in  the  cost  function  used  to 
generate  the  optimal  control.  For  single -loop  control  problems  with 
linear,  wide  bandwidth  manipulators,  this  weighting  is  purposely  selected 
to  yield  T^  of  approximately  0.1  sec  to  represent  this  inherent  limita¬ 
tion.  As  will  be  seen  later,  this  weighting  tends  to  set  the  frequency 
range  over  which  the  pilot-vehicle  system  may  approximate  the  crossover 
model  in  a  single  input,  single  output  system.  When  everything  is  taken 
into  account  in  an  effective  pilot  describing  function  Yp  the  direct 
neuromuscular  lag  represented  by  T^  will  be  cancelled  by  other  quantities, 
although  the  total  effective  time  delay  may  reflect  some  neuromuscular 
lag. 


The  remaining  elements  of  the  human  operator  are  adaptive  to  the 
system  characteristics  and  to  changes  in  the  explicit  human  operator  limi¬ 
tations  described  above.  Estimation  of  the  delayed  state  vector  is  accom¬ 
plished  via  a  Kalman  filter.  This  delayed  state  estimate  is  fed  to  a 
least-mean-squared  predictor  to  yield  the  estimated  state  vector,  x(t). 
The  optimal  gain  matrix,  L,  is  generated  by  solving  the  optimal  regulator 
problem  for  a  quadratic  cost  function  of  the  form 


J(u)  =*  E  x  J  (y'Qy  +  u'Ru  +  u'Gu)  dtj’ 


(15) 


Because  the  cost  functional  weightings  preordain  the  details  of  the  con¬ 
troller  gain  matrix,  L,  the  selection  of  weightings  is  critical  to  the 
model's  success.  This  is  particularly  the  case  when  the  model's  purpose 
is  to  simulate  human  operator  responses.  For  simple  single- loop  control 
situations,  excellent  agreement  with  experimental  measurements  has  been 
obtained  with  a  cost  functional  of  the  extremely  simple  form: 
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(16) 


J(u)  -  E  |T1imoo  |  (e2  +  Gc2)  dt 

•  • 

where  e  is  the  compensatory  system  error  and  c  -  u  is  the  operator's  con¬ 
trol  rate.  The  value  of  G  is  selected  as  described  above  to  yield  an 
appropriate  neuromuscular  delay,  T^.  For  more  complex  situations,  the 
relative  weights  are  determined  based  either  on  maximum  allowable  devia¬ 
tions  or  limits,  or  from  a  knowledge  of  human  preferences  and  capabili¬ 
ties.  This  is  similar  to  the  technique  suggested  by  Ref.  51,  wherein  the 
weighting  on  each  quadratic  term  is  simply  the  inverse  of  the  square  of 
the  corresponding  allowable  deviation.  The  solutions  for  this  modified 
Kalman  filtering  prediction  and  optimal  control  problem  are  given  by,  for 
example,  Refs.  38,  39,  42,  47,  and  52. 

Because  the  power  spectral  densities  of  the  observation  and  motor 
noises  scale  with  mean-squared  values  of  variables  within  the  system,  an 
iterative  solution  method  is  required  to  achieve  a  solution  with  pre¬ 
scribed  noise  ratios.  This,  and  other  complications ,  make  application  of 
the  OCM  dependent  on  computer  implementation.  Although  open  literature 
implementations  exist,  notably  PIREP  (Ref.  47),  they  are  not  readily 
available  or  suitable  for  PC  style  calculations.  Here  we  will  describe  an 
implementation  of  the  OCM  using  Program  CC  (Ref.  53),  which  is  anticipated 
to  greatly  increase  the  availability  and  understanding  of  the  OCM. 

Transfer  functions  of  the  operat'  ;  dynamic  characteristics  obtained 
using  the  OCM  are  intrinsically  very  high  order.  For  example,  for  single 
loop  systems  as  treated  later  the  human  transfer  function  denominator  is 
of  order  2n  +  5,  where  n  is  the  order  of  the  system  which  models  the 
controlled  element  plus  the  driving  noise.  This  contrasts  markedly  with 
such  simple  descriptions  of  human  dynamic  behavior  as  the  crossover  model 
(Ref.  1).  Both  the  crossover  model  and  the  OCM  can  descrit  experimental 
data  fairly  well,  so  it  follows  that  OCM  results  should  be  susceptible  of 
considerable  simplification.  This  possibility  will  be  explored  below. 
The  approach  to  simplification  is  to  include  in  the  Program  CC  implementa¬ 
tion  of  the  OCM  a  novel  transfer  function  approximation  of  the  OCM's  pilot 
describing  function.  This  is  expected  to  greatly  enhance  understanding 
and  interpretation  of  OCM  results 
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Additional  features  and  numerous  applications  have  appeared  in  the 
years  following  the  OCM's  introduction,  notably  (Refs.  54-56).  The  addi¬ 
tional  features  included  in  Program  CC's  implementation  are  visual  indif¬ 
ference  thresholds  and  fractional  attention  parameters.  Notably  absent 
is  the  use  of  pseudo-noise  to  induce  low  frequency  phase  droop,  and  the 
optimization  of  fractional  attention  for  multi- input  problems.  Article  B 
will  present  a  single  axis  version  of  the  OCM,  including  a  description  of 
time  and  frequency  domain  performance  measures  and  a  novel  transfer  func¬ 
tion  approximation  of  the  algorithmic  model's  pilot  describing  function. 
To  later  handle  divided  attention  situations  attentional  fractions  are 
incorporated  in  the  model.  These  permit  a  relatively  simple  extension  to 
cover  the  multi-axis  case,  which  is  done  at  the  end  of  the  article. 
Article  C  provides  an  outline  of  the  implementation  of  the  OCM  using 
Program  CC.  This  is  more  thoroughly  treated  in  Volume  III.  In  Article  D 
the  OCM  is  exercised  in  a  single  axis  example  which  is  used  to  compare 
with  the  crossover  model. 

B.  DETAILS  OF  THE  HUMAN  OCM 

1.  Single  Axis  Problem  Statement 

In  full  attention  operations  the  human  operator  manually  controls  a 
system  with  the  objective  of  minimizing  the  system  error  in  some  sense. 
For  a  simple  mathematical  formulation  we  will  presume  that  imprecisely 
stated  objective  will  be  satisfied  by  minimizing  a  mean  square  error. 
This  is  the  situation  for  many  standard  tracking  tasks.  A  state  space 
model  of  the  controlled  system  has  n  states,  1  input,  2  outputs,  and  2 
observations : 

x  -  Ax  +  Bu  +  Ew  (17) 

y  -  Cx  +  Du  (18) 

yp  -  y( t  -  r)  +  vy  (t  -  r )  (19) 

The  2  outputs  are  error  ye  and  error  rate  <fe ,  and  the  observations  are 
delayed  by  the  human's  assumed  visual  delay  of  r  -  .2  seconds.  The  source 
of  the  error  is  modelled  by  a  white  noise  source  w(t)  with  intensity  Vw , 
passed  through  a  shaping  filter  which  has  its  dynamics  included  in  the  A 


The  observation  noise  is  modelled  by  a  white  noise  source  Vy(t)  with 
intensity  Vy  -  diag(vu^  vy2^-  where: 


The  observation  noise  ratio  is  py ^  -  .01.  The  attentional  fraction  is 
0  <  f  <  1,  with  the  remainder  devoted  either  to  other  control  tasks  or  to 
non-control  tasks  such  as  communication.  (In  the  crossover  model  for 
divided  attention  the  attentional  fraction  was  devoted  as  rj ,  which  is 
equivalent  to  the  f  used  here.  The  notational  change  is  made  to  keep 
these  developments  consistent  with  the  literature  from  which  they  were 
derived.)  No  control  action  results  from  errors  below  the  indifference 
threshold  T^,  which  is  based  on  human  eye  perceptual  levels  (e.g.,  .05° 
and  ,l°/sec).  The  mean  square  error  of  the  output  is  ay^<  which  is 
increased  due  to  the  indifference  threshold  via  the  erfc  function,  which 


ranges  from  0  to  1 : 

E(<7y^ i  Ti) 


B) 


Prob  (|yi|  >  Ti) 


It  is  important  to  note  that  Vy^  is  proportional  to  the  mean  square  error, 
which  makes  an  iterative  KBF/linear  predictor  solution  necessary. 

According  to  the  human  OCM,  the  human  operator  exerts  control  actions 
which  minimize  the  performance  index: 

J  ‘  T1!*.  {  f  Jo  ♦  G“2)  «  }  (22) 

The  control  action  is  a  linear  combination  of  the  estimated  state  of  the 
system.  The  linear  combination  is  determined  via  a  Linear  Quadratic  Regu¬ 
lator  (LQR)  problem,  and  the  state  estimation  is  determined  via  a  Kalman 
Bucy  Filter  (KBF)  and  a  linear  predictor. 

The  LQR  problem  is  solved  by  augmenting  the  input  with  an  integrator 
and  minimizing  an  equivalent  performance  index  (using  p  -  &) : 

^o  ”  A0x0  +  B0p  (23) 
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A  B 


xo 


The  solution  is: 


Ao 


0  0 


M  “  '[Ll  1-2 1  xo 
L  -  G' 1BqK0 

0  -  AqK0  +  K0A0  +  Q0  -  KqBqG^B^Ko 


(24) 

(25) 

(26) 
(27) 


The  weight  G  on  the  control  rate  &  is  chosen  in  an  iterative  fashion  so 
that  L2^  -  Tjj  -  .1,  i.e.,  so  that  the  state  feedback  around  the  integra¬ 
tor  models  the  neuro-muscular  mode  with  a  time  constant  of  .1  seconds. 

The  KBF/linear  predictor  problem  estimates  the  state  of  the  con¬ 
trolled  system  augmented  with  the  neuromuscular  mode  and  noise  driven  by 
an  additional  noise  source: 


*1  ”  Alxl  +  Blua  +  wi;  y  - 


> 

X 

> 

W 

A  B  ’ 

xp  - 

;  wL  = 

;  A!  - 

-vua. 

0  'L2. 

W1 


0  1 
Ll2 


;  Ci  -  (C  D) 


evwe  ' 


L2vuaL2 


(28) 

(29) 


(30) 


(31) 


The  white  noise  source  w^(t)  has  intensity  Wj  as  shown  above  and  incor¬ 
porates  the  motor  noise  vu  (t)  ,  which  is  inserted  before  the  neuro- 

a 

muscular  mode,  and  has  its  intensity  defined  as  a  ratio: 


(32) 


with  the  noise  ratio  pu  =  .01  or  .003.  It  is  not  physically  possible  to 

a 

experimentally  verify  this  motor  noise  model,  but  it  is  commonly  accepted 
for  reasons  of  mathematical  tractability . 

The  KBF  computes  p(t),  the  linear  mean-square  estimate  of  x^(t  -  r), 
based  on  observations  yp (a)  for  a  <  t . 
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p(t)  -  Aip(t)  +  Hx[yp(t)  -  Cxp(t) ]  +  BXua(t) 
Hi  -  SiCiVy-1 

0  -  Ai^i  +  £XAi  +  Wi  —  EiCiVy'^CiEi 


(33) 

(34) 

(35) 


The  linear  predictor  updates  p(t)  to  obtain  xx(t),  the  linear  mean- 
square  estimate  of  xi(t)  based  on  observations  yp(a)  for  a  <  t.  Note  that 
p(t  -  r)  *  xX(t) . 

|(t)  -  A^(t)  +  BXua(t)  (36) 
xX  -  ?(t)  +  eAir  (p(t)  -  i(t  -  r ) ]  (37) 


A 

The  LQR  control  weights  are  applied  to  xX(t),  but  only  to  the  portion 
which  estimates  the  system  states: 

ua(t)  -  -L*xx(t);  L*  -  (L2‘1L1  0)  (38) 

A  block  diagram  of  the  human  optimal  control  model  is  presented  in 
Fig.  17. 

2.  Time  Domain  Performance 

Time  domain  performance  measures  are  the  mean  square  errors  and  the 


optimal  i 

cost.  Define 

A 

the  estimation  error  ex  -  xx  —  xx ,  then: 

Ei  - 

E(  exei )  - 

eAlTZieAlr  +  T  eAla  WleAla  da 

J  0 

(39) 

A 

Xi  = 

A  A  ( 

E(xxxi )  = 

£  e(Al  -  eAlrH1VyHieAiTe(Al  ~  BlL*>'d* 

(40) 

Xl  - 

E{xxxi )  = 

El  +  Xl 

(41) 

Y  - 

Elyy'  ) 

ClXiCi 

(42) 

ua  ~ 

E(uaua}  = 

L*XXL* ' 

(43) 

The  KBF/linear  predictor  solution  is  iterative  due  to  the  definition 
of  these  noise  intensities.  Use  the  mean  square  value  Y  and  Ua  to  update 
the  noise  intensities,  and  stop  the  iteration  when  the  achieved  noise 
ratios  are  within  about  .1  power  dB  of  the  desired  values  (i.e.,  within 
■01  of  20  logio [ p ] ) • 
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PREDICTOR-  CONTROLLER 


Strictly  speaking,  the  performance  index  J  is  infinite,  because  the 
motor  noise  feeds  directly  into  Circumvent  this  problem  by  using  u 

instead.  The  optimal  cost  is: 

•  •  •  A 

U  -  E{uu'}  -  LXj^L'  +  (0  L2)E1(0  L2)'  (44) 

J  -  Tr(YQ)  +  Tr(UG)  (45) 

3.  Frequency  Domain  Performance 

Classical  studies  of  human  operator  models  use  transfer  function 
models.  The  controlled  system  has  just  the  error  as  its  output,  i.e., 
Yc(s)  -  ye/u,  and  the  pilot  describing  function  is  defined  just  from  the 
error,  i.e.,  Yp(s)  -  (u/ye)  +  s(u/£e).  Frequency  domain  performance 


measures  are  based 

on  Yp  and  YpYc ,  for  example 

“c 

-  system  crossover  frequency 

(46) 

<f>H 

-  phase  margin 

(47) 

Te 

=  (m/2  -  <£pm)/wc  -  effective  delay 

(48) 

^pc 

=  Arg[ Yp( jwc) ]  -  wcre  -  arctan[ucrN] 
pilot  phase  compensation 

(49) 

Due  to  the  observation  delay  the  optimal  solution  incorporates  a 
linear  predictor.  It  is  therefore  not  possible  to  express  Yp  as  a  simple 
transfer  function,  and  this  has  been  a  major  impediment  to  understanding 
the  OCM  from  a  classical  point  of  view. 

The  following  new  result  gives  an  approximate  transfer  function  for 
Yp.  Use  a  Padt  pproximation  for  the  delay,  where  Up(t)  =  ua  (t  -  r): 

&p(t)  =  ApXp(t)  +  BpUa( t)  (50) 

up ( t )  ”  CpXp(t)  +  DpUa(t)  (51) 

It  then  follows  that  the  combined  KBF/linear  predictor/neuro-muscular 
system  is  approximated  by: 

*2  =  A2x2  +  B2(y  +  vy)  +  E2vu 


A  ^ 


u  -  C2x2 


(52) 


Ai  -  HiCi  -BiDpL*  BiCp  0  1 

-e“AlTH1C1  Ax  -  B]L*  0  0 


a2  “ 

0 

-BpL* 

AP 

0 

L 

0 

-L2DpL*  L2Cp  - 

-l2. 

f 

P 

'  o' 

X1 

eAlTHl 

0 

x2  “ 

Xp 

B2  - 

0 

e2  - 

0 

.  u. 

.  0  . 

L2. 

(53) 


(54) 


C2  -  (0  0  0  1) 


(55) 


The  derivation  is  tedious  but  straightforward  (replace  f  with  ft^) •  A 
state  space  realization  for  Yp  is: 

*2  “  a2x2  +  (b2i  +  A2B22)Ye  (56) 

u  -  C2x2  +  C2B2zye  (57) 

The  2nd  subscript  indicates  columns  of  the  input  matrix,  and  is  based  on 
the  state  space  derivative: 

S7(sl  -  a)'^  /3  -  7/3  +  7 ( s I  -  a)'l  a/3  (58) 


The  state  space  realization  for  Yp  can  be  converted  to  a  transfer  function 
by  standard  techniques  such  as  Fadeeva's  method  or  eigenvalue  and  gener¬ 
alized  eigenvalue  problems.  The  order  of  Yp(s)  will  be  2n  +  5,  and  will 
contain  many  exact  and  approximate  pole/zero  cancellations. 

4.  Multi-Axis  Extension 

The  human  OCM  formulation  naturally  extends  to  the  multi- -axis  case. 
Here  we  are  concerned  with  the  case  of  multiple  independent  axes.  The 
human  operator  scans  several  error  indicators,  each  of  which  is  controlled 
by  only  one  input.  Further,  the  operator  is  instructed  to  minimize  the 
errors  in  all  axes  and  not  to  give  preference.  These  conditions  can  be 
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forced  experimentally,  and  are  a  very  good  approximation  of  a  pilot  separ¬ 
ately  controlling  longitudinal  and  lateral  motion. 


The  assumptions  of  the  human  OCM  problem  for  these  conditions  are  (1) 
the  single  axis  tasks  are  separately  optimized,  and  (2)  the  attentional 

fractions  are  allocated  so  as  to  minimize  the  following  total  cost: 

N 

"axes 


^total 


± _  T 

2  Jaxis 
oc . 
ci 


<fi) 


(59) 


The  axes  are  normalized  with  respect  to  the  mean  square  command  oc^, 
defined  to  be  due  to  the  driving  noise.  Note  that  these  assumptions 
allow  the  operator  to  pay  more  attention  to  an  axis  which  is  difficult  to 
control.  The  attention  is  apportioned  across  axes,  so  the  same  attention 
within  an  axis  is  paid  to  the  error  and  error  rate.  The  attentional  sum 
across  axes  is  <  1,  depending  on  whether  a  non-control  task  is  required. 


The  multi-axis  problem  can  be  set  up  and  solved  as  a  multivariable 
problem.  PIREP,  one  of  the  few  currently  available  implementations  of  the 
OCM,  solves  the  multivariable  problem  for  given  set  of  attentional  frac¬ 
tions,  and  then  optimizes  over  the  fj/s  using  a  gradient  technique  which 
iteratively  computes  at  each  step  Jtotal  ancl  an  estimate  for  3f/3Jtota^. 
PIREP' s  implementation  is  very  involved,  and  is  not  suitable  for  personal 
computers.  The  multi-axis  example  of  Appendix  A  is  based  on  PIREP. 

A  different  solution  to  the  multi-axis  problem  which  requires  far 
fewer  calculations  is  presented  here.  This  approach  is  exemplified  in 
Appendix  B  using  the  same  data  as  Appendix  A.  The  single  axis  problems 
are  separately  solved,  and  then  the  final  optimization  over  the  f^'s  is 
computed  by  hand.  The  final  optimization  is  based  on  the  empirical  obser¬ 
vation  that: 

^axis^  ~  ai  +  ^i  f°r  ”  ^■•■•••^axis  (60) 

n  _  1 


This  approximation  has  been  observed  to  hold  (for  the  experimental  systems 
used  by  Dander  (Ref.  57)  in  a  range  varying  from  f^  >  .5  to  f^  >  .1, 
depending  on  the  amount  of  lead  required  by  the  operator.  Each  single 


axis  problem  is  solved  for  enough  values  of  f^  to  be  able  to  estimate  the 
a^  and  b^  coefficients  and  to  determine  the  range  of  validity.  The  final 
optimization  over  the  fj/s  can  then  be  explicitly  solved  for  the  important 
cases  of  2  and  3  axes. 


The  problem  is  solved  using  standard  Lagrangian/Hamiltonian  techniques. 
For  2  axes: 

1/fl  -  1  +  A27a!  (63) 

1/f  2  “  1  +  •/al/a2  (64) 

and  for  3  axes: 

1/fl  -  1  +  J&2/a\  +  7a 3/al 
l/f2  -  1  +  Ja\/ a2  +  Ja 3/a2 

l/f3  =  1  +  J al/a  3  +  7a2/a3 

A  final  check  should  be  made  to  see  if  the  fj/s  lie  in  range  for  which  the 
approximations  using  a^  and  b^  are  valid.  If  not,  a  gradient  search  tech¬ 
nique  using  piecewise  linear  approximations  can  be  implemented,  but  we 
have  yet  to  find  this  necessary. 

C .  IMPLEMENTATION 

The  human  OCM  is  a  variant  of  the  standard  Linear  Quadratic  Gaussian 
(LQG)  problem,  with  the  following  distinctions.  The  algorithm  for  solving 
the  human  OCM  is  summarized  in  Fig.  18. 


(65) 

(66) 
(67) 
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Figure  18.  Iterative  Algorithm  for  Solving  the  OCM 
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•  The  performance  index  weights  6,  not  u.  The  LQR  problem  is 
solved  by  augmenting  the  system  with  an  integrator,  and  the 
feedback  around  the  integrator  models  a  neuro-muscular 
mode . 

•  The  observations  are  delayed,  hence  the  need  for  a  linear 
predictor . 

•  The  KBF  uses  the  controlled  system  augmented  with  the 
neuro-muscular  mode,  not  the  system  alone  or  the  system 
with  the  integrator. 

•  The  observation  and  motor  noise  intensities  are  defined  as 
a  ratio  of  the  respective  mean  square  outputs  and  inputs. 

An  iteration  is  required  to  achieve  this  ratio. 

•  An  additional  iteration  is  required  to  minimize  the  multi¬ 
axis  problem  with  respect  to  the  attentional  fractions. 

Program  CC  contains  all  of  the  required  low  and  high  level  routines 
for  solving  the  OCM  in  a  command  driven,  user  friendly  environment.  It  is 
significant  that  Program  CC  is  implemented  on  a  personal  computer,  thus 
greatly  increasing  the  availability  of  the  OCM.  Program  CC  is  used  for 
design  and  analysis  of  a  wide  range  of  classical,  digital,  sampled-data , 
state  space,  and  optimal  control  systems.  The  state  space  optimal  control 
part  of  the  package  are  used  to  solve  the  OCM  problem,  the  state  space 
realization  for  Yp  can  be  converted  to  a  transfer  function,  and  then  Bode 
plots  (also  Nichols,  Nyquist,  root  loci,  and  Bode-root-loci  plots)  can  be 
used  to  analyze  Yp  just  like  any  other  classical  compensator.  The  inter¬ 
active  graphics  makes  the  classical  analysis  particularly  insightful. 

The  low  level  numerical  routines  which  are  required  to  solve  the  OCM 

are  : 

1.  Block  matrix  manipulation 

2.  Riccati  equation 

3.  Lyapunov  equation 

4.  Matrix  exponential 

5.  erfc  function 

6.  iteration  confrol 

Elaborating  further:  (1)  Block  matrix  manipulation  is  symbolically 
implemented,  so  that  equations  can  be  entered  and  solved  in  a  very  natural 
fashion.  State  space  systems  or  quadruples  are  treated  as  a  single  data 
element,  which  greatly  decreases  bookkeeping.  (2)  Riccati  equations  are 
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used  to  solve  the  LQR  and  KBF  problems.  Potter's  method  is  used,  which 
computes  the  stable  invariant  subspace  of  a  Hamiltonian  matrix.  Schur 
vectors  are  used  for  the  required  matrix  decomposition,  which  is  numer¬ 
ically  better  conditioned  then  the  more  standard  eigenvector  methods. 

(3)  Lyapunov  equations  are  used  to  compute  mean  square  errors.  A  variant 
of  the  Bartel-Stewart  method  is  used.  The  Bartel-Stewart  method  decom¬ 
poses  the  linear  Lyapunov  matrix  equation  into  a  form  which  can  be  recur¬ 
sively  solved.  The  variant  uses  Schur  vectors  instead  of  eigenvectors. 

(4)  The  matrix  exponential  is  needed  as  part  of  the  linear  predictor.  The 
matrix  exponent  is  first  scaled,  a  Pade  approximation  is  computed,  and 
then  the  result  is  unsealed.  More  information  on  matrix  algorithms  is 
available  from  Thompson  (Ref.  53)  and  Golub  and  Van  Loan  (Ref.  58). 

The  many  commands  needed  for  the  solution  are  combined  into  several 
macros  and  one  user-def ined-command.  A  macro  is  a  sequence  of  commands 
which  allows  nesting  and  parameter  insertion.  A  user-def ined-command  is 
a  separate  program  callable  from  within  Program  CC,  which  is  this  case 
computes  the  erfc  function  and  controls  the  LQR  and  KBF  iterations. 

There  are  many  parameters  in  the  OCM  which  define  human  behavior,  yet 
there  are  generally  accepted  values  for  each.  It  is  easiest  to  use  the 
default  values,  but  any  other  values  may  be  substituted.  Considerable 
work  may  be  needed  to  characterize  the  controlled  system  and  the  driving 
noise,  but  the  solution  of  the  OCM  is  for  the  most  part  a  turn- the -crank 
type  operation. 

The  OCM  problem  stretches  the  limits  of  personal  computers,  even  for 
relatively  small  order  problems.  There  is  no  problem  with  memory  size, 
disk  size,  and  numerical  accuracy;  but  there  realistically  is  a  problem 
with  computation  time,  becoming  quite  noticeable  for  systems  above  15^ 
order,  with  the  number  of  LQR  and  KBF  iterations  being  the  bottleneck. 
The  number  of  iterations  can  be  reduced  by  using  good  starting  value  for 
the  noise  intensities,  which  is  feasible  if  the  same  or  similar  problems 
have  been  solved.  Low  order  effective  models  of  the  controlled  system  can 
be  used  for  most  flying  qualities  problems. 
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Solving  the  multi -axis  as  several  single  axis  problems  requires  fewer 
total  calculations  (i.e.,  flops),  as  now  conservatively  estimated. 
Several  of  the  individual  steps  such  as  Riccati  solutions  are  n-^  opera¬ 
tions,  so  call  the  entire  OCM  solution  an  n^  operation  with  a  large  multi¬ 
ple.  If  the  3  axis  multivariable  problem  is  solved  in  one  step  this  takes 
27n^  operations,  and  if  there  are  5  iterations  for  the  f^  gradient  search 
routine  then  the  total  is  135n^.  On  the  other  hand  if  the  single  axis 
problems  are  solved  this  requires  3n^  operations,  and  each  axis  must  be 
solved  say  5  times  for  different  fj/s  for  a  total  of  15n^  operations;  a 
nine -fold  savings . 

D.  AN  EXAMPLE  FOR  Yc  =  Kc/s 

To  give  an  idea  of  the  procedure  and  some  of  the  results  that  can  be 
obtained,  an  example  with  a  very  simple  plant  comprising  a  pure  integrator 
is  presented.  Although  this  is  a  first-order  approximation  to  many  prac¬ 
tical  manual  control  systems,  such  as  aircraft  roll  control,  it  is 
selected  here  to  permit  a  fine-grained  focus  on  the  OCM  version  of  the 
operator's  dynamics  without  distortion  of  the  numerical  results  from 
higher  order  controlled  element  complications.  The  controlled  system  and 
driving  noise  are: 

ye(s)  =  Yc(s)  ( u ( s )  +  Yw(s)w(s) ]  (68) 

Yc(s)  “  p  Yw(s)  =  (69) 

Note  that  the  driving  noise  is  added  to  the  input  of  Yc .  Use  the  follow¬ 
ing  parameters  (as  in  Ref.  44): 

W  T  TN  Pyi  Py2  Pya  f  Ti  l2 

8  8  .15  .1  .01  .01  .003  100 

The  command  and  macro  calls  for  Program  CC  are: 

COg-l/s  enter  Yc 

CC>gl-3/(s+2 )  enter  Yw 

CC>@ocmyc2 , g, gl ,p40  create  s.s.  p^p  “  Yc 

CC><3ocmal  1 ,  p40  solve  OCM 

CC><<;ocmg  ,  g2  compute  g2  -  Yp 
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The  time  domain  results  are: 
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.00017  .12  3.1  3.9  4.8  242  .16 


The  pilot  describing  function  is: 

Y  f  .  179(0) (2) (3. 25) (6 .39) (12, 4) (12. 6) [-,866, 23.1] 

P*1  ;  "  (0)(1.99)(2)(6.46)(12.4)(12.4)[  .367,23.3]  (42 . 5)  (70) 

where  the  notation  (a)  -  (s  +  a)  and  [f,wn]  -  (s^  +  2fcons  +  w^).  With  a 
ninth-order  denominator  this  transfer  function  is  of  very  high  degree 
indeed.  Explanations  for  each  of  these  factors  is  now  provided. 

Two  orders  stem  from  the  ratio  of  second-order  Pade  factors  repre¬ 
senting  the  pure  time  delay.  This  is  reflected  in  exact  form  in  the 

numerator  [-.866,  23.1],  which  represents  0.075  sec  of  the  given  r  «  0.15 

sec  pure  time  delay.  The  denominator  quadratic  [.367,23.3]  is  the  other 
half  of  this  Pade  approximate  as  it  is  shifted  from  its  open  loop  value  of 

[.866,  23.1]  by  the  loop  closure.  At  low  frequencies  [.367,23.3]  can  be 

replaced  by  a  pure  time  delay  of  2 [.367, 23. 3]  -  0.032  sec.  Thus  0.107  sec 
of  the  initially  prescribed  time  delay  of  0.15  sec  is  "recovered"  by  the 
shifted  version  of  the  Pade  approximation.  If  the  high  frequency  lag  at 
(42.5)  is  treated  as  a  pure  delay  at  low  frequencies  it  will  add  0.023 
sec,  giving  a  total  of  0.130  sec  for  the  effective  human  operator  pure 
time  delay. 

The  terms  associated  with  the  driving  noise  Yv  and  with  the  neuro¬ 
muscular  lag  at  Tfj  sec  will  be  exactly  cancelled;  in  this  example  these 
terms  are  (2) (12 ,4)/(2) (12 .4)  .  These  exact  cancellations  can  be  explained 
by  the  way  in  which  the  (si  —  A^)"^B^  block  appears  in  the  linear-predic¬ 
tor  portion  of  the  OCM  in  Fig.  17.  The  driving  noise  terms  cancel  because 
the  input  to  this  block  does  not  excite  these  terms.  The  neuromuscular 
mode  cancels  because  this  mode  is  redundant  with  the  sTfj  +  1)*1  block 
appearing  above  the  vu  input. 

An  exact  cancellation  also  occurs  in  this  example  at  the  controlled 
element  dynamics:  (0)/(0).  The  Yc  dynamics  occur  in  both  the  KBF  and 
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linear-pred-ctor  portions  of  the  OCM.  This  redundancy  gives  rise  to  can¬ 
cellations.  Despite  the  exact  cancellation  in  this  and  other  examples, 
we  conjecture  that  in  general  this  cancellation  is  only  approximate. 


Of  the  remaining  terms  in  the  describing  function,  (6 . 39 ) ( 12 . 6)/ 
[ (6 . 46) ( 12 . 4) ]  are  close  approximations  and  can  be  dropped  in  any  rational 
simplifying  procedure.  The  (12.4)  derives  from  the  desire  to  set  up  a 
first-order  neuromuscular  lag  of  T^  =  .08  sec  (see  Fig.  17).  In  the  OCM 
this  is  supposed  to  be  accomplished  by  the  selection  of  the  weight  G  on 
the  control  rate  in  the  performance  index  of  Eq .  22.  Interestingly,  the 
neuro -muscular  mode  does  not  survive  the  process  of  reducing  all  the  OCM 
elements  to  a  transfer  function.  Clearly  the  selection  of  G  has  a  dif¬ 
ferent  effect  than  was  intended  by  the  OCM's  originators  in  this  respect. 


When  all  of  these  cancellations  and  near  cancellations  are  taken 

account  the  describing  function  for  pilot  dynamics  for  frequencies 

than  about  15  rad/sec  is: 

Y  (s)  -  4.17(3^25)  _  130S 

Vsj  (1.99) 


into 

less 


(71) 


This  result  agrees  remarkably  well  in  form  with  the  crossover  model: 


^p^c 


~T„S 


(72) 


The  close  to  cancelling  low-frequency  lag-lead,  (  3 . 25)/(  1 . 99)  in  the  OCM 
can  be  rationalized  as  a  slight  improvement  on  the  simple  crossover  model 
in  that  a  tiny  bit  of  the  so-called  low-frequency  phase  droop  will  also 
be  picked  up.  This  is,  of  course,  handled  in  the  extended  crossover  and 
the  structural  -  isomorphic  models  by  the  "alpha"  or  "low-frequency  trim¬ 
ming"  terms  (Ref.  59). 

If  the  control  rate  weight  doesn't  really  "adjust  the  neuro-muscular 
time  constant",  what  does  it  do?  For  this  case,  at  least,  we  can  specu¬ 
late  that  the  value  of  G  and  the  use  of  control  rate  in  the  performance 
index  implicitly  "specifies"  a  frequency  range  in  the  region  of  crossover 
where  the  OCM  will  behave  similarly  to  the  crossover  model.  Other 
examples  we  have  worked  tend  to  support  this  speculation,  although  still 
others,  sucn  as  cases  which  involve  co  trolled  elements  of  the  form 
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77 - 7 >  don't  work  out  as  well.  Consequently  this  is  currently  an  open 

1 s  >  wn  J 

issue.  But  this  Program  CC-based  version  of  the  OCM,  with  its  features 
to  determine  the  dynamic  form  of  the  pilot's  describing  function  and  to 
reduce  it  to  essentials,  permits  a  mucft  wider  cadre  of  researchers  to 
easily  examine  such  issues . 

E.  CONCLUDING  REMARKS 

The  algorithmic  and  computational  advantages  of  the  optimal  control 
model  make  it  extremely  valuable  as  a  means  to  make  quantitative  estimates 
of  the  human  operator's  dynamic  response  in  control  tasks  for  which  the 
model  is  appropriate.  Besides  the  need  to  simplify,  as  illustrated  above, 
there  are  three  other  aspects  which  give  some  difficulty.  The  first  is 
philosophical  and  relates  to  the  explicit  requirement  that  the  human  oper¬ 
ator  description  contain  a  complete  internal  model  of  the  human's  intrin¬ 
sic  characteristics  and  the  system  dynamics  and  disturbances.  Thus,  for 
the  state  estimation  to  be  accomplished,  the  A,  B,  C,  D,  and  E  matrices 
plus  the  system  disturbances  and  the  human  time  delay,  observation  noise, 
and  motor  noise  must  all  be  known.  Further,  for  the  controller  equaliza¬ 
tion  adjustments,  the  A  and  B  matrices  plus  the  weights  in  the  cost  func¬ 
tional  are  needed.  All  of  this  amounts  to  an  essentially  complete  "know¬ 
ledge"  by  the  human  of  the  man-machine  system  characteristics.  Internal 
models  have  a  long  history  in  psychology  for  several  purposes.  For 
instance,  their  elaboration  and  refinement  have  served  as  a  useful  con¬ 
struct  for  the  development  of  skill  by  dint  of  training.  In  fact,  even 
the  simple  crossover  model  can  be  interpreted  as  an  implicit  internal 
model  of  the  human  and  controlled  element  dynamic  characteristics  in  the 
crossover  region.  The  key  problem  is  thus  not  with  the  concept  of  an 
internal  model,  but  rather  its  degree  of  perfection,  especially  in 
extremely  complex  systems  where  the  required  internal  model  is  equally 
complicated.  This  philosophical  point  is  more  sophistic  than  practical  as 
long  as  the  OCM  is  considered  to  be  only  a  mimic  rather  than  an  indication 
that  all  good  human  operators  incorporate  a  Kalman  Bucy  filter  as  a  com¬ 
ponent  of  their  neurological  apparatus! 


The  second  difficulty  is  that  of  attempting  to  identify  the  under¬ 
lying  model  parameters  from  experimental  data.  Not  only  is  this  inverse 
problem  fundamentally  difficult,  but  the  optimal  control  model  reviewed 
here  suffers  from  overparameterization.  Thus,  from  an  identification 
viewpoint,  the  observation  and  motor  noises  are  not  resolvable,  and  the 
feedback  matrix  and  the  observer  gain  matrix  can  only  be  determined  up  to 
a  similarity  transformation  of  the  model  (Ref.  60). 

The  third  problem  area  is  specification  of  the  cost  function.  The 
teleological  character  of  the  linear  quadratic  optimal  model  is  imperfect 
because  the  performance  criterion  must  be  shaped  to  the  task.  As  a  prac¬ 
tical  matter,  this  has  seldom  posed  a  serious  problem  when  the  model  has 
been  applied  by  an  experienced  practitioner.  Nonetheless,  an  aura  of 
artistry  is  present  in  this  requirement. 

In  the  structural  isomorphic  model,  a  very  large  number  of  experimen¬ 
tally  observed  phenomena  are  accounted  for.  Since  its  inception,  a  great 
deal  of  effort  has  been  devoted  to  similarly  account  for  human  operator 
behavior  with  the  algorithmic  model.  This  has  required,  in  the  main, 
adjustments  in  the  cost  function  or  in  those  properties  associated  with 
the  human  operator's  limitations,  such  as  normalized  observation  or  motor 
noise.  The  model  has  proved  to  be  quite  flexible  in  accommodating  most  of 
the  many  behavior  changes  desired.  Table  5  summarizes  some  procedures  and 
techniques  which  have  been  found  suitable  to  accomplish  this  accommodation 
(Refs.  36-38,  41-45).  Thus  advanced  modeling  features,  can  be  handled 
with  the  OCM.  Consequently  both  the  structural  and  algorithmic  forms  of 
pilot  model  are  now  quite  mature  and  can  be  used  in  a  complementary 
fashion  to  solve  pilot -vehicle  analysis  problems  and  to  help  resolve  data 
interpretation  issues. 
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TABLE  5.  PROCEDURES  FOR  ADJUSTMENT  OF  THE  ALGORITHMIC  MODEL 


FEATURES  TO  BE  MODELED 

SUITABLE  PROCEDURES  AND  MEANS 

Effective  time  delay 
a  ccoraoda  t ion 

Least  squares  prediction  applied  to  output  of  Kalman  estimate 
of  delayed  states 

Basic  crossover  behavior 

Use  of  control  rate  weighting  in  distinction  to  control 
weighting  in  cost  function 

Effective  neuromuscular 

iag  Tn 

Select  ratio  of  control  weighting  to  control  rate  weighting 
(e.g.,  "g")  in  cost  function 

Selection  of  cost 
function  weights  on 
scates  and  control 

Choose  weights  to  be  inverse  of  squares  of  the  respective 
maximum  allowable  values 

Remnant 

Observation  noise  covariances  scaled  with  mean-squared  state. 
Residual  (non-scaled)  observation  noise  component  to  account 
for  imprecision  due  to  lack  of  references. 

Motor  noise  to  reflect  inability  to  generate  control  motions 
precisely. 

Residual  motor  noise  to  reflect  human's  introduction  of  noise 
into  an  undisturbed  system. 

Low-f req uen cy  phase  lag 

Use  larger  motor  noise  level  than  actually  present  in  deter¬ 
mining  Kalman  filter  gains 

Perceptual  and 
indifference  thresholds 

Scale  observation  noise  inversely  with  equivalent  gain 
(random  input  describing  function  for  threshold) 

Scale  observation  noise  inversely  with  attentlonal  frac¬ 
tion  <  f  ^ )  of  each  display,  subject  to  the  constraint  that 

<r£i>  +  fmargin  <  fl  >  °- 

Different  noise  levels  for  foveal  and  parafoveal  viewing. 

Workload  (  a 1 1 e n t Iona  1 ) 

Attentlonal  workload  effects  evaluated  by  examining  perform¬ 
ance  as  a  function  of  the  reserved  workload  margin,  f 

°  ma r  g 1 n 

Mo  t i on  cues 

Add  model  of  human  motion  sensory  apparatus  (e.g.,  vestibular 
system,  proprioception)  to  state  and  output  equations. 
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SECTION  V 


ESTIMATION  OF  PILOT  RATINGS 

In  the  previous  sections  we  have  described  techniques  for  the  estima- 
tion  of  pilot  dynamics  in  closed- loop  tasks.  Because  aggressively 
performed  closed- loop  tasks  are  ordinarily  critical  from  the  standpoint  of 
pilot  compensation  or  skill  required,  and  are  usually  high  workload  flight 
phases,  they  tend  to  be  dominant  discriminators  in  flying  qualities 
assessments.  The  assessments  themselves  are  provided  by  pilot  comments 
and  associated  ratings,  such  as  the  Cooper-Harper  scale  (Ref.  61)  shown  in 
Fig.  19.  It  is  apparent  that  pilot  compensation  (equalization)  and 
effort  (workload)  and  task  performance  are  major  constituents  of  the 
rating  scheme.  When  the  task  variables  (effective  vehicle  dynamics, 
forcing  functions  and  disturbances,  etc.)  are  well-defined  the  pilot- 
vehicle  system  dynamic  models  presented  in  Sections  II -IV  can  be  used  to 
make  quantitative  estimates  of  pilot  compensation  and  task  performance. 
Workload,  on  the  other  hand,  is  much  more  difficult  to  quantify.  Still, 
we  can  expect  some  connections  between  the  subjective  pilot  ratings  and 
comments,  which  are  subjective  workload  indices,  and  the  pilot  and  pilot- 
vehicle  system  dynamics  and  performance.  These  connections  are  intrin¬ 
sically  empirical.  They  are  also  awkward  theoretically  because  the  rating 
scale  is  ordinal.  Consequently  averages,  standard  deviations,  etc.  are 
not  legitimate  statistics,  although  this  has  never  stopped  flying  qual¬ 
ities  engineers  from  using  them!  (Fortunately  the  scales  seem  to  be  close 
to  interval  in  some  ranges  or,  for  the  purist,  data  can  be  converted  to  an 
underlying  interval  scale  wherein  all  the  parametric  statistics  can  be 
applied  and  then  converted  back  --  see  Ref.  62.) 

Our  goal  in  this  section  is  to  summarize  the  available  connections 
between  pilot  and  pilot-vehicle  system  dynamics  and  pilot  ratings.  There 
are  fundamentally  two  approaches  which  have  been  used  with  some  success. 
The  first  directly  associates  pilot  and  system  dynamic  and  performance 
characteristics  with  the  pilot  rating  via  a  functional  relationship.  Such 
functionals  have  been  developed  for  use  with  both  the  classical  and  OCM 
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Aircraft 

Characteristics 


Demands  on  the  Pilot  Pilot 

in  Selected  Task  or  Required  Operation*  Rating 


Excellent 

Highly  Desirable 

Pilot  Compensation  Not  a  Factor  for 

Desired  Performance 

1 

Good 

Negligible  Deficiencies 

Pilot  Compensation  Not  a  Factor  for 

Desired  Performance 

2 

Fair  -  Some  Mildly 
Unpleasant  Deficiencies 

Minimal  Pilot  Compensation  Required  for 
Desired  Performance 

3 

Minor  But  Annoying 
Deficiencies 

Desired  Performance  Requires  Moderate 

Pilot  Compensation 

4 

Moderately  Objectionable 
Deficiencies 

Adequate  Performance  Requires 
Considerable  Pilot  Compensation 

5 

Very  Objectionable  But 
Tolerable  Deficiencies 

Adequate  Performance  Requires  Extensive 
Pilot  Compensation 

6 

Major  Deficiencies 

Adequate  Performance  Not  Attainable 

With 

Maximum  Tolerable  Pilot  Compensation. 
Controllability  Not  in  Question 

7 

Major  Deficiencies 

Considerable  Pilot  Compensation  Is 

Required  for  Control 

8 

Major  Deficiencies 

Intense  Pilot  Compensation  Is  Required  to 
Retain  Control 

9 

Major  Deficiencies 

Control  Will  Be  Lost  During  Some  Portion 
of  Required  Operation 

10 

‘Definition  of  required  operation  involves  designation  of  flight  phase  and/or 
subphase  with  accompanying  conditions. 


Figure  19.  Cooper  -  Harper  Handling  Qualities  Rating  Scale 


versions  of  pilot  models.  The  actual  connections  which  have  been  estab¬ 
lished  are  based  on  specific  tasks  and  circumstances. 


The  second  approach  is  more  clinical  in  style.  It  takes  into  account 
the  pilot  and  pilot-vehicle  system  characteristics  in  terms  of  their 
implications  for  control.  A  list  of  assessment  features  is  considered  to 
reveal  symptoms  of  flying  qualities  problems.  Some  quantitative  aspects 
can  be  set  forth,  but  others  are  only  qualitative.  Consequently  this 
,  appfpa^  f  i,s  more  the  basis  for  a  pseudo  pilot  commentary  rather  than  a 
means  to  make  numerical  rating  estimates  directly.  Of  course,  if  the 
"commentary"  is  sufficiently  complete  it  can  be  converted  to  a  rating  by 
working  through  Fig.  19.  The  clinical  technique  is  especially  useful  to 
define  possible  flying  qualities  problems  and  key  effective  airplane 
dynamic  parameters,  or  as  a  means  of  interpreting  experimental  data. 

The  two  approaches  described  are  currently  most  useful  for  single 
axis  situations.  Multi-axis  rating  estimates  can  be  developed  from  single 
axis  results  using  a  "product  rule,”  which  is  described  in  the  third 
article.  OCM-based  multi-axis  results  can  also  be  the  basis  for  direct 
estimates  of  multi-axis  ratings.  This  is  summarized  in  the  third  article 
as  well,  and  is  developed  in  Appendix  A. 

A.  PILOT  RATING  FUNCTIONALS 


A  direct  approach  is  to  formulate  a  functional  which  incorporates  the 
pilot  and  system  dynamic  and  performance  quantities  which  are  presumed  to 
underlie  the  pilot  rating.  A  general  form  which  explicitly  contains  some, 
and  implicitly  contains  all  of  the  desired  features  is  given  by, 


Mission/Task 

Performance 


Pilot  Workload 


Dominant  Aircraft  Pilot  Activity  Pilot  Equalization 

Motion  Quantities  (Scale  of  (Dynamic  Quality 

and  Task  Measures  Pilot  Effort)  of  Pilot  Effort) 
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r2 

•2 
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w) 

;  •••  J 

"lk 

(73) 


80 


The  subscript  notation  used  is  intended  to  imply  that  if  I  motion  and  task 
measures  are  controlled  by  keK  pilot  loops  actuating  jfj  control  points. 
This  functional  form  is  general  enough  to  include  the  existing  (e.g., 
Refs.  62-72)  approaches  to  quantitative  flying  qualities  rating  functions. 
The  key  closed- loop  system  quantities  in  the  rating  functional  are 
measures  of  mission/task  performance.  These  are  conveniently  described  by 
a  set  of  dominant  weighted  aircraft  motion  deviations  and  total  task 
accuracy  or  error  indications  (represented  by  the  q^)  . 

~2 

The  pilot  activity  component  of  pilot  effort,  6^  (either  force  or 

^  #  2 

displacement ,  as  pertinent  to  the  manipulator  involved)  and  fij  are  partic¬ 
ularly  dependent  on  the  level  of  pilot  gain.  For  a  given  gain,  these  will 
increase  directly  with  gust  disturbance  spectrum  amplitude  and  remnant 
amplitude.  Accordingly,  both  the  mission/task  and  pilot  activity  quan¬ 
tities  will  reflect  turbulence  and  remnant  levels. 

The  pilot  equalization  component  of  pilot  workload  is  represented  in 
Eq.  18  by  the  slope  (in  dB  per  octave  or  decade)  of  the  pilot's  amplitude 
ratio  evaluated  at  a  particular  frequency  (generally  near  crossover). 
This  is  by  no  means  the  only  measure  available  to  describe  the  dynamic 
quality  of  the  pilot's  effort,  others  (e.g.,  Refs.  64-68)  use  pilot  lead 
time  constants,  a  desirable  alternative  for  particular  situations  with  a 
sufficient  data  base.  Then  the  rating  functional  takes  the  very  useful 
form  illustrated  in  Fig.  20.  At  present,  adequate  functions  of  this  form 
exist  for  precision  hover  tasks  (Refs.  62-68),  pitch  attitude  control 
(Ref.  66),  and  roll  attitude  control  based  on  Refs.  69  and  70.  In  addi¬ 
tion,  the  Ref.  71  data  provide  a  base  for  a  multiloop  functional. 

The  technique  pioneered  in  Ref.  64  actually  used  the  pilot  rating 
functional  as  a  performance  index,  as  well  as  a  rating  estimator.  That 
is,  the  pilot  model  parameters  (Fig.  20a)  were  adjusted  to  minimize  R,  the 
pilot  rating  functional. 

The  follow-on  work  of  Ref.  68,  which  was  dedicated  to  experimentally 
verifying  the  Ref.  64  result,  produced  a  "modified"  pilot  rating  func¬ 
tional  for  the  R2  +  R3  component,  as  shown  in  Fig.  21.  The  correlation 
between  predicted  and  actual  ratings  shown  in  Fig.  22  are  reasonably  good. 
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a  )  Pilot  /  Vehicle  Model 


R  1  I  +  R  i  +  R2  R3 
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b)  Pilot  Rating  Functional 
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J  =  I  +  PERF  ♦  R2  ♦  R3 
P ERF,  Rz  ,  R3  -  os  defined  in  b) 

1)  KPi  ,  TLj  ,  KPg  ,  7j lq  found  lo  minimize  J 

2)  KPj  ,  Kpe  simultaneously  increosed  by  20%. 

If  system  asymptotically  stable ,  performance  and  rating  computed  based  on 
values  found  in  I).  If  system  not  stable ,  go  to  3). 

3)  TLe ,  TLt  held  fined  ot  values  found  in  I).  Kpi  ,  Kpg  simultaneously  decreased 
until  values  found  for  which  a  20.3%  increase  would  result  in  an  unstable 
system  while  a  20%  increase  would  result  in  a  stable  system.  These  values 
are  used  to  compute  performance  and  rating. 

cj  Minimum  Pilot  Rating  Concept  For  Pilot  Mode!  Parameter  Adjustment 
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Actual  Pilot  Rating 


1.8 

(TLe  +  Tl'e)  (sec) 


Modification  for  Hover  Paper  Pilot  Rating 
Functionals  (Ref.  68) 


Figure  21.  "Paper  Pilot"  Rating  Functionals 


I  23456789 
Pilot  Rating  Computed  Using  Measured  System  Properties 
and  the  Modified  Paper  Pilot  Rating  Functional  With”Equiva!enl',TL')l 


Figure  22.  Comparison  of  Actual  Pilot  Ratings  with  Ratings  Computed 
Using  Modified  Hover  "Paper  Pilot"  Rating  Functional  (Ref.  68) 
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In  Ref.  66  the  task  was  changed  to  pitch  attitude  control  and  the 
resulting  pilot  rating  functional  evolved  to: 


R  - 


R1 

0.1 


. 974— a 


R2 

2 . 5Tl  +  1.0 


where 


a  -  o&/o\  -  ratio  of  error  variance  to  input  variance 

Tl  -  pilot  lead,  seconds 


and 


R  <  10  ,  0  <  R2  <  3.25  ,  0  <  Ri 


(If  a  value  of  a  <  .974  cannot  be  obtained,  R  -  10) 


Yet  another  way  to  estimate  pilot  ratings  is  to  use  correlations 
developed  for  the  algorithmic  pilot  model.  (Refs.  72  through  75.)  This 
pilot  rating  estimation  procedure  is  based  on  the  hypothesis  that  the 
pilot  rating  for  a  particular  task  and  set  of  vehicle  dynamics  can  be 
correlated  with  the  numerical  value  of  the  index  of  performance  (minimum 
values  of  the  OCM  Cost  Function)  resulting  from  the  optimal  pilot  modeling 
procedure.  As  indicated  in  Fig.  23,  this  has  worked  fairly  well  for  some 
single-axis  cases  (e.g.,  helicopter  hover  and  longitudinal  approach). 


The  extension  of  the  OCM  performance- index-based  pilot  rating  esti¬ 
mating  procedure  to  the  multi-axis  case  has  been  addressed  as  part  of  this 
project.  The  basic  developments  are  given  in  Appendix  A  and  Volume  III. 
Reference  76  is  the  primary  source  of  connected  single  and  multi-axis 
rating  data,  although  no  pilot  dynamic  information  is  available.  Conse¬ 
quently  the  OCM  was  used  to  establish  pilot  and  system  dynamics  estimates. 
For  the  rating  estimates  only  the  performance  index  is  needed.  The  appro¬ 
priate  performance  index  for  each  single  axis  was  chosen  to  be, 


'  axis ; 


(74) 


And  for  a  multi-axis  task,  the  objective  function  used  was 


J  task 


'axes 


\ 


^axis  i 


(75) 
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VALUE  OF  INDEX  OF  PERFORMANCE 


Model  Parameters  for  Hover  Task  of  Miller  and  Vinje  (1968) 


Parameter 


Value 


Time  delay  r 

Neuromuscular  time  constant  TN 
Visual  thresholds 

"Full-attention''  noise-signal  ratio  lor  observation  noise 
Fraction  ol  attention  on  control  task  fe 
Noise-signal  ratio  lor  motor  noise 


0.2  s 
0.2  s 
None 
0.0025 

0.25  —  1.0  (configuration  dependent) 
0.003 


Index  ol  Performance 
j  =  Ef-  !  l'  [yT(t)Qy(t)  +  u?(t)Ru(t)]dt 

It—*  t  'o 


y,  -  S 
yi  =  * 
h  »  * 

u,  »  commanded  control 


q„  -  (1/0.0873)’  s*/rad* 
=  (1/3.5t»/ft* 
q„  -  (1/3.51*  s* /It* 
r„  =  (1/0. 138F/IF 


Figure  23.  Pilot  Rating  vs.  Value  of  Model  Index  of  Performance  (Ref.  72) 
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The  justification  of  this  selection  involves  three  considerations.  The 
first  relates  to  the  selection  of  equal  (unity)  weighting  on  each  Jaxis  i 
in  the  definition  of  Jtask  *-n  multi-axis  tasks.  This  decision  was  based 
on  the  instructions  given  to  the  subjects  in  the  Ref.  64  experiment.  They 
were  to  attempt  to  minimize  the  errors  in  all  controlled  axes.  That  is, 
they  were  instructed  that  no  axis  was  to  be  given  preference,  which  would 
then  define  primary  and  secondary  sub- tasks. 

Secondly,  the  normalization  of  the  mean- square  error  with  the  mean- 
square  command  deals  nicely  with  the  fact  that  different  units  and 
different  command- signal  strengths  were  used  in  the  axis. 

Finally,  the  interpretation  of  g^  requires  some  discussion.  In  the 
OCM,  the  selection  of  g^  defines  the  frequency  range  over  which  the  open- 
loop  system  amplitude  ratio  approximates  a  K/s-like  form.  In  connection 
with  the  OCM  it  is  often  cited  that  g^  is  selected  to  yield  a  desired 
neuromotor  time  constant,  Tn,  in  the  pilot's  describing  function  obtained 
from  the  model.  But,  as  indicated  in  Section  IV,  when  the  total  pilot 
describing  function,  Yp ,  is  actually  constructed  from  its  various  elements 
in  the  OCM,  the  Tn  established  in  this  fashion  is  canceled  by  a  directly 
compensating  lead,  leaving  the  actual  estimated  Yp  with  no  (Tn  s  +  1)_1 
lag.  Still  it  has  been  convenient  to  adjust  g^  in  this  fashion  even 
though  the  lag  will  later  disappear.  In  this  vein,  the  value  of  the 
desired  neuromotor  time  constant  used  is  either  0.1  sec,  or  the  Tn  that 
yields  the  lowest  error  (e.g.,  best  performance),  whichever  is  greater . 
Notice  that  after  Tn  is  determined  in  the  above  fashion,  this  "operating 
point"  is  associated  with  some  weight  g^  in  Jax£s  i-  This  value  may  also 
infer  the  subject's  subjective  trade  between  performance  (a£/  and  workload 
(o$)  ■  And  since  pilot  lead  and  are  correlated,  this  procedure  maxi¬ 
mizes  the  possibility  of  relating  the  resulting  value  of  Jtask  t0  the 
subjective  rating  of  the  task. 

Shown  in  Fig.  24  is  the  correlation  between  Jfask  as  modeled,  and 
the  subjective  ratings  of  the  task.  The  correlation  between  Jtask  and  P0R 
from  the  single-axis  results  appears  to  hold  for  the  multi-axis  results  as 
well.  This  result  seems  to  indicate  that  the  ratings  reflect  the  actual 
oerformance  and  workload  (stick  rate)  in  the  overall  task.  The  results 


-2.0  -1.0  0.0 
Log  J tosk 


Figure  24.  Pilot  Rating  vs.  Performance  Index  for  Dander 
Single  and  Multi-Axis  Tasks  (See  Appendix) 

also  tend  to  support  the  hypothesis  that  determining  the  weightings  g^  in 
the  manner  discussed  leads  to  the  "correct"  relative  weightings  on  control 
rate  in  the  axis,  and  the  relative  weight  between  control  rate  and  normal¬ 
ized  error.  Because  the  multi-axis  correlations  follow  the  same  trend  as 
the  single-axis  data  this  study  indicates  that  the  objective  function  for 
multi-axis  situations  can  be  extrapolated  (or  calibrated)  from  single  axis 
correlations . 

B.  WORKLOAD,  ATTENTIONAL  DEMANDS,  AND  THE 
PRODUCT  RULE  FOR  MULTI -AXIS  RATINGS 

There  is  a  strong  connotation  of  increasing  pilot  effort  and  workload 
in  the  phrases  of  the  Cooper-Harper  Scale  (Ref.  61)  which  invoke  levels  of 
"pilot  compensation,"  but  workload  is  difficult  to  define  and,  con- 
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sequently,  Co  quantify.  A  general  definition  that  can  be  measured  and 
predicted  is  workload  margin,  defined  as  the  ability  (or  capacity)  to 
accomplish  additional  (expected  or  unexpected)  tasks.  The  pilot  opinion 
rating  scale  satisfies  this  definition  up  to  its  "uncontrollable"  limit 
point.  It  is,  therefore,  a  key  workload  measure,  easy  to  obtain  in  some 
experimental  circumstances. 

Auxiliary  tasks  have  been  developed  that  satisfy  the  workload  margin 
definition  given  above  and  that  permit  more  objective  measurements.  One 
such  task  provides  a  complementary  pair  of  measures  suitable  for  integrat¬ 
ing  many  workload  concepts  and  factors  into  one  basic  context.  These  are 
the  "attentional  demand"  and  the  "excess  control  capacity." 

The  attentional  demand  and  excess  control  capacity  measures  have  been 
connected  with  pilot  rating  in  a  multiaxis  experiment  using  the  so-called 
cross -coupled  subcritical  task  (see,  for  example.  Refs.  15  and  62).  A 
block  diagram  of  the  general  experimental  setup  is  shown  in  Fig.  25.  The 
pilot  first  performs  the  primary  task  alone,  attempts  to  achieve  satisfac¬ 
tory  levels  of  performance,  and  provides  a  Cooper-Harper  pilot  rating. 


eC 

Figure  25.  Single-Loop  Primary  Task  with  Secondary 
Cross-Coupled  Loading  Task 
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The  secondary  subcritical  tracking  task  is  then  connected  in  order  to 
"load"  the  pilot.  The  difficulty  of  the  secondary  task  is  made  propor¬ 
tional  to  primary  task  performance  via  the  cross-coupling.  Thus,  when  the 
pilot  keeps  primary  task  performance  less  than  a  criterion  value  (based 
on  the  runs  with  the  primary  task  alone),  the  secondary  task  difficulty  is 
automatically  increased  by  increasing  the  rate  of  divergence  of  the  sec¬ 
ondary  task  instability.  Conversely,  when  the  pilot  becomes  so  busy  with 
the  secondary  task  that  the  primary  task  error  becomes  larger  than  the 
criterion  value,  the  secondary  task  difficulty  is  automatically  decreased. 
The  final  "score"  is  ys ,  the  stationary  value  of  the  secondary  unstable 
pole  (7)  in  rad/sec.  The  scores  obtained  from  this  cross -coupled  second¬ 
ary  task  represent  its  difficulty;  consequently,  they  also  represent  the 
"degree  of  ease"  of  the  primary  task  or  the  excess  control  capacity  avail¬ 
able  with  respect  to  the  primary  task.  The  ys  scores  can  be  appropriately 
scaled  into  proportional  workload  indices  by  normalizing  them  with  respect 
to  the  maximum  sidetask  score  attainable  under  full  attention  conditions 
(no  primary  task).  In  this  case,  fs  approaches  7C ,  the  "critical  task" 
score.  The  attentional  demand  of  the  primary  task  is  then  given  by 


The  attentional  demand,  AD,  is  a  dimensionless  decimal  fraction  that 
can  be  equated  with  the  average  primary  control  task  attentional  dwell 
fraction,  rj .  Its  complement,  the  excess  control  capacity,  which  measures 
the  average  fraction  of  time  available  for  other  than  the  primary  task, 
is 

As 

XSCC  =  —  (76) 
As 

If  the  side  task  is  taken  to  be  a  surrogate  for  all  of  the  managerial 
functions,  XSCC  will  be  just  the  average  managerial  task  dwell  time  frac¬ 
tion  1-r;. 

Achieving  the  critical  limiting  score  in  the  cross-coupled  secondary 
task  indicates  a  condition  of  maximum  available  excess  control  capacity; 
the  secondary  task  is  a  "critical"  task  in  this  limiting  case.  The  criti¬ 
cal  task  provides  a  divergent  controlled  element  form  that  tightly  con¬ 
strains  allowable  pilot  equalization  near  the  region  of  gain  crossover  so 
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that  the  pilot's  effective  time  delay,  re,  is  the  sole  determinant  of 
system  stability.  Thus,  pilot  activity  that  demands  an  increase  in  re  on 
the  whole  task  will  prevent  the  attainment  of  the  pilot's  critical  limit¬ 
ing  score  on  the  cross -coupled  secondary  task. 

Secondary  scores  obtained  for  a  variety  of  primary  controlled  ele¬ 
ments  are  presented  in  Ref.  62.  Figure  26  shows  how  the  scores  for  the 
best  gain  configurations  of  each  controlled  element  compare  with  the 
Cooper-Harper  ratings.  In  Fig.  26  a  score  of  As  -  0  corresponds  to  100 
percent  of  the  pilot's  attention  being  devoted  to  the  primary  task  or  no 
excess  control  capacity;  whereas,  a  limiting  score  (As  -  5.5)  means  that 
no  attention  is  required  to  maintain  primary  task  performance  or  that  100 
percent  of  excess  control  capacity  is  available. 

These  relationships  show  that  subjective  pilot  ratings,  can  be  asso¬ 
ciated  closely  with  the  objective  measures  of  workload  provided  by  the 
attentional  demand  and  the  excess  control  capacity.  The  lower  (better) 
values  of  pilot  rating  correspond  to  low  attentional  demands  and  large 
excess  capacity  to  perform  other  functions.  More  difficult  effective 
vehicle  dynamics  that  receive  poorer  pilot  ratings  of  their  flying  quali¬ 
ties,  require  much  more  of  the  pilot's  attention  and  hence  leave  less 
capacity  for  other  tasks. 


The  excess  control  capacity  concept  also  provides  a  potential  basis 
for  estimating  ratings  for  multiloop  situations  Ref.  77.  First,  assume 
that  the  relationship  between  pilot  rating  and  excess  control  capacity, 
An  »  A  S/Ac  given  by  Fig.  26,  is  applicable  to  each  axis  in  a  multi-axis 
situation.  Then,  single-axis  capacity,  or  attention,  values  can  be  com¬ 
bined  to  yield  the  combined  axis  value  by  a  multiplication  process,  i.e., 
the  multiaxis  excess  capacity,  Ap^ ,  is  given  by  the  product  of  the  excess 
capacities  for  the  individual  axes: 


m 

nAn^ 


(77) 
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For  R  -  A  +  as  a  linear  tit  ot  the  Fig.  26  data,  the  multiaxis  rating 
will  be  given  by, 

m  m  _  a! 

Rm  -  A  +  BAnm  -  A  +  Biun.  -  A  +  Bn  { — — J 

1  m  (78) 

Rm  -  A  +  —  n  (Rj.  -  A) 

Combined  ratings  are  always  greater  than  (or  equal  to)  individual  ratings, 

since  combined  An's  are  always  less  than  any  individual  An.  Also,  the 

m 

maximum  value  of  Rm  never  exceeds  A,  i.e.,  for  large  R^  <  A,  II 
(R  -  A)  -  0. 

The  logical  value  for  A  is  10.0,  and  B  is  determined,  using  the 
empirical  data,  to  be  equal  to  -8.3.  As  depicted  in  Fig.  27,  this  results 
in  a  good,  overall  fit  to  all  of  the  available  multiaxis  rating  data 
(Ref.  77).  Notice  that  in  its  final  form  the  multi-axis  rating,  Rm,  can 
be  computed  directly  from  the  single -axis  ratings  R^ .  Measures  or  com¬ 
putation  of  excess  control  capacity  or  attentional  demand  are  not 
required. 

C.  THE  CLINICAL  APPROACH  TO  RATING  ESTIMATION 


The  treatment  above  has  the  great  merit  that,  when  appropriate 
measures  and  experimental  correlates  are  available,  a  set  of  pilot  rating 
estimates  can  be  made  using  relatively  simple  formulas.  The  detailed 
reasons  for  the  rating  estimates  are  inherently  buried  in  the  empirical 
data  which  serve  as  bases  for  the  correlations.  In  other  words,  the  pilot 
commentary  and  reasons  behind  whatever  the  rating  estimate  comes  out  may¬ 
be  quite  obscure.  To  alleviate  this  difficulty,  and  to  provide  an  alter¬ 
native  for  situations  where  the  data  i.-,e  is  insufficient  or  non-existent, 
a  clinical  approach  is  indicated.  Ht  che  characteristics  exhibited  by 
the  pilot  and  pilot-vehicle  system  dynamics  are  examined  for  "symptoms”  of 
potential  problems.  These  are  then  reflected  into  a  summary  of  properties 
whic  amount  to  a  pilot  commentary  expressed  in  technical  terms. 

Consider,  for  the  most  elementary  situation,  that  the  crossover  model 
is  used  in  accomplishing  a  pilot-vehicle  analysis  for  a  given  set  of 


26.  Calibration  of  Pilot  Rating  with  Attentiona!  Workload  and  Excess  Cont 


Predicted  Multi-Axis  Rating 


Figure  27.  Correlations  Obtained  with  Product  Method 
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effective  aircraft  dynamics.  The  data  directly  available  from  the 
analysis  includes  an  estimate  of: 

•  the  stability-limited  maximum  crossover  frequency,  wu; 

•  pilot  lead  equalization  required  in  the  region  of  crossover 
to  make  good  the  crossover  law  (measured  in  terms  of  pilot 
amplitude  ratio  slope,  [ d j Yp | d$/d  'log  w]wc 

•  the  nominal  full-attention  crossover  frequency,  wc; 

The  two  crossover  frequencies  are  closely  related,  i.e. 

fr’cA’u 

No  pilot  lead  0,78 

Low-frequency  pilot  lead  0.66 

If,  in  addition,  equivalent  forcing  function  information  is  available  the 
system  steady-state  performance  can  be  determined  easily  (e.g.  by  using 
Figs .  9  or  13a) . 

As  might  be  expected,  the  most  important  pilot  dynamics  correlates 
with  pilot  rating  are  pilot  gain  and  pilot  lead.  Empirical  connections 
between  these  are  given  in  Fig.  28.  For  a  particular  controlled  element 
there  is  an  optimum  controlled  element  gain  which  depends  on  the  manipu¬ 
lator  dynamics,  controller  sensitivity,  control  harmony  among  axes,  etc. 
No  theory  yet  exists  to  establish  this  optimum  gain,  so  it  must  be  deter¬ 
mined  empirically.  Then,  curves  such  as  those  shown  in  Fig.  28  can  be 
used  to  assess  any  rating  decrements  from  the  optimum.  By  virtue  of  the 
wc-Kc  independence  property  any  change  in  Kc  will  be  countered  by  a  change 
In  Kp  to  keep  the  crossover  frequency  approximately  constant.  However, 
either  too-sluggish  (Kc  too  small,  Kp  too  large)  or  too-sensitive  condi¬ 
tions  can  give  rise  to  major  decrements.  This  can  be  greater  than  6 
rating  points  even  for  the  Yc  -  Kc/s  controlled  element  dynamics.  As  can 
be  appreciated  from  Fig.  28  the  optimum  is  quite  broad  (changes  of  plus  or 
minus  50%  in  either  direction  are  less  than  1  rating  point  for  even  the 
narrowest  U-shaped  curve) ,  so  once  the  controlled  element  sensitivity  is 
properly  adjusted  minor  controlled  element  gain  changes  are  not  major 
factors  in  pilot  rating. 

The  pilot  lead  equalization  required  to  make  good  the  crossover  model 
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has  a  major  effect  on  the  pilot  rating.  For  example,  Fig.  28  indicates 
that  the  difference  between  a  Yc  -  Kc/s  controlled  element,  which  requires 
no  pilot  lead,  and  Yc  -  Kc/s^,  which  demands  +1  lead  units,  is  a  pilot 
rating  decrement  of  about  3  Cooper -Harper  rating  points.  Considered  as 
idealised  systems  these  correspond,  respectively,  to  "rate  command"  and 
"acceleration  command"  effective  vehicle  characteristics.  Re-examining 
Fig.  26,  the  rate  command  system  data  point  shows  a  pilot  rating, 
PR  -  2  1/2  with  an  attentional  demand  of  0.2  while  the  acceleration 
command  system  data  point  has  PR  -  6  and  a  control  task  dwell  fraction  of 
about  0.65.  The  primary  reasons  for  rating  shifts  for  these  data  are  the 
amount  of  lead  required  and  the  reduction  in  system  performance  (the 
attainable  crossover  frequency  for  the  acceleration  case  is  less  than  that 
for  the  rate  command  situation  because  of  the  increased  re  due  to  the  need 
for  lead  generation).  In  any  event,  even  a  best  gain  acceleration  command 
system  will  be  Level  2  (3  1/2  <  PR  <  6  1/2)  from  a  flying  qualities  stand¬ 
point.  From  the  descriptive  adjectival  phrases  of  Fig.  19  this  level  of 
low  frequency  lead  generation  would  therefore  be  interpreted  as  "consider¬ 
able  pilot  compensation"  required  to  achieve  adequate  performance. 

There  are,  of  course,  factors  other  than  pilot  lead  and  gain  adjust¬ 
ment  that  affect  the  pilot  rating.  In  general,  flying  qualities  ratings 
tend  to  be  given  on  a  global  basis  which  may  include  several  maneuvers  in 
a  task  complex.  Both  open-loop  (unattended)  and  closed-loop  (attended) 
piloting  operations  will  be  considered  in  the  rating.  In  the  present  work 
we  are,  of  course,  concerned  primarily  with  the  closed-loop  piloting 
aspects.  In  fact,  for  stability  and  control  flight  testing  the  important 
connection  is  with  "flying  qualities  while  tracking"  aspects  and  other 
precision  and/or  aggressive  tasks  which  involve  tight  closed-loop  pilot- 
vehicle  control. 

The  unattended  category  can  be  the  major  factor  In  determining  the 
acceptable  values  of  very  low  frequency  divergences  such  as  the  spiral  or 
a  divergent  phugoid.  It  can  also  be  decisive  in  setting  the  nature  of  the 
"hold"  characteristic  built  into  the  stability  augmentation  system.  For 
example,  for  many  up  and  away  tasks  the  attended  longitudinal  pitch 
attitude  system  should,  ideally,  require  no  pilot  lead  equalization.  In 
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terms  of  ideal  effective  airplane  dynamics  this  would  approximate 
Yc  -  Kc/s  in  the  region  of  crossover.  But  for  unattended  operations  a 
rate  command  system  is  not  ideal  in  that  an  attitude  -  stable  platform  is 
de°ired.  Thus  a  rate-command/atritude-hold  system  has  superior  pilot 
ratings  to  rate -command/rate -hold . 

Another  major  facet  in  nearly  unattended  operations  is  "command- 
ability",  the  ability  of  the  airplane  to  respond  in  a  precise,  orderly, 
and  predictable  manner  to  highly  skilled,  precognitive  pilot  command 
inputs.  These  inputs  are  pure  commands,  functions  of  time  alone,  and,  as 
such,  are  basically  open-loop  in  character.  Typical  examples  are  turn 
entries,  step-like  (for  attitude  command  systems)  or  pulse-like  (for  rate 
command  systems)  inputs  to  adjust  attitude,  etc.  These  maneuvers  may  have 
to  be  fine-tuned  at  the  end  via  closed-loop  control,  but  for  an  ideal 
vehicle  and  a  skilled  pilot  this  will  not  ordinarily  be  necessary.  Again, 
to  the  extent  that  this  feature  of  the  airplane's  characteristics  enter 
into  the  rating  game,  closed- loop  dynamics  considerations  are  not  explic¬ 
itly  involved. 

An  important  distinction  about  the  unattended  factors  in  the  current 
context  is  that  they  may  set  a  base  for  the  pilot  rating  which  does  not 
depend  on  closed-loop  factors.  This  base  can  itself  shift  as  the  divided 
attention  requirements  shift.  For  example,  if  managerial  tasks  take  up 
almost  all  the  available  time  the  effective  vehicle  dynamics  in  the 
unattended  state  may  have  to  be  highly  automated  even  including  path, 
altitude,  or  position  control.  In  any  event,  the  closed-loop  effects 
should  be  thought  of  as  increments  from  the  base  level  determined  from  the 
unattended  operation  requirements. 

Table  6  presents  listing  of  primary  factors  to  which  the  pilot  is 
sensitive  and  which,  accordingly,  underlie  pilot  rating.  Except  for  pilot 
lead  and  gain  variation  from  optimum  these  factors  are  not  individually 
quantifiable  in  ratings  terms.  On  th^.  other  hand,  witv>  ^^dern  flight 
control  system  technology  most  of  them  can  be  modified  by  design.  Conse¬ 
quently  these  system  aspects  can  be  profitably  compared  in  competing 
system  studies,  and  also  serve  as  a  useful  checklist  for  interpreting 
mannad  simulation  or  flight  test  results.  As  remarked  earlier,  a  pseudo 
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pilot  commentary  can  be  constructed  by  considering  them. 

In  Table  6  both  items  under  "Unattended  Operations"  and  the  "Pilot 
Lead"  and  "Pilot  Gain/Optimum"  parts  of  the  "Attended  Operations"  list 
have  been  covered  above.  The  remaining  items  will  be  discussed  below. 
Some  of  the  considerations  can  be  developed  from  the  crossover  model, 
while  .•'i»ers  will  require  application  of  the  structural  -  isomorphic  pilot 
moiti.  in  some  form  or  other. 


TA£LE  6.  PILOT-VEHICLE  SYSTEM  FACTORS  IN  PILOT  RATING 


ATTENDED  OPERATIONS 

Pilot  Lead 

Pilot  Gain/Optimum 

Urgency  Adjustment  Gain  Tolerance 
Without  Changing  Closed-Loop 
Dynamic  Form 

Stability  Margin  Gain  Tolerances 
Including  Total  Available  Gain  Range 

Neuromuscular  System  Coupling 

Attentional  Demands/Excess 
Control  Capacity 

Closed-Loop  System  Performance 


•  UNATTENDED  OPERATIONS 

Allowable  Fluctuations  in 
Pi  lot -Control  -  Precision  Demands 

Equilibrium/Trim  Properties  (Effective 
"Hold"  Characteristics) 
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The  "Urgency  Adjustment  Gain  Tolerance"  factor  can  best  be  understood 
by  considering  two  limiting  cases  of  controlled  element.  For  the  first, 
consider  a  Kc/s  controlled  element  form.  From  the  crossover  model  The 
pilot  dynamic  characteristic  for  this  system  will  be  a  pure  gain  plus 
effective  time  delay.  The  closed-loop  system  for  this  case  can  support  a 
range  of  pilot  gains  which  correspond  to  crossover  frequencies  from  zero 
to  an  octave  or  so  below  wu  with  only  minor  changes  in  the  basic  dynamic 
form  of  the  closed-loop  system.  In  terms  of  pilot-vehicle  system  input/ 
output  characteristics  this  will  be  approximately, 

M(s)  •_  1 

I  ( s  )  (  s/ojc  +  1 ) 

As  the  pilot  urgency  or  aggressiveness  modifies  his  gain,  u'c  will  increase 
or  decrease,  with  the  dominant  closed- loop  system  time  constant,  1 /coc , 
waxing  and  waning  in  corresponding  fashion.  Thus  there  is  a  very  wide 
range  _>f  excellent  closed- loop  dynamic  response  properties  available  to 
the  pilot  which  is  easily  adjusted  in  direct  proportion  by  his  effort. 
In  the  words  of  Fig.  19  "pilot  compensation  is  not  a  factor  for  desired 
performance",  and  the  configuration  will  highly  rated.  For  the  other 
extreme  imagine  a  set  of  effective  airplane  characteristics  which  has 
dynamics  in  the  region  of  crossover  which  require  precise  adjustment  of 
the  pilot's  lead-lag  equalization  and  gain  to  close  the  loop  in  a  stable 
manner.  The  pilot  can  exert  closed-loop  control,  but  the  dynamic  quality 
and  even  closed-loop  system  stability  require  that  his  describing  function 
be  precisely  tuned  to  offset  the  controlled  element  deficiencies.  The 
pilot's  compensation  in  this  case  will  range  from  "considerable"  to 
"intense",  and  the  configuration  will  be  rated  very  poorly. 

The  "Stability  Margin  Gain  Tolerances"  factor  is  most  easily  des¬ 
cribed  when  the  p i lot - veh ic le  system  is  conditionally  stable.  In  this 
situation  the  system  becomes  unstable  if  the  gain  is  either  too  low  or  too 
high.  When  the  pilot  lead-lag  equalization  is  adjusted  to  maximize  this 
range  (which  will  ordinarily  provide  crossover  model-like  features  in  the 
nominal  crossover  region),  there  is  a  "total  available  gain  range"  (TAGR) 
through  which  the  pilot  can  maintain  some  semblance  of  closed-loop 
control.  Clearly,  the  more  narrow  this  range  becomes  the  more  difficult 
the  pilot's  adjustment  and  the  worse  his  rating  will  become. 
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The  "Neuromuscular  System  Coupling"  factor  can  become  important  when 
the  low-frequency  effective  airplane  dynamics  are  excellent  but  the 
closed- loop  system  gain  margin  in  the  region  of  the  neuromuscular  actua¬ 
tion  mode  is  reduced.  This  was  touched  on  briefly  in  Section  III,  and  is 
covered  in  Refs.  12  and  13.  The  resulting  closed- loop  system  instability 
is  high  frequency,  2-3  hz.  It  is  one  explanation  for  "roll  ratchet". 

The  "Attentional  Demands/Excess  Control  Capacity”  factor  is  primarily 
related  to  divided  attention  operations.  When  the  control  task  itself  is 
responsible  for  using  most  of  the  pilot's  excess  control  capacity  the 
reasons  for  this  are  invariably  due  to  factors  already  covered.  When  the 
managerial,  communication,  planning,  and  other  non-control  tasks  consume 
too  much  of  the  pilot's  available  attention  pilot  ratings  will  suffer. 
The  obverse  of  this  is  that  the  effective  vehicle  dynamics  must  very  good 
in  order  to  require  a  minimum  of  attention.  The  ratings  for  control  alone 
should,  in  general,  be  superior  and  the  unattended  operations  factors 
would  be  good  as  well. 

The  last  Table  6  factor  to  be  discussed  is  "Closed-Loop  Performance". 
Many  facets  of  task  performance  stem  directly  from  mission  requirements 
and  are  hence  mission-specific.  The  status  of  the  pilot -vehicle  system 
relative  to  mission  requirements  is  the  very  first  thing  the  pilot 
assesses  before  he  even  establishes  more  detailed  ratings.  Average  error 
performance  in  command  and  regulation  tasks  can  be  calculated  with  all  the 
pilot  models  once  these  inputs  are  defined.  These  estimates  can  serve  as 
one  basis  for  flying  qualities  assessments. 

There  are  other,  more  general,  closed- loop  dynamic  performance 
aspects  which  should  also  be  considered  in  flying  qualities  assessments. 
Four  of  these  are  listed  in  Table  7.  The  first  two  are  simple  statements 
of  closed-loop  dynamic  response  quality.  They,  in  essence,  suggest  that 
there  be  one  dominant  closed- loop  mode  which  can  have  a  damping  ratio 
greater  than  0.35  to  0.5.  (It  could  also  be  a  dominant  first-order  mode.) 
The  requirement  to  avoid  a  closed- loop  mid- frequency  droop  is  tantamount 
to  the  one  dominant  mode  prescription,  for  the  droop  will  show  up  as  an 
additional  minor  mode  with  a  longer  time  constant.  The  Neal-Smith 
criteria  (Ref.  78),  for  example,  call  specific  attention  to  the 
mid- .frequency  droop  a. id  require  that  it  be  less  than  3  dB  to  achieve 
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TABLE  7.  DESIRABLE  CLOSED  LOOP  DYNAMIC  FEATURES 


•  ADEQUATE  CLOSED-LOOP  DAMPING,  fCL  -  0.35-0.50 


•  AVOIDANCE  OF  CLOSED-LOOP  MID- FREQUENCY  DROOP 


•  MULTILOOP  CONTROL  VIA  SERIES  STRUCTURE 
FOR  SINGLE  CONTROL 


•  FREQUENCY  SEPARATION  OF  INNER,  OUTER  LOOPS, 

E .  G .  ,  wr  .  **  2 . 3  cjr  —  0.5-1. 0 
ci  ‘-o 


•  SIMPLE  CROSSFEEDS  TO  DIRECTLY  NEGATE 
SUBSIDIARY  RESPONSES 


•  CONTROL  HARMONY 


Level  1  ratings.  By  way  of  example  a  3  dB  mid- frequency  droop  can  be 
associated  with  the  presence  of  a  minor  mode  comprising  a  single  dipole 
pair  in  the  closed- loop  pilot -vehicle  system  (with  the  | zero/pole |  <  1.41) 
supplementing  the  major  dominant  mode. 

The  remaining  two  desirable  closed-loop  dynamic  features  are  asso¬ 
ciated  with  multi -loop,  single  control  axis  situations.  Common  examples 
of  this  include:  1)  the  control  of  altitude  wherein  altitude  error  is  the 
outer  loop  feedback  and  pitch  angle  is  an  inner  loop;  and  2)  hover  con¬ 
trol,  as  shown  in  Fig.  20.  "Desirable"  aspects  of  such  systems  include 
the  qualitative  feature  that  a  "series"  (rather  than  parallel)  closure  of 
the  outer  loop  is  possible  in  the  presence  of  an  inner  loop  system  which 
serves  both  independently  and  as  a  means  to  equalize  the  outer  loop.  Thus 
the  pilot  closure  of  a  pitch  attitude  loop  satisfies  an  attitude  control 
function  and  gives  rise  to  an  effective  outer,  altitude  control,  loop 
which  needs  very  little  if  any  further  pilot  equalization.  This  is  sup¬ 
ported  in  a  more  quantitative  sense  by  the  suggestion  for  the  separation 
of  crossover  frequencies  in  multiloop  systems  with  series  pilot  elements. 
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The  desirable  crossfeed  feature  listed  accounts  for  the  possibility 
of  pilot- induced  crossfeeds  to  reduce  or  eliminate  subsidiary  modes  or 
response  quantities.  A  common  example  of  this  is  an  aileron  to  rudder 
crossfeed  for  turn  coordination.  The  last  feature  on  the  list  of  desires 
is  control  harmony,  which  relates  to  multi-axis  control  conditions.  Force 
and  position  gradients,  pre-loads  and  centering  springs  and  other  manipu¬ 
lator  features  between  elevator  and  aileron  need  to  be  in  proper  balance 
so  that  the  effective  controlled  element  gains  in  each  axis  are  near 
optima,  interaxis  crosstalk  is  minimized,  etc.  Just  as  with  the  setting 
of  controlled  element  optimum  gain,  control  harmony  is  a  subject  of 
experimental  determination. 

As  a  consequence  of  flying  qualities  analyses  using  pilot-vehicle 
analysis  to  examine  the  factors  of  Tables  6  and  7,  the  analyst  can  develop 
a  set  of  conclusions  and  arrive  at  a  wide  variety  of  issues  and  possible 
problems.  Table  8  illustrates  the  type  of  problems  that  might  be 
uncovered  by  such  examinations  for  the  case  of  longitudinal  attitude  and 
path  control . 
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TABLE  8.  TYPICAL  PILOT  CENTERED  PATH  REGULATION  PROBLEMS 


ATTITUDE  CONTROL 

•  INADEQUATE  BANDWIDTH 

•  INNER-OUTER  LOOP  EQUALIZATION  CONFLICT 

•  LOW  STATIC  GAIN 

•  OVER-SENSITIVITY  TO  GAIN/ EQUALIZATION 


PATH  CONTROL 

•  PERFORMANCE  REVERSALS 

•  INADEQUATE  BANDWIDTH 

•  INADEQUATE  SEPARATION  OF  PATH  AND 
ATTITUDE  RESPONSES 

•  DIFFICULT  OR  CONFLICTING  CROSSFEEDS 

•  EXCESSIVE  DEPLETION  OF  SAFETY  MARGINS 

•  LOW  (HIGH)  EFFECTIVE  PATH  GAINS 
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Abstract 

In  missions  where  high  cognitive  and  managerial 
requirements  are  placed  on  the  piloc,  or  where 
failures  may  significantly  degrade  one  or  more  air¬ 
craft  control  axes,  the  pilot  must  divide  his 
attention  among  several  tasks.  The  pilot  and  sys¬ 
tem  behavior  in  such  divided  attention  conditions, 
and  the  combination  of  pilot  ratings  from  single  to 
multi-axis  conditions  are  treated  using  classical 
and  optimal  control  models  for  the  human  in  a  com¬ 
plementary  fashion.  Ic  is  shown  that  the  crossover 
frequency  and  closed-loop  system  performance  for  a 
given  axis  under  divided  attention  will  be  less 
than  full  attention  values  while  Che  remnant  and 
the  phase  margin  will  be  greacer,  and  chat  che 
model-based  crends  are  consistent  with  experiment. 
Also,  the  optimal  control  performance  index  used  in 
the  pilot  and  system  behavioral  modeling,  when 
"calibrated"  with  single-axis  correlations,  shows 
potential  for  che  development  of  subjective  rating 
estimates  for  multi-axis  casks. 

Introduction 

The  piloc  of  a  modern  high-performance  aircraft 
muse  perform  both  control  and  managerial  functions 
in  most  mission  phases.  In  some,  such  as  final 
approach  and  landing  or  a  variety  of  air-to-air  or 
air-to-ground  tracking  tasks,  the  control  function 
is  paramount  and  requires  taosc  of  the  piloc's 
available  attention.  In  che  past  the  more  diffi¬ 
cult  of  these  control  Casks  have  tended  Co  expose 
whatever  unfavorable  effective  aircraft  dynamic 
(flying  qualicy)  problems  were  present.  Conse¬ 
quently,  piloc  modeling  for  the  consideration  of 
critical  flying  qualities  has  tended  to  be  focused 
on  full  attention  or  nearly  full  acention  control 
operations.  Further,  the  vast  majority  of  piloc 
dynamics  and  pilot  racing  data  used  to  develop  fly¬ 
ing  quality  boundaries  have  been  obtained  for  situ¬ 
ations  where  che  dynamic  properties  of  one  control 
axis  are  varied  while  maintaining  the  dynamics  of 
ocher  axes  at  "good"  levels. 

With  modern  missions  in  which  high  cognitive 
and  managerial  requirements  are  placed  on  che 
piloc,  or  where  failures  may  significantly  degrade 
one  or  more  aircraft  control  axes,  Che  piloc  must 
divide  his  attention  among  several  axes  or,  more 
generally,  between  control  and  ocher  tasks.  In 
such  situations  the  vehicle  dynamics  in  each  axis 
being  controlled  by  the  piloc  must  be  superior  Co 
those  which  would  be  suitable  If  che  pilot  could 
devote  almost  full  attention  to  just  one  axis. 

This  state  of  affairs  brings  two  key  Issues  Co 
the  fore.  The  flrsc  is  pilot  and  system  behavior 
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In  divided  attention  or  task  interference  situa¬ 
tions,  where  several  tasks  compete  for  the  pilot's 
attention.  Critical  questions  In  these  situations 
relate  to  the  standards  that  should  be  sec  for  both 
desirable  and  minimum  flying  qualicy  levels  in  a 
given  axis  when  these  standards  are  conditioned  by 
the  ocher  control  and  managerial  tasks  at  hand. 
The  second  is  the  subjective  assessment,  or  racing, 
of  effective  vehicle  dynamics  in  multi-axis  control 
or  multiple  cask  conditions. 

This  paper  will  treat  both  of  these  issues  for 
multi-axis  control  casks.  It  begins  with  a  review 
of  piloc  modeling  for  divided  attention  conditions, 
developing  the  differences  to  be  expected  in  pilot 
and  pilot-vehicle  system  behavior  for  single-axis 
(full  attention)  and  multi-axis  control  tasks.  The 
major  data  source  for  which  pilot  ratings  for  both 
single  and  multi-axes  are  available  is  then  used  as 
a  basis  for  che  estimation  of  piloc  and  system 
behavior  data.  The  estimates  are  obtained  using  an 
opcimal  control  model  and  serve  to  augment  the 
piloc  racing  "data  set"  by  providing  estimates  for 
these  quantities  (which  were  not  available  from  Che 
original  experiments).  The  estimates  are  compared 
with  the  theoretical  expectations  and  with  the 
limited  experimental  evidence  available.  The  esti¬ 
mation  of  piloc  ratings  is  then  treated  using  an 
extension  of  previously  available  correlations. 
Important  feacures  in  these  developments  are  the 
complementary  connections  of  classical  theory  and 
opcimal  concrol  modeling  techniques. 

Divided  Attention  and  Task  Interference  Modeling 

The  classical  theory  of  pilot-vehicle  system 
dynamics  is  based  on  the  well-kncvn  "Crossover 
Model"  (see,  e.g.,  Ref.  1).  In  this  theory  the 
human  pilot's  behavior  is  characterized  in  mathe¬ 
matical  terms  as  a  random-input  describing  func¬ 
tion,  Yp,  plus  pilot-induced  "noise"  expressed  as  a 
power  spectral  density,  ^,n,  commonly  referred  to 
as  "remnant."  In  a  specified  task  with  the  effec¬ 
tive  vehicle  dynamics  described  by  the  transfer 
function  Y  ,  che  crossover  model  open-loop  describ¬ 
ing  function,  C,  takes  Che  form  shown  in  Fig.  1. 
This  form  is  established  by  appropriate  adjuscmenc 
of  Yp  by  the  piloc. 

In  the  crossover  model  there  are  three  key 
variables  —  the  crossover  frequency  ,  che  effec¬ 
tive  system  latency  r,  and  the  remnant  power  spec¬ 
tral  density,  *nn.  The  crossover  frequency  is 
sometimes  loosely  referred  to  as  che  piuot-vehlcle 
system  "bandwidth."  Ic  has  the  usual  feedback 
system  physical  interpretation  as  the  metric  which 
divides  the  world  of  the  control  system  into  two 
frequency  regimes:  below  i^.  the  benefits  of 
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feedback  are  present  (e.g.,  output  follows  input, 
error  is  reduced,  etc.)  while  above  the  system 
becomes  essentially  open-loop.  The  system  latency 
theoretically  appears  as  a  pure  time  delay, 
although  it  is  actually  a  low  frequency  approxi¬ 
mation  to  all  the  high  frquency  (i.e.,  well  above 
•q.)  net  lags  and  delays  in  the  system.  It  includes 
time  delays  and  neuromuscular  lags  contributed  by 
the  pilot  as  well  as  net  lags  from  manipulators  and 
other  manual  controller  elements,  higher  frequency 
effective  aircraft  dynamics,  etc.  The  effective 
aircraft  dynamics  include  aircrafc  plus  stability 
augmentation  plus  display,  etc.,  dynamics.  By 
using  the  effective  system  latency  as  a  normalizing 
var'able  che  crossover  model  dynamic  properties  can 
be  given  in  more  general  terms  as  either  a  normal¬ 
ized  crossover  frequency,  Tiq.,  or  a  phase  margin, 

-  V  2  -  rsj.. 

The  remnant  power  spectral  density  depends  on 
the  amount  of  pilot -gene rat ed  lead  required  to  make 
good  the  crossover  model  form,  the  nature  of  the 
manipulator,  and  other  features  of  the  system  and 
task.  It  arises  primarily  from  fluctuations  of 
attention  and  other  time-variations  internal  to  the 
pilot.  1  When,  as  is  usually  the  case,  these  time- 
variations  are  random  the  power  spectral  density  is 
continuous  and,  more  importantly,  the  remnant  spec¬ 
tral  density  will  scale  with  che  mean-square 
error.  1  11  This  is  depicted  in  Fig.  1  by  the  docted 
line  introducing  e?  inco  Che  remnant  block. 
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general  name  (“Perceptual  Scanning")  for  the 
divided  attention  or  task-interference  remnant  to 
account  for  the  divided  attention. 

Divided  attention  operations  can  be  considered 
as  an  extension  to  the  well-established  theory  of 
display  scanning  and  signal  sampling/ 
reconstruction.  1 »  3 »  S .  6  mosz  important  and 

fundamental  model  change  is  the  remnant.  A  theo¬ 
retical  model  for  remnant  due  to  quasi-random  cask 
swi tching/sampling  with  a  constant  average  dwell 
time  on  che  task  at  hand,  ,  and  a  sampling/ 
switching  interval  among  tasks  varying  about  some 
mean  value  Tg  was  derived  by  Clement^  and  validated 
for  visual  inputs.  6  This  remnant  form  is 


w<-> 

•  Tse2  (1  -  n)(l  -  6)  units2  \ 

f  /■“ff d  I2  1  rad/sec 

*  |i  +  (-y)  1  CO 

where  e2 

is  che  mean-square  of  Che  signal 
s  anpl  ed 

TS 

is  che  nean  sampling  incerval 

n 

is  effeccive  concrol  dwell  fraccion  * 

v*. 

Td 

is  effeccive  dwell  incerval  on  che  con- 
c  rol  cask  considered 

6 

is  normalized  lower  bound  on  che  douarn 
of  Ts:  T0/Ts  <  1 

This  power  spectrca  is  valid  at  frequencies  well 
above  the  low-pass  breakpoint,  2/7^,  defined  by  che 
control  dwell  time.  The  variability  of  the  switch¬ 
ing  interval  is  rot  purely  random,  but  is  subject 
to  a  lower  bound,  TQ.  This  is  represented  in  Eq.  1 
by  the  term  1  -  d.  The  numerical  values  in 
Table  l,  based  on  the  Ref.  6  experiments  involving 
the  visual  modalities  and  attention  switches 
between  displays,  provide  some  appr2ciation  for  the 
magnitudes  of  some  of  the  quanities  in  Eq.  \. 


TA2LE  1.  TYPICAL  VALUES  FOR  DIVIDED 
ATTENTION  REMNANT  QUANTITIES 


Figure  I.  Pilot-Vehicle  System  for  Divided 
Attention  Control  Tasks 

When  the  pilot's  full  attention  is  focused  on  a 
single-axis  control  task  the  crossover  frequency  of 
the  p  ilct -vehicle  system  is  maximized,  and  the  rem¬ 
nant  has  a  minimum  value  which  is  dependent  on  the 
details  (pilot-lead,  an ni pul  a  cor  characteristics, 
etc.)  af  the  particular  tasx.  When  more  chan  one 
axis  is  Co  be  controlled,  or  when  managerial  tasks 
are  present,  or  when  both  multi-axis  control  and 
managerial  operations  are  to  be  accomplished,  the 
pilot  must  snare  his  attention  among  the  several 
casks.  In  these  cases  che  crossover  model  will 
still  hold  for  the  control  activities,  ouu  with 
some  mo  di  f  ica  c  i  ons .  The  possible  nature  of  the 
changes  is  partially  exhibited  in  Fig.  I  by  the 
presence  oi  a  perceptual  describing  function,  Y^ , 
an  accompanying  additional  latency,  ,  and  a  more 


ITEM 

NOMINAL 

RANCH 

6  -  T0/Ts 

1/2 

r* 

O 

Miniauo  Tj  (sec) 

0.4  -  0.5 

Tj  for  Velocity 
Detection  (sec) 

O 

1 

O 

T^  for  Acceleration 

Detection  (sec) 

The  major  low  frequency  remnant  variations  are 
with  the  controL  dwell  fraction,  n,  and  the  mean- 
square  Input  to  the  pilot  (shown  as  mean-square 
error,  e?,  in  Fig.  1).  Perceptual  Scanning  remnant 
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is  reduced  by  increasing  the  dwell  fraction,  going 
to  zero  when  the  divided  attention  control  task 
becomes  full  attention.  (In  this  limiting  case 
full  attention  remnant  will  still  be  present;  this 
is  not  included  in  Eq.  1  because  Che  perceptual 
scanning  remnant  is  typically  ouch  larger  for  the 
divided  attention  conditions  of  interest  here.  It 
can  easily  be  included  because,  as  already  noted 
the  same  scaling  with  e?,  with  a  different  propor¬ 
tionality  factor,  applies  for  che  full-attention 
remnant.)  Just  as  for  full-accent  ion  remnant,  the 
increase  of  remnant  power  with  Che  magnitude  of  the 
pilot's  input  stems  from  time  variations,  in  this 
case  the  time  modulation  of  attention  Inherent  in 
the  perceptual  scanning  process. 

The  implications  of  divided  accention  have  been 
explored  theoretically  and  experimentally  for 
almosc  two  decades.  (The  classical  theory 

described  to  this  point  stems  from  Ref.  5,  and 

includes  Refs.  1,  3-8).  A  quantitative  example  of 
the  effects  of  divided  attention  on  task  perform¬ 
ance  is  presented  in  Fig.  2,  taken  from  Ref.  8. 

Here  the  forcing  function  is  white  noise  through  a 
third-order  Buccerworch  filcer  with  normalized 
breakpoint  t  •  0.25.  Full  attention  is  the 

lowest  curve.  The  divided  actencion  conditions  are 
then  shown  as  a  family  with  task  dwell  fraction,  n, 
as  the  parameter.  In  this  example  the  normalized 
concrol  cask  dwell  interval  is  sec  ac  Td/  t  *  1.5, 

and  the  divided  attention  remnant  normalized  lower 
bound  on  the  scanning  interval,  6  *  0.5. 
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Figure  2  illustrates  Che  pcofound  effects  of 
divided  attention  on  task  performance.  The  percep¬ 
tual  scanning  penalizes  the  performance  in  three 
important  ways: 

•  by  reduction  of  the  crossover  frequency 
(and  a  concomitant  increase  in  the  phase 
margin) 

o  by  increase  in  the  remnanc 

o  by  the  introduction  of  a  new  kind  of 
"stability"  constraint  —  "instability  in 
the  mean-square” 

While  all  of  these  ace  apparent  in  the  figure,  the 
third  point  requires  some  additional  explanation. 
The  full  attention  stability  limit  for  the  cross¬ 
over  model  function,  given  by 


is  approached  by  the  full  attention,  n  “  1  curve. 
In  general  the  mean-squared  error  for  the  feedback 
system  shown  in  Fig.  1  is  given  by, 


Figure  2.  effect  of  Divided  Attention  on  Task  Performance 
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where  e^  is  the  error  component  correlated  with  the 
forcing  function  and  e^  is  due  to  the  remnant. 
Because  the  remnant  component  is  proportional  to 
the  total  mean-squared  error, 


Ref.  9  to  reflect  task  interference  are  centered  on 
Che  increased  "measurement  noise”  contaminating  the 
subject's  observations  of  sensed  quantities. 

Specifically,  in  the  optimal-control  modeling 
approach,  10  the  human  is  considered  to  observe  an 
array  of  system  responses,  which  then  comprise  the 
subject's  observation  vector,  y  (t).  These  measure¬ 
ments  are  (each)  contaminated  with  (white)  observa¬ 
tion  noise,  Vy^,  or 


e^  +  F 


(«) 


so  that 


yP(°  ’  yPtrue(t)  +  vy 


(6) 


e 


2 


Finally,  the  noise-intensity  matrix  of  v  is 
assuned  diagonal,  with  elements  equal  to  Vj^,  where 
(5)  i  corresponds  to  the  lcf>  observed  variable  (e.g., 
displayed  error  in  one  axis). 


As  F  — »  1,  which  is  associated  with  the  limit  on  For  a  single-axis  control  task,  and  display,  V.. 
the  inequality  constraint,  F  «  e2/e2  <  1,  the  blow-  iS  taken  to  be  ”  11 
up  in  the  mean-square  occurs. 


The  performance  data  shown  in  Fig.  2  indicate 
that,  as  tasks  require  the  pilot's  attention  the 
dwell  fraction  for  the  primary  cask  must  decrease, 
and  che  optimum  operating  point  for  the  pilot 
shifts.  This  is  further  illustrated  in  Fig.  3, 
which  shows  an  increase  in  the  minimum  mean-squared 
error  and  in  phase  margin  and  a  decrease  in  the 
normalized  crossover  frequency  as  che  task  dwell 
fraction  decreases. 


Figure  3.  Effect  of  Divided  Accenci on  on  Phase 
Margin  and  Crossover  Frequency  for  Minimum 
eZ  and  F  «  0.5  and  l  (from  Ref.  8) 


Divided  attention  effects  can  also  be  modeled 
with  optimal  control  models  for  pilot  behavior. 
Reference  9  is  an  early  emp irical/analy c leal  scudy 
of  cask  interference  which  Cook  this  approach.  In 
this  study,  the  simultaneous  concrol  of  several 
independent  axes  was  the  situation  creacing  che 
Cask  Interference.  From  two  to  four  axes  were  con¬ 
trolled,  each  with  an  associated  display.  The 
results  were  compared  to  single-axis  control 
resuits.  A  oocel  for  task  interference  was  hypoth¬ 
esized,  and  this  model  also  focuses  on  remnant- 
related  parameters.  The  key  parameters  proposed  in 
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-  P„ 


‘Vi 


NCa,!’  Tt> 


(7) 


where  .  is  the  res  value  of  the  response  vari¬ 
able.  >rpi>  and  N(  cy  .T*)  is  the  random-input 
describing  function  of  a  dead  zone  of  half-width 
Tj.  This  dead-tone  model  represents  the  effect  of 
perception  or  indifference  thresholds  in  che 
subjects'  observations.  Finally  P0  is  Che  nominal, 
full-attention  noise-to-signal  ratio  for  the 
observed  variable,  y4.  This  (f ull-accencion) 
noise-to-signal  ratio  has  been  empirically  shown  to 
be  relatively  constanc  over  a  range  of  casks  and 
controlled  element  dynamics,  ar.d  is  usually 
measured  (for  foveal  viewing  in  the  case  of 
visually  sensed  responses)  to  be  approximately 
-20  dB  (i.e.,  .01  units  of  normalized  power  per 

rad/sec,  defined  over  positive  frequencies).1*'^ 
Note  that  all  the  above  is  consistent  with  the 
"standard"  optimal-control  modeling  approach 
described  in  Ref.  10,  for  example. 

Now  the  interference  model  of  Ref.  9  staces 
chat  che  effect  of  task  interference  is  reflected 
in  an  increase  in  the  effective  PQ,  the  noise  to 
signal  ratio  defining  Che  observation  noise  inten¬ 
sity.  Uith  che  description  of  cask  interference, 
the  noise  intensity  becomes 

‘Vi 

N(cy.,  Ti>  (8) 

where  is  che  fractional  attention  devoced  to 

cask  i  (or  more  precisely,  reponse  i),  and 
0  <  fj  <1.  Also,  since  the  human's  capacity  is 
finite,  ffi  ■  1.0  over  all  tasks  and  displays. 

This  model  was  validated  in  Ref.  9  by  comparing 

model-based  tracking  errors,  piloC  describing  func¬ 
tions,  and  remnant  (reflected  to  che  displayed 

error)  wich  experimental  results.  Because  the 
filter  dynamics  which  are  a  portion  of  the  pilot 
(model)'  dynamics  depend  on  observation  noises, 
variation  in  these  noises  manifest  themselves  in 
changes  in  pilot  (model)  dynamics  and  hence  in 

piloc-vehicle  closed-loop  dynamics.  (Note  chat  in 
Ref.  9,  the  threshold  effeccs  were  modeled  as  addi¬ 
tional  additive  observation  noise  that  did  not 
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scale  with  a y  .  The  threshold  effects  considered 
in  this  reference  were  primarily  due  to  parafoveal 
viewing.  The  use  of  the  dead-zone  describing  func¬ 
tion  did  not  appear  until  after  Ref.  9  was 
published  although  its  use  is  basically  an  exten¬ 
sion  of  results  presented  in  Ref.  9.) 

Attention  will  now  turn  to  the  application  of 
the  optimal-control  theoretic  model  for  task  inter¬ 
ference  to  the  modeling  of  other  experimental 
situations  involving  single  and  multi-axis  tasks. 
These  are  the  multi-axis  experiments  of  Dander 
(Ref.  11),  which  are  the  only  daca  set  currently 
available  in  which  pilot  ratings  are  available  for 
single  axis  and  multi-axes  made  up  of  the  single 
axis.  The  model  results  will  provide  an  "estimated 
pilot  behavior  and  system  performance  daca  base," 
which  is  otherwise  noc  available.  These  results 
can  then  be  compared  with  Che  trends  from  the 
divided  attention  theory,  and  also  with  some  of  the 
experimental  results  of  Ref.  6. 

The  Experimental  Data 

Dander11  has  reported  experimental  findings 
from  single-axis  as  well  as  multi-axis  tracking 
experiments.  The  only  experimental  results  were  in 
Che  form  of  subjective  ratings,  using  the  Cooper 
scale,  but  both  single-axis  and  multi-axis  ratings 
of  the  same  dynamics  were  obtained. 

Dynamically  independent  axes  were  considered  in 
Che  tracking  experiments.  Single-axis  as  well  as 
simultaneous  two-  and  three-axis  tracking  casks 
uere  performed,  and  the  casks  were  subjectively 
rated.  For  Che  two-  and  three-axis  experiments,  an 
overall  rating  was  given,  noc  a  rating  of  an  indi¬ 
vidual  axis  in  Che  cask.  Although  some  of  che  two- 
axis  cases  will  be  discussed  here,  the  ratings  for 
these  cases  were  considered  by  Dander  to  be  incon¬ 
sistent.  This  conclusion  was  reinforced  by  an 
independent  analysis.  12  Consequently,  the  two-axis 
ratings  will  noc  be  used  in  che  analyses  presented 
here. 

The  controlled-elesent  dynamics  in  each  axis 
were  varied  Co  obtain  a  wide  variety  of  combina¬ 
tions  of  configurations.  The  dynamics  for  each 
axis  uere  initially  selected  to  yield  Level  i,  2, 
and  Level  3  racings  in  the  single-axis  tracking 
task.  So  three  different  dynamic  elements  were 
selected  for  each  axis,  yielding  a  CotaL  of  nine 
different  dynamic  elemencs.  These  transfer  func¬ 
tions  are  listed  in  Table  2. 


TABLE  2.  DANDER  CO  STROLLED -ELEMENT  DYNAMICS 


Axis  Transfer 

Functions 

Good 

Dynamics  (H) 

Intermediate 
Dynanics  (M) 

Bad 

Dynanics  (L) 

a  *.(.04)(0.») 

l/T,  -  0 

l/T,  -  -.5 

l/T,  -  -.3 

(l/T»K5 .)(.?,.  k  4 ) 

«  0. 3(0.1) 

C*  -  .84 

C,  ■  -.3 

C$  -  --84 

fl  10.(0..) 

«Bp  ’  (J.O)ic'*  .51 

Cfl  -  -5 

ta  *  --2 

ta  '  --5 

The  configurations  to  be  considered  herein, 
along  with  the  racings  obtained  in  the  Ref.  11 
experiments  are  listed  in  Table  3,  with  configura¬ 
tion  symbols  used  in  later  graphical  presentations. 


TABLE  3.  CONFIGURATIONS  CONSIDERED 


Conf igurat ion 
Identlfyer 

No.  of 
Axes 

Subjective  Graphical 

Racings  Symbol 

e-H 

1 

2.5,2. 3, 3. ,3. ,3. 

O 

e-M 

1 

4.  ,4. 5, 4. 5, 5.  ,5. 

A 

8-L 

1 

6.  ,6.  ,6.5,7. ,8. 

y 

1 

2.5,3. ,3. ,4.  ,4. 

□ 

4rM 

1 

3.5,4. ,4. ,5,5. .5. 

o 

$-L 

1 

5. 0,5. 5, 6. 0,6. 5 

D> 

e-H 

1 

3. 0,3. 5, 3. 5, 3. 5, 3. 5 

<3 

Sp-HH 

2 

(2. ,3. .3. ,3. ,3.5)* 

a 

0$-KL 

2 

(6.  ,6.  ,6.5)* 

O' 

6jrML 

2 

(6. ,7. ,7.)* 

£ 

B-rLH 

2 

(4. ,4. ,4.)* 

v' 

B^e-KHH 

3 

5.5,6. ,6.  ,6. 

• 

6*3-HMH 

3 

5.  ,5.5,6. 

▲ 

e?9-HLH 

3 

7.  ,9. 

▼ 

epe-MHH 

3 

5. 5, 6. 5, 6. 5 

■ 

epe-MMH 

3 

7.5 

♦ 

epe-MLH 

3 

8.  ,8.5 

► 

Bsd-LHH 

3 

9. ,10. 

< 

8$3-LMH 

3 

9.0 

m 

epe-LLH 

3 

8. 5, 9. 5 

w 

a 

Data  considered 

inconsistenc-not  specifically 

used 

The  command  signals  in  each  axis  consisted  of 
sums  of  sinusoids,  generating  random-appearing 
signals.  The  amplitudes  and  frequencies  of  the 
sinusoids  used  in  each  of  the  three  command  signals 
are  listed  in  Table  4  below,  along  with  the  approx¬ 
imate  maximum  command-signal  amplitude.  Note  chat 


TABLE  4.  DANDER  COMMAND  SICNAL  CONTENT 


Axis 

Frequency 

Unsealed 

(r/s) 

Amplitude 

0C 

0.367 

4.80 

0.483 

5.50 

(^max 

•  ±0.75  in. ) 

0.816 

1.32 

1.44 

1.64 

*c 

0.367 

4.80 

0.483 

2.20 

(  ^max 

*  ±45  deg. ) 

0.816 

1.15 

1.32 

1.64 

Sc 

0.483 

0.92 

0.816 

0.33 

(^max 

••±1.8  units) 

1.32 

0.33 

I 


« 
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the  units  of  displayed  commands  are  inches  of 
command  bar  displacement  on  the  display  (8C),  angle 
of  rotation  of  the  command  bar  (^),  and  "display 
units"  (approximately  7  degs  of  rotation  of  a 
needle  on  a  circular  dial  •  1  unit)  (£^).  Also  it 
is  noted  that  most  of  the  power  in  Che  signals  is 
in  frequencies  below  0.5  rad/sec.,  with  the  maximum 
frequency  in  the  signals  equal  to  1.32  rad/sec. 

Optimal-Control  Modeling  of  Che  Experiment 

The  modeling  process  new  involves  representing 
the  task  in  terms  of  Che  parameters  in  the  optimal 
control  framework.  10  These  include  the  following: 

1.  An  appropriate  objective  function  represen¬ 
tative  of  the  task  (or  a  generic  task 
cricical  to  the  performance  of  the  actual 
cask). 

2.  The  value  for  the  neuromotor  time  constant 
in  each  axis. 

3.  The  parameters  that  determine  the  observa¬ 
tion  noise  intensities: 

a.  Nominal  noise-Co-signal  ratio  (PQ) 

b.  Perception/indifference  threshold  (T^) 

c.  Atcentional  allocation  (may  be 
optimized). 

4.  Vector  of  human  operator  observations  in 
the  cask. 

5.  Observation  time  delay  (routinely  selected 
as  r  »  0.2  sec). 

6.  Neuromotor  noise-to-s ignal  ratio, 

7.  Characterization  of  the  command  signals. 


Finally,  the  interprecation  of  g^  requires  some 
discussion.  In  the  OCM,  the  selection  of  g^ 
defines  the  frequency  range  over  which  the  open- 
loop  system  amplitude  ratio  approximates  a  K/s-like 
form.  In  connection  with  the  OCM  it  is  often  cited 
that  g^  is  selected  to  yield  a  desired  neuromotor 
time  constant,  t,,,  in  the  pilots’s  describing  func¬ 
tion  obtained  from  the  model.  But,  when  the  total 
pilot  describing  function,  Y  ,  is  actually  con¬ 
structed  from  its  various  elements  in  the  OCM, 
Che  tjj  established  in  this  fashion  is  canceled  by  a 
directly  compensating  lead,  leaving  the  actual 
estimated  Y  with  no  (\,  s  t  l)-1  lag.  Scill  it 
has  been  convenient  Co  adjust  g^  in  this  fashion 
even  though  the  lag  will  later  disappear.  In  this 
vein,  che  value  of  the  desired  neuromotor  time 
constant  used  is  either  0.1  sec,  0£_  the  tn  that 
yields  the  lowest  error  (e.g.,  best  performance), 
whichever  is  greater.  Notice  that  after  t,,  is 
determined  in  the  above  fashion,  this  "operating 
point"  is  associated  with  some  weight  g^  in 
Jaxis  f  This  value  may  also  infer  the  subject’s 
subjective  trade  between  performance  (<0  and  work¬ 
load  (o'j).  And  since  piloc  lead  and  oj  are  corre¬ 
lated,  this  procedure  maximizes  the  possibility  of 
relating  the  resulting  value  of  Jtaslc  to  the  sub¬ 
jective  racing  of  the  task.  More  about  this  later. 

The  parameters  related  to  observation  noise 
were  selected  as  follows.  For  the  experimental 
set-up  considered,  and  using  visual  perception 
thresholds  of  0.05  deg  and  0.1  deg/sec  of  visual 
arc  and  arc  rate  at  the  pilot’s  eye,  the  dead-zone 
widchs  for  the  perception  thresholds  are 

8  —  0.03  in;  0.05  in/sec 

^  —  1.5  deg;  3.0  deg/sec 

E  —  0.14  unit;  0.28  units/sec 


For  these  experiments,  the  appropriate  objec 
cive  function  for  each  axis  was  chosen  to  be 


•^axis^ 


(9) 


And  for  a  multi-axis  Cask,  the  objective  function 
for  che  task  was 

Naxes 

Jtask  ”  S  ^axis; 

i  (10) 


The  nominal  full-ac:encion  foveal-viewing  value  of 
-20  dB  was  selected  as  the  basic  noise-co-signal 
ratio,  PQ,  for  both  perceived  error  and  error  rate 
in  each  axis.  And  finally,  for  multi-axis  tasks, 
the  attention  was  optimized  among  the  observations 
for  each  axis,  with  che  fractional  attention  f^ 
(see  Eq.  8)  constrained  to  be  equal  for  both  error 
and  error  rate  in  each  axis.  This  optimization  was 
performed  using  che  criceria  that  f^'s  were  to  be 
chosen  such  thac  Jtas^  was  minimized,  subject  to 
che  constraint  that 


I  £i-i 

axis  (11) 


The  justification  of  this  selection  involves  three 
considerations.  The  first  relaces  to  the  selection 
of  equal  (unity)  weighting  on  each  Jax£s  ^  in  che 
definition  of  dtas^  In  multi-axis  tasks.  This 
decision  was  based  on  the  instructions  given  to  che 
subjects  in  Che  experiment.  They  were  to  attempt 
to  minimize  the  errors  in  all  controlled  axes. 
Thac  is,  they  were  instructed  thac  no  axis  was  to 
be  given  preference,  which  would  then  define 
primary  and  secondary  sub-tasks. 


For  further  discussion  of  the  actentional  alloca¬ 
tion,  see  Ref.  13. 

The  observations  assumed  available  to  che  sub¬ 
jects  were  error  and  error  races  in  each  axis.  The 
displays  were  either  integrated  or  considered 
sufficiently  close  together  such  thac  additional 
observation  errors  due  to  scanning  or  parafoveal 
viewing  was  neglected. 


Secondly,  the  normalization  of  the  mean-square 
error  wich  the  mean-square  command  deals  nicely 
with  che  fate  that  different  units  and  different 
command-signal  itrengchs  were  user  in  the  axes. 
This  will  be  discussed  furcher  witr.  regard  to  sub¬ 
jective  racings. 


The  neuromotor  nolse-co-signal  ratio  was  estab¬ 
lished  as  -20  dB.  In  addition  to  neuromocor  noise 
reflecting  a  source  of  remnant,  the  use  of  this 
noise  has  also  been  found  to  be  necessary  for  the 
OCM  to  yield  gain  and  phase  margins  appropriate  for 
human  operator  modeling.  For  the  relation  between 
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this  noise  and  stability  margins  of  LQG  control¬ 
lers,  the  reader  is  referred  to  Ref.  14. 

Finally,  the  command  signal  for  each  axis  was 
modeled  as  a  second-order  Markov  process,  using  Che 
f ilcer 


And  the  pilot  phase  compensation  is  then 
*p(  v^)  -  arg  Yp(j  Sfc) 


The  phase  margin  is  determined  from 


(14) 


Fc(s) 


s2  +  2(.7)(.5)s  +  (.5)2 


(12) 


^  "  "  -  »r2  YpV-i  “fc> 

while  effective  pilot  latency  is 


(15) 


Note  that  the  break  frequency  of  0.5  rad/s  corre¬ 
sponds  to  maximum  command  power  below  0.5  rad/sec. 
The  gain  Ac  was  selected  such  that  the  maximum 
command  amplitude  listed  in  Table  4  was  twice  che 
standard  deviation  (2c)  of  the  random  command 
generated  with  che  above  filter. 

Results  of  the  Model  Analysis 


\  *  (  s/2  -  4^)/  ujc 


(16) 


Finally,  closed-loop  bandwidth  is  defined  as 


“S', 


“J  closed-loop 
phase  -  -90° 


(17) 


The  above  OCM  set-up  data  were  then  inserted 
into  a  version  of  the  PIREPS  computer  program*^  to 
obtain  OCM-based  estimates  for  the  pilot  dynamic 
behavior  and  various  syscem  performance  metrics. 
The  model  results  obtained  include  both  cime-  and 
frequency-domain  quantities.  The  time-dooan 

results  include  root-mean-square  tracking  error, 
c£,  control  manipulation  defleccion  race,  a^,  the 
objective  function  magnitude,  J  and  the  opti¬ 

mized  accencional  allocation,  f^,  for  each  axis 
(for  multi-axis  tasks).  The  frequency-domain 

results  include  pilot  describing  function 

6i/ci(j  a),  and  remnant,  ^n(w),  along  with  the 
concrolled-elemenc  Cransfer  function  in  Bode  form, 
c^/  5i  (jw).  From  this  information,  the  crossover 
frequency,  uCi,  pilot  phase  compensation  at  cross¬ 
over,  $p(uc),  phase  margin,  4^,  effective  pilot 
latency,  Tr  ,  and  closed-loop  bandwidth  may  be 
determined.  These  parameters  are  defined  graphi¬ 
cally  on  some  example  results  in  the  Appendix. 
This  phase  of  che  analysis  is  similar  to  that  of 
Refs.  15  and  16.  The  pilot's  describing  function 
(in  each  axis)  corresponds  to  the  fora 


V-> 


e)  “,\>bs 
3  “  ’-n  +  1 


Vj“) 


(13) 


The  key  resulcs  are  summarized  in  Tables  5,  6, 
and  7.  Table  5  liscs  the  single-axis  resulcs. 
When  the  single  axis  estimates  are  converted  to 
classical  forms  the  results  show  trends  of  racing 
with  increased  effective  lacency,  te ,  and  reduced 
crossover  frequency,  ,  which  are  consistent  with 
similar  sets  of  vehicle  dynamics  considered  in 
Refs.  1  and  17.  Tables  6  and  7  summarize  the  two- 
and  three-axis  results,  respectively.  These 

results  are  discussed  further  below. 

Since  one  of  the  fundamental  effects  of  task 
interference  is  Increased  remnant,  this  area  is  the 
first  to  be  discussed.  Shown  in  Fig.  4  are  the 
modeling  resulcs  for  the  remnant  ^n,  normalized 
with  mean-squared  error  for  two  cases  considered. 
This  figure  does  indeed  indicate  an  increase  in 
remnant  at  low  frequencies  for  the  two  axis  task 
over  the  single  axis.  Note  chat  this  tvo-axis  task 
involves  poor  dynamics  laterally,  with  the  same 
good  dynamics  in  pitch  as  che  single  axis  case.  As 
a  result  of  the  interference  from  difficult  lateral 
dynamics,  the  attention  allocated  to  the  pitch  axis 
reduces  from  100  percent  (t-itch  s  1.0)  to  20  per- 
cent  <£pitch  ‘  °-2)- 


I 


\ 


TABLE  5.  SINGLE  AXIS  SUMMARY 


case  * 

u.(r/s ) 

“3W 

*1 <  Te  > 

ac/ci 

J 

PCR 

3.2 

3.7 

37”(  .289) 

Uu  .2° 

.188 

.103 

0.46 

2.5-3 

3.1 

3.5 

28. 4(  .346) 

42.6 

.295 

.193 

.124 

4-5 

®L 

3.1 

3.5 

25. 5( .  363) 

44.8 

.413 

.273 

.245 

b-i 

3.1 

3.6 

37. 2(  .  297) 

50.7 

.202 

.119 

.055 

2.5-4 

<t>. 

3.1 

3.5 

29 . 5( . 3  35) 

64.5 

.281 

.176 

.110 

3.5-5 

♦l 

2.8 

3.1 

2  1. 5(  .427) 

64.0 

.523 

.383 

.422 

5-6.5 

(  v12) 

2.9 

3.6 

43. 4( .280) 

38.9 

.23- 

.125 

.071 

3-3.5 

*  Tn  ’  0.1s, 

excep  c 

as  noted 
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TABLE  6.  TV.'0-AXIS  SUMMARY 


case/axis* 

“c 

“Sw 

°t/  °i 

a‘i/ Si 

Actn 

Jtask 

POR 

8$,HH/  8 

2.75 

3.3 

41.51  .308) 

39.0 

.230 

.111 

.5 

.140 

(2.9) 

n 

2.7 

3.3 

42. 2( .308) 

43.7 

.242 

.126 

.5 

8b, HL/  8 

2.15 

2. 85 

49.21.331) 

33.9 

.315 

.127 

.2 

.646 

(6.2) 

/♦ 

2.75 

2.95 

22. 1( .431) 

62.8 

.594 

.421 

.8 

8b, LH/  8 

3.0 

3.4 

26.31.371) 

44.1 

.472 

.235 

.8 

.430 

(4.0) 

h 

2.2 

2.87 

51.01.309) 

34.1 

.306 

.134 

.2 

*  xn  in  each 

axis  as  in  Table  5 

TABLE  7.  THREE-AXIS  SUMMARY 


» 


case/a xi s** 

"c 

“3W 

bn(  ) 

4?c 

aj  0i 

Ob/gi 

a  t  c  n 

J  t  a  s  V. 

POR 

0?3,HH:  8 

2.5 

3.  1 

45.  1  .314) 

36.6 

.  262 

.118 

.33 

.  300 

5.5-5 

0  5  3 , HHK / b 

2.45 

3.  1 

46.1.313) 

39.  1 

.273 

.130 

.33 

8  b  S ,  K  H  H  /  S 

2.0 

2.8 

54.21.312) 

J  7 ,  2 

.330 

.126 

.33 

0  b  3,  KLH  /  0 

1.9 

2.6 

53.11. 339) 

32.  5 

.368 

.  134 

.  1  5 

.935 

7. -9.0 

8bS,HLH/ b 

2.65 

2.  9 

23.51.438) 

61.3 

.644 

.445 

.67 

0  b  3 , HLH / 3 

1.75 

2.6 

60.71.292) 

16.9 

.389 

.140 

.  IS 

0  b  3 ,  L  H  H  /  0 

2.9 

3.3 

26.81. 380) 

42.8 

.  505 

.244 

.57 

.539 

9.0-10.0 

0  b  3 ,  L  H  H  /  b 

1.  9 

2.6 

56 . 5( . 308) 

27  .  9 

.360 

.140 

.  1  5 

0  b  3 , LHH / B 

1.63 

2.6 

60.21.319) 

4 . 4 

.  394 

.132 

.18 

0  b  3 ,  M  M  H  /  0 

2.6 

3.0 

32.31.395) 

36.0 

.430 

.165 

.40 

.631 

7.5 

0  b  S , MMH / b 

2.7 

3.  1 

34.51.366) 

58  .  5 

.  4  20 

.186 

.40 

0  b  3 ,  M  M  H  /  3 

1  .  7 

2.5 

59.51.318) 

7.3 

.380 

.131 

.20 

0  b  3 ,  M  H  H  /  0 

it 

* 

* 

* 

.370 

.211 

.48 

.426 

5.5-6. 5 

8  b  3 ,  .M  H  H  /  b 

* 

* 

* 

* 

.293 

.13  1 

.27 

0  b  3,  MHH  /  3 

* 

★ 

It 

* 

.349 

.  138 

.25 

8  b  3 , HMH / 0 

* 

* 

* 

* 

.295 

.  126 

.25 

.42  1 

V_A 

o 

1 

O' 

o 

0b3,HMH/b 

* 

* 

it 

it 

.381 

.210 

.46 

0b3, HMH/ 3 

* 

* 

it 

it 

.333 

.137 

.29 

0  b  3 , LMH / 8 

* 

* 

* 

* 

.533 

.322 

.54 

.  859 

9.0 

0  b 8 , LMH / b 

* 

* 

* 

* 

.490 

.242 

.28 

SbS.LMH/ 3 

★ 

* 

* 

* 

.388 

.140 

.18 

0b3, MLH / 0 

* 

* 

* 

it 

.508 

.257 

.25 

1.23 

in 

CO 

i 

O 

CO 

0  b  3 ,  M  L  H  /  b 

* 

* 

* 

it 

.695 

.468 

.59 

0bB,MLH/3 

* 

* 

* 

it 

.4  28 

.142 

.  1  3 

0b8, L  LH / 8 

* 

* 

it 

it 

.682 

.386 

.39 

1.74 

8. 5-9. 5 

8b3,LLH/b 

* 

* 

it 

a 

.788 

.5  12 

.  5  1 

0  b  3 , LLH / 6 

* 

* 

* 

* 

.472 

.143 

.10 

-  no:  determined 
Tn  in  each  axis  as  in  Table  5 
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All  Remnant  Reflected  to  Error 

o-yyj  *-20d8  Each  Obs  (  plus  threshold ) 
o-yuj  *-20d8  Eoch  Input 
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Figure  4.  Model-Based  Error  Remnant  Speccrua 
Task-Interference  Effect 


One  reason  the  controlled  elements  of  Fig.  4 
were  chosen  for  single  axis,  dual  axis  comparisons 
is  the  possibility  to  coapare  the  estimation 
results  with  similar  data.  Reference  6  considered 
a  cwo-axis  case  with  one  axis  involving  an  unstable 
controlled  element  similar  to  the  lateral  dynamics 
for  the  two-axis  case  modeled  here.  Figure  5  shows 
chat  the  remnant  shown  for  the  pitch  axis  in  the 
two-axis  task  in  Fig.  4  compares  favorably  with 
thac  from  Ref.  6  for  the  nacural-scanning  case.  To 
the  extent  chat  these  remnant  data  from  Ref.  6  are 
a  reasonable  surrogate  for  the  Dander  case  this 
close  correspondence  is  compelling.  However,  the 
key  result  involving  remnant  is  thac  an  increase  in 
remnant  is  Che  model-based  result  of  task  inter¬ 
ference,  as  expected  from  the  theory  reviewed 
earlier. 


Figure  5.  Comparison  of  Model-Estimated  and 
Ref.  6  Normalized  Remnant  Spectra 


Another  effect  of  interference  noted  was  up. 
regression.  Shown  in  Fig.  6  are  some  resulcs  from 
Ref.  6  along  wich  some  model-based  results  gener¬ 
ated  In  this  study.  The  ratio  of  crossover  fre¬ 
quency  to  chat  for  a  single-axis  task  is  plotted 
versus  effective  dwell  fraction,  ne ,  defined  as 


Figure  6.  Crossover  Regression  with  Task 
Interference  (after  Ref.  6) 


where  II  *  0.566  fron  Ref.  6,  and  n*  actual  dwell 

fraction.  If  dwell  fraction  is  taken  to  approxi¬ 

mate  fractional  attention,  l±,  as  defined  herein 
(or  vice  versa),  one  can  compare  model-based 

results  to  those  of  Ref.  6.  The  modeling  results 

are  la  general  agreement  with  the  experimental  data 
and  exhibit  the  same  trends  as  expected  from  the 
theory.  Therefore,  we  can  say  that  the  cask  inter¬ 
ference  effect  of  regress  ion  is  also  observed  in 
the  modeling  resulcs. 

Another  comparison  becween  the  theory  and  Che 
modeling  results  discussed  herein  is  revealed  in 
rig.  7.  The  results  in  this  figure  are  consistent 
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Figure  7.  Trend  of  Optimal  Operating  Poincs 
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wich  the  parametric  results  shown  in  Figs.  2  and 
3.  The  minima  in  Fig.  2  represent  the  optiaum  per¬ 
formance  at  for  a  given  dwell  fraction,  n.  Again 
interpreting  n  to  be  analogous  to  fractional  atten¬ 
tion,  f ^ ,  the  model  based  resulcs  for  the  Dander 
configuration  in  Fig.  7  corresponding  to  "best" 
dynamics,  for  example,  are  the  model-based  resulcs 
that  show  trends  comparable  to  the  minimum  mean- 
squared  error  points  of  Figs.  2  and  3.  Llsced 
alongside  the  data  points  in  Fig.  7  is  the  frac¬ 
tional  attention  used  for  control  of  the  "best" 
axis  dynamics  ( 8H)  as  well  as  some  results  for  the 
"worst"  axis  dynamics  (9L).  Noce  that  as  frac¬ 
tional  attention  (dwell  fraction)  decreases  for  a 


given  controlled  element,  (optimum)  tracking  per¬ 
formance  is  degraded,  ng  regresses,  and  phase  mar¬ 
gin  increases  —  all  generally  consistent  with  the 
divided  attention  theory. 


Pilot  Rating  Estimation 


Given  the  consistencies  between  full-  and 
divided-attention  behavior  described  above  the  next 
concern  is  the  estimation  of  pilot  workload  asso¬ 
ciated  with  the  multi-axis  control  tasks.  This  is 
quantified  in  the  Dander  data  set  by  Cooper  pilot 
ratings.  It  is  well  known  (e.g.,  Refs.  1,  12,  17) 
thac  the  subjective  ratings  of  flying  qualifies 
provided  by  pilots  depend  on  task  performance  and 
che  pilot  equalization  needed  to  offset  any  delete¬ 
rious  vehicle  dynamic  characteristics.  In  tasks 
with  a  significant  closed-loop  tracking  concent 
these  features  can  be  reasonably  represented  by 
normalized  tracking  error  and  pilot  lead  equaliza¬ 
tion. 


Pilot  lead  equalization  will  skew  che  pilot's 
outpuc  ("concrol  input"  in  the  OCM)  spectrum  toward 
higher  frequency  contenc.  Consequently,  le3d 
equalization  when  reflected  into  a  component  of  an 
Integral  performance  index  will  be  manifest  as 
increased  concrol  input  rate.  Also,  the  Cask  per¬ 
formance  in  terms  of  mean-squared  error  is  a 
natural  component  of  a  system  performance  metric  or 
index.  With  these  two  facts  in  mind,  Hess10  sug¬ 
gested  thac  che  combination  of  error  "-'red  and 
piiot-oucpuc-(concrol  input)-race-squared  might  be 
used  as  a  mecric  giving  an  indication  of  che  sub¬ 
jective  racing.  With  a  particular  weighting  (g)  on 
the  control  rate  term  chis  is  just  the  optimal  con¬ 
trol  performance  index  in  the  OCM.  Using  the 
limited  experimental  data  available  Hess10  was  able 
to  show  chat  a  reasonable  linear  correlation  did 
indeed  exist  between  che  OCM  performance  index  J 
and  the  pilot  racing,  and  the  basis  for  this  corre¬ 
lation  was  extended  further  in  Ref.  19.  Indeed, 
one  reason  chat  che  presenc  study  was  undertaken 
was  the  hope  that  Chis  relationship  could  be 
expanded  to  che  multi-axis  situation. 


To  determine  if  the  (appropriately  chosen) 
objective  function  did  in  fact  correlate  with 
multi-axis  as  well  as  single-axis  ratings,  the 
Dander  resulcs  are  again  considered.  In  the  task 
objective  functions  chosen  in  this  modeling  work, 
che  error  was  normalized  with  the  variance  of  the 
command  signal  in  the  respective  axis.  Not  only 
does  this  offset  the  effects  of  different  units  in 
che  axes,  but  It  also  cancels  che  effect  of  command 
signal  strengen  on  the  racings.  17  Finally,  wich 
the  model-based  determination  o-  .  che  weighc  on 
control  race  in  che  objective  function,  one  may 
cc-sider  che  resultant  weight  as  "normalizing" 


control  rate.  These  normalized  quantities,  mean- 
square  error  and  concrol  race,  are  tabulated  in 
Tables  5,  6,  and  7,  along  with 

Shown  in  Fig.  8  is  Che  correlation  between 
J[3Sk»  3S  D°dele:l>  and  Che  subjective  ratings  of 
Che  task.  The  correlation  between  Jtasjc  and  POR 
from  the  single-axis  results  appears  Co  hold  for 
the  multi-axis  results  as  well.  This  result  would 
seem  to  indicate  Chat  the  ratings  reflect  the 
actual  performance  and  workload  (stick  rate)  in  the 
overall  task,  or  thac  the  subject  does  not  accept 
Che  fact  thac  the  performance  should  be  worse  since 
che  task  is  more  difficult.  In  other  words,  he 
does  noc  accept  worse  performance  and  workload  In  a 
multi-axis  cask  compared  to  single-axis  results. 


-20  -1.0  0.0 
Lc<3  Jtosk 


Figure  8.  "Cost  Function"/Racing  Correlation 


The  results  also  tend  co  support  che  hypothesis 
thac  determining  the  weightings  ga  in  che  manner 
discussed  leads  to  the  "correct”  relative  weight¬ 
ings  or.  control  rate  in  che  axis,  and  the  relative 
weight  between  control  rate  and  normalized  error. 

Conclusions 


This 

study  has 

illustrated 

chat : 

l) 

When  the 

piloc  must 

divide 

across  several  control  axes  che  crossover  frequency 
will  be  less  than  full  actencion  values  while  Che 
remnant,  closed-loop  system  performance  (error)  and 
phase  margin  will  be  larger. 

2)  Pilot  modeling  estimates  using  the  optimal 
control  model  adjusted  to  take  into  account  divided 
attention  or  cask  interference  among  several  con¬ 
trol  axes  exhibit  characteristics  consistent  with 
those  listed  above. 

3)  Comparison  of  the  estimates  and  available 
data  show  similar  trends. 

4)  The  objective  function  (optimal  concrol 
performance  Index)  used  In  the  pilot  and  system 
behavioral  modeling,  when  "calibrated"  with  single¬ 
axis  correlations,  show  potential  for  the  develop¬ 
ment  of  subjective  racings  estimates  for  multi-axis 
casks. 
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Appendix 

Included  in  Figs.  A. 1  and  A. 2  are  the  open-  and 
closed-loop  Bode  plots,  and  the  pilot  describing 
function  Y  (ju),  respectively.  Shown  on  these 
figures  are  the  frequency -domain  measures  deter¬ 
mined  from  the  modeling  process.  The  measures, 
specifically,  are  gain  crossover  frequency, 
phase  margin,  4^,  closed-loop  bandwidth,  BW,  and 
pilot  phase  compensation,  (including  the  phase 

due  to  Tq-js  at  i^.  and  (1  +  j^u^.)-^). 
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